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Abstract 

How can one prove that a given state is entangled? In this paper we review different methods that have been proposed 
for entanglement detection. We first explain the basic elements of entanglement theory for two or more particles and 
then entanglement verification procedures such as Bell inequalities, entanglement witnesses, the determination of 
nonlinear properties of a quantum state via measurements on several copies, and spin squeezing inequalities. An 
emphasis is given to the theory and application of entanglement witnesses. We also discuss several experiments, 
where some of the presented methods have been implemented. 
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1. Introduction 

Entanglement was first described by Einstein, Podolsky, and Rosen [l| and Schrodinger Q as a strange 
phenomenon of quantum mechanics, questioning the completeness of the theory. Later, Bell recognized that 
entanglement leads to experimentally testable deviations of quantum mechanics from classical physics Q. 
Finally, with the advent of quantum information theory, entanglement was recognized as a resource, enabling 
tasks like quantum cryptography [4( , quantum teleportation [5| or measurement based quantum computation 
. Together with the rapid experimental progress on quantum control, this lead to a rapidly growing interest 
in entanglement theory and many experiments nowadays aim at the generation of entanglement. 

Indeed, in the last years an enormous progress on the generation of entanglement has been achieved. For 
instance, six or eight ions have been entangled @, H| , photons have been used to demonstrate entanglement 
between six particles or ten qubits 9, 10] and in diamond nuclear and electronic spins have been entangled 



These systems allow for individual addressing of the parties. In other types of systems only collective 
measurements are possible. Through spin squeezing, entanglement of 10 7 atoms were created in cold atomic 
clouds and large scale entangling operations were realized in optical lattices of 10 5 two-state atoms (l3j . 
As the underlying techniques of quantum control improve continuously, it can be expected that in the near 
future even larger systems can be entangled. 

In any of these experiments, typical questions arise: How can one be sure that entanglement was indeed 
produced? How can one detect the presence of entanglement? Can we quantify the entanglement in the 
experiment? These questions are difficult to answer and many possible ways to tackle this problem have 
been proposed. These range from Bell inequalities, entanglement witnesses and spin squeezing inequalities to 
entropic inequalities, the measurement of nonlinear properties of the quantum state and the approximation 
of positive maps. 
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In general, one could imagine several desirable properties for an entanglement verification procedure. 
First, a crucial requirement is that the scheme under consideration must be easy to implement. Depending 
on the type of experiment, some measurements are easy to carry out, and some are not. In experiments using 
photons, the count rates are often not very high, making it impossible to characterize the experimental state 
completely. Then, one has to find the measurements that allow to conclude much about the entanglement 
content of the given state. In addition, the scheme should be robust against noise, and it should also detect 
weakly entangled states. 

Second, most experiments nowadays aim at the generation of entanglement between more than two parti- 
cles and for this case different entanglement classes exist. Therefore, an entanglement detection scheme has 
to be capable of distinguishing between these different classes. In an TV-qubit experiment it is not sufficient 
to say that some qubits are entangled, one has to prove that all of the iV qubits are entangled. 

Third, the conclusion that the experimentally generated state was entangled should not depend on some 
assumptions about the state. Most experiments aim at the generation of a special state, and clearly one can 
use this to design appropriate measurements. But the final conclusion that the state was entangled, must 
not rely on assumptions concerning a special form of the state or its purity. 

In this paper, we review different methods for the characterization of entanglement in experiments. This 
includes the mere detection of entanglement, but also its experimental quantification, the estimation of the 
state fidelity and the characterization whether a prepared state is useful for some task or not. We also explain 
several experiments, where these methods have been used. Due to their higher experimental relevance, we 
especially explain methods that can be used for the characterization of entanglement between more than 
two particles. A special emphasis is given on the method of entanglement witnesses, as they are the most 
frequently used tool nowadays, moreover, they have also applications in other areas of quantum information 
theory. 

Throughout this paper we restrict our attention to the case of discrete systems; mainly we consider the case 
of qubits. It should be stressed that entanglement in continuous variable systems like harmonic oscillators 
or light modes is also under intensive research. However, as the underlying theory of Gaussian states is 
significantly different from the case of qubits, we have not addressed these systems here. Moreover, for the 
entanglement theory of continuous variable systems already several good reviews exist 1J, [l5j, [l6|, [l7| . 

As entanglement is a central topic in quantum information theory, there are some other excellent re- 
view articles that are related to some questions addressed here. For the case of bipartite entanglement, a 
recent exhaustive review was written by the Horodecki family [lit ] and entanglement measures have been 
reviewed in detail by Virmani and Plenio (l9| . Some older proposals for entanglement detection schemes 
have been summarized by Terhal [20| and a survey about conceptual differences and possible problems 
of entanglement verification schemes has been given by van Enk, Liitkenhaus and Kimble [2lll . Moreover, 
several other overview articles concerning bipartite and multipartite entanglement [22j, 

HESSE!!!], 

en- 
tanglement measures [53, 0] or Bell inequalities [1^, [13, EH have been published. Concerning the experi- 
mental technique s, p hoton experiments have been reviewed in Refs. [32, and ion trap experiments in 

Refs. (aESlHlIIIai- 

This review is organized as follows: The Sections 2-4 give an overview about theoretical aspects of 
entanglement, while the Sections 5-8 explain different entanglement detection methods in detail. 

In Section 2, we give an introduction into entanglement theory of bipartite systems. We introduce the 
notions of entanglement and separability, explain several separability criteria and the phenomenon of bound 
entanglement. We also introduce entanglement witnesses as a theoretical concept and explain how to con- 
struct and optimize them. In Section 3, we discuss entanglement between more than two particles. We explain 
the different classes of multipartite entanglement and introduce the most relevant families of multiparticle 
entangled states, such as GHZ states, cluster states or Dicke states. In Section 4, we introduce entanglement 
measures. We explain their basic properties and give the definitions of the most popular measures. 

In Section 5, we review Bell inequalities as the oldest tool to verify entanglement. We explain how Bell 
inequalities bound the set of correlations originating from local hidden variable models and how the violation 
of a Bell inequality implies the quantum phenomenon of entanglement. We also explain parametric down- 
conversion as a process to generate entangled photons, leading to a violation of Bell inequalities and also an 
experiment, where this process has been used to demonstrate violation of a Bell inequality under locality 
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conditions. 



In Section 6, we review entanglement witnesses as a tool for the detection of entanglement. We show how 
entanglement witnesses can be implemented with simple local measurements and that they are especially 
useful for the detection of multipartite entanglement. We also explain several experiments, where witnesses 
have been used in various experimental situations. Finally, we explain how entanglement witnesses can be 
useful in other parts of quantum information theory, such as quantum cryptography. 

In Section 7, we review other entanglement detection schemes that have been implemented or proposed. 
This includes entropic inequalities, estimates of the concurrence and positive maps via measurements on 
several copies, variance based criteria, and nonlinear entanglement witnesses. 

Finally, in Section 8 we give an overview over entanglement detection with collective observables. This 
concerns spin squeezing inequalities and optical lattices, but also entanglement verification schemes in spin 
models, where global observables such as the magnetic susceptibility of the energy can be used to derive 
statements about the entanglement in the system. 

2. Bipartite entanglement 

In this Section, we will explain the basic notions of bipartite entanglement. As many experiments nowa- 
days aim at the generation of multiparticlc entangled states, we will explain the theory of multipartite 
entanglement in the following section. However, the study of bipartite entangled states will already enable 
us to introduce the central concepts of entanglement detection. 

2.1. Entanglement of pure states 

Let us assume that we are given two quantum systems. The first one is owned by one physicist, called 
Alice, and the second one by another one, called Bob. The physical states of Alice's system may be described 
by states in a Hilbert space Ha of dimension d,A, and in Bob's system in a Hilbert space Kb of dimension 
ds- The composite system of both parties is then described by vectors in the tensor-product of the two 
spaces H = Ha ®Hb- Thus, any vector in Ha ®Hb can be written as 



with a complex oIa x ds matrix C = (ctj). To keep the notation simple, we often write tensor products of 
vectors as \a) (g> \b) = \a)\b) = \ab). Now one can define separability and entanglement for pure states. 

Definition 1 (Entanglement for pure states). A pure state \ip) £ H is called a product state or 
separable if we can find states \(j) A ) £ Ha and \<p B ) £ Hb such that 



holds. Otherwise the state \ip) is called entangled. 

Physically, the definition of product states means that the state is uncorrelated. A product state can 
thus easily be prepared in a local way: Alice produces the state \(j) A ) and Bob produces independently the 
state \4> B ). If Alice measures any observable A and Bob measures B, then the probabilities of the different 
outcomes factorize. Thus, the measurement outcomes for Alice do not depend on the outcomes on Bob's 
side. 

Before proceeding to the definition of entanglement for mixed states we shall mention a very useful tool 
in the description of entanglement for bipartite systems. This is the so-called Schmidt decomposition. 

Lemma 2 (Schmidt decomposition). Let = Y^i^Li c ij\ a i^j) *= Ha®Hb be a vector in the tensor 
product of two Hilbert spaces. Then there exists an orthonormal basis |a,) of Ha and an orthonormal basis 
\/3j) oIHb such that 



AaAb 




(1) 



\f) = \4> A ) ® \<t> B ) 



(2) 



R 




(3) 
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holds, with positive real coefficients Afc. The Xk are uniquely determined as the square roots of the eigenvalues 
of the matrix CC\ where C = (c^ ) is the matrix formed by the coefficients in Eq. (1). The number R < 
mm{dA,d,B} is called the Schmidt rank of If the are pairwise different, then also the \c*k) and 
are unique up to a phase. 

The proof of this statement can be found in many textbooks 4(| 41, 42, Note that pure product 



states correspond to states of Schmidt rank one. Therefore, the Schmidt decomposition can be used to 
decide whether a given pure state is entangled or separable. 

2.2. Entanglement of mixed states 

In a more general situation one does not know the exact state of a quantum system. It is only known 
that it is, with some probabilities Pi, in one of some states \4>i) G Ti. This situation is described by a density 
matrix 



= y^Pil&K&l; with = 1 and > 0. (4) 



In a given basis this density matrix or state is represented by a complex matrix. This matrix is positive 
semidefinitj^l and hermitian, since all the operators \<pi){4>i\ are positive and hermitian. Due to the condition 
on the pi also Tr(g) = 1 holds. Conversely, any positive semidcfmite matrix of trace one can be interpreted 
as a density matrix of some state. This leads to a geometrical picture of the set of all states as a convex set. 
This means that given two states Q\ and g2, their convex combination g = agi + (1 — ct)g2 with a £ [0; 1] 
is again a state. This property holds also for combinations of more than two states. Given some pi > with 
YliPi = 1 then the convex combination YliPiSi °f some states is again a state. We call coefficients pi > 
with the property X^Pi = 1 often convex weights. 



Now one can define separability and entanglement for mixed states according to Ref. [44j . The idea is the 
same as for the pure state case: A state is separable, if it can be produced locally, otherwise it is entangled. 

Definition 3 (Entanglement for mixed states). Let g be a density matrix for a composite system. 
We say that g is a product state if there exist states g A for Alice and g B for Bob such that 

g=g A ®g B . (5) 

The state is called separable, if there are convex weights pi and product states gf £g> gf such that 

Q = Y.P^t®Qf (6) 

i 

holds. Otherwise the state is called entangled. 

Physically, this definition discriminates between three scenarios. First, a product state is an uncorrelated 
state, where Alice and Bob own each a separate state. For non-product states there are two different kinds 
of correlations. Separable states are classically correlated. This means that for the production of a separable 
state only local operations and classical communication (LOCC) are necessary. Alice and Bob can, by 
classical communication, share a random number generator that produces the outcomes i with probabilities 
Pi. For each of the outcomes, they can agree to produce the state gf ® gf locally. By this procedure they 
produce the state g — J2iPi@t ® Qf ■ This procedure is not specific for quantum theory, which justifies the 
notion of classical correlations. Otherwise, if a state is entangled, the correlations cannot originate from 
the classical procedure described above. In this sense entangled states are a typical feature of quantum 
mechanics. 

For our later discussion it is very important to note that the set of separable states is a convex set. This 
is clear from the definition of separability, obviously a convex combination of two separable states is again 
separable. Furthermore, the definition of separability implies that any separable state can be written as a 
convex combination of pure product states. Thus, the set of separable states is the so called convex hull of 



1 A hermitian matrix M is called positive semidefinite (M > 0), iff the eigenvalues are non negative. Equivalently, M > if 
{ip\M\ip) > for all \ip). 
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Fig. 1. (a) Schematic picture of the set of all states as a convex set and the set of separable states as a convex subset, (b) 
Different schematic picture of the same sets. Here, it is stressed that the extremal points of the separable states (the pure 
product states), are also extremal points of the set of all states, hence they are located on the border of the total set. 

the pure product states. Further, any separable state can be written as a convex combination of maximally 
d\d 2 B pure product states. This follows from Caratheodory's theorem, which states that any element of a 
rf-dimensional convex set can be written as a convex combination of d+ 1 extremal points of this set [451 ]. For 
the special case of two qubits, however, this bound can be improved and any separable state can be written 



as a convex combination of four product states only [46|, |47[ ■ The set of separable states is also shown in 
Fig. 10 

Given the definition of entanglement and separability, it is very natural to ask whether a given density ma- 
trix is separable or entangled. This is the so-called separability problem. There are several criteria known that 
imply separability or entanglement of a state. However, up to now, no general solution for the separability 
problem is known. 

2.3. Separability criteria 

In this Section, we will present some criteria for bipartite entanglement. As it is not our aim to discuss 
all of them in detail, we focus our discussion on the most important ones. 

2.3.1. The PPT criterion 

Let us start with the criterion of the partial transposition. In order to formulate this, note that we can 
expand any density matrix of a composite quantum system in a chosen product basis as 

N M 

8 = EE QijM\i)(j\®\k)(l\- (7) 

i,j k,l 

Given this decomposition, one defines the partial transposition of g as the transposition with respect to 
one subsystem. Thus, there are two partial transpositions: The partial transposition with respect to Alice 
is given by 

N M 

T "= EE ^«K>0'l®l fc K'l (8) 

i,j k,l 

and similarly we can define q Tb by exchanging k and I instead of the i and j. Note that the partial trans- 
position is related to the usual transposition by g T — (q Ta ) Tb and thus q Tb = (q Ta ) t ■ 

It is worth mentioning that the partial transposition depends on the product basis in which it is performed. 
But one can show that its spectrum does not depend on the basisO] We say a density matrix g has a positive 



2 A quantitative study on the shape of these sets in space where the coordinates are the density matrix elements is given in 
Ref. Ul. 



3 Note that this also holds for the full transposition. 
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partial transpose (or: the matrix is PPT) if its partial transposition has no negative eigenvalues, i.e., it is 
positive semidefinite: 

q Ta > <=> q Tb > 0. (9) 

If a matrix is not PPT, we call it NPT. Now we can formulate the PPT criterion (also called Peres-Horodecki 
criterion), originally introduced in Ref. [49l |. 

Theorem 4 (PPT Criterion). Let g be a bipartite separable state. Then g is PPT. 

Proof. This fact follows directly from the definition of separability in Eq. (6) since for a separable g = 
EkPkQk ^el 3 ™ have g TA =J2 kPk (g£) T ®gZ = J2 kPk g£ ®gi>0. □ 

This theorem provides a very strong criterion for the detection of entanglement. For a given density matrix 
one can easily calculate the partial transpose and compute its spectrum. If one finds negative eigenvalues 
one can conclude that the state is entangled. Given this result, the question arises if this criterion is also 
sufficient for separability, i.e., whether g TA > implies separability. As it was shown already shortly after 
the discovery of the PPT criterion, this is the case only in low dimensional systems: 

Theorem 5 (Horodecki Theorem). If g is a state ina2x2or2x3 system, then g TA > implies that 
g is separable. In other dimensions this is not the case. 

Proof. For the proof of the 2 x 2 or 2 x 3 case, see Ref. [Ho|. For the first counterexample in a 2 x 4- or 
3 x 3-system see Ref. [5l| . We will give an example in the context of bound entanglement later in Section 
2.4. □ 

Although the PPT criterion does not constitute a necessary and sufficient criterion, it is the most popular 
criterion. This is, beside its simplicity, due to several reasons. First, the fact that it provides a complete 
characterization of entanglement for two-qubit systems makes it very appealing, since two-qubit systems 
are the most studied bipartite systems. Second, it has been shown that the amount of violation of the PPT 
condition can be used to quantify entanglement [H, (see also Section 4.2.4). As we will see in Section 
3.2.2, however, the PPT criterion is of limited use for the investigation of multipartite entanglement. 



2.3.2. The CCNR criterion 

Since the PPT criterion does not detect all states, the question arises, how one can prove that a state is 
entangled, if the PPT criterion fails. For this problem, many criteria have been proposed. Here, we want to 
explain the computable cross norm or realignment criterion (CCNR criterion) [54U55j|. since it is simple and 
strong. 

In order to formulate the CCNR criterion, we utilize the Schmidt decomposition in operator space. For a 
density matrix g, such a decomposition is 

g = J2^G£®G k 3 , (10) 
fc 

where the A& > and G k and G k are orthonormal bases of the observable spaces of Ha and H b ■ Such a 
basis consists of d 2 hermitian observables which have to fulfill Tr(G k Gf) — Tr{G k Gf) = Ski - Such observ- 
ables are often called local orthogonal observables (LOOs) [56]. For instance, for qubits the (appropriately 
normalized) Pauli matrices together with the identity form a set of LOOs. 

Eq. (10) is nothing but the Schmidt decomposition from Sec. 2.1, but applied to the space of observables F1 
Therefore, the A& are (up to a permutation) unique and if the A& are pairwise different, the G k and G k are 
also unique (up to a sign). The A& can be computed as in the Schmidt decomposition: First, one decomposes 
g = fikiG^^iGf with arbitrary LOOs G k and Gf, then, by performing the singular value decomposition 
of Hki one arrives at Eq. (10), the A& are the roots of the eigenvalues of the matrix fi^ . One can formulate: 

Theorem 6 (CCNR criterion). If the state g is separable, then the sum of all A& in Eq. (10) is smaller 
than one: 

Hence, if Xk > 1 the state must be entangled. 

4 Note that one can also define a Schmidt rank for density matrices [H?! . this is, however, not related to Eq. (10). 
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Proof. For a pure product state, g = \a)(a\ (£> \b)(b\ the Schmidt decomposition is directly given, and the 
statement is trivial. For the general case g — X)fePfcl a fe)( a fcl ® \bk)(bk\, note that the sum over the Schmidt 
coefficients defines a norm on the set of density matrices. Then we have due to the triangle inequality 
ll^ll — SfcPfelll a fc)( a fe| ® |&fc)(fefe||| = 1- Note that the question, whether other functions of the Schmidt 
coefficients can be used for entanglement detection, has been studied in Ref. 

Let us discuss this criterion. First, the two different names come from the fact that this criterion has been 
discovered in two different forms. It has been first found while investigating cross norms of density matrices 
[H^ |. the above formulation was taken from (H3 |. Independently, it has been formulated in a way that one 
realigns the entries of a density matrix and calculates then the trace norm of this new matrix [Hit ]. This 
second formulation is somehow more complicated, it allows , however , to see formal similarities with the PPT 
criterion [(3(3] or to generalize it into different directions 61, 62, 63|. 



The remarkable fact is that the CCNR criterion allows to prove the entanglement for many states where 
the PPT criterion fails. Combined with its simplicity, this makes it a useful tool for the analysis of entan- 
glement. However, it does not detect all entangled states of two-qubits [64[. Therefore, one may view it as 
complementary to the PPT criterion. 

2.3.3. Other criteria 

Beside the PPT and CCNR criterion, there are many other approaches to derive separability criteria. Not 
all of them will be needed for the future discussion. However, in order to give the reader a coherent overview 
we would like to give a small list: 

(i) The rang e criterion was one of the first criteria for the detection of states for which the PPT criterion 
fails 51]. It states that if a state g is separable, then there is a set of product vectors |aj&j) such 
that the set {\aibi}} spans the range of g as well as the set {|q*6,)]T^1 spans the range of g TA . This 
condition allowed to discover many states which are entangled, but not detected by the PPT criterion 



51|, |65|, |66j . However, it cannot be used if some state is affected by noise: Then, the density matrix 
and its partial transpose will usually have full rank, hence the condition in the range criterion is 
automatically fulfilled. 

(ii) The PPT criterion is an example of a criterion using a positive, but not completely positive map. These 
objects are defined as follows. Let Hb and He be Hilbert spaces and let B(Hi) denote the linear 
operators on it. A linear map A : B(Hb) — > B(Hc) is called positive if it maps hermitian operators 
onto hermitian operators, fulfilling A(A^) = A(A)^ and it preserves the positivity, i.e., if X > then 
A(X) > 0. Note that the second condition implies that it maps valid density matrices onto density 
matrices, up to normalization. A positive map A is called completely positive when for an arbitrary 
Hilbert space Ha the map Ia ® A is positive, otherwise, A is positive, but not completely positive. 
Here, I4 denotes the identity on B(Ha)- For example, the transposition map T is positive, but not 
completely positive: while X > implies X T > 0, the partial transposition does not necessarily 
preserve the positivity of the state. 

From this, other entanglement criteria similar to the PPT criterion can be formulated from other 
positive, but not completely positive maps. For any separable state g and any positive map A we have 

(]U<aA)(0)>O. (12) 

Furthermore, it has been shown in Ref. [50( that a state g is separable if and only if for all positive 
maps A the relation in Eq. (12) holds. In this sense, the separability problem is equivalent to the 
classification of all positive maps. Of course, in order to develop a separability criterion, in Eq. (12) 
only the positive, but not completely positive maps arc of interest. 

This problem was considered in the mathematical literature for a longer time jf}?],!!!, 69|, 7(J, however, 



it has not been solved yet. From the perspective of quantum information theory, the classification of 
positive maps was under intensive research and has led to many new positive, but not completely 



positive maps, resulting in strong separability tests [71|, l72|, |73j, [7 



The symbol \a*) denotes the vector resulting when all coefficients of \a) in a certain basis are complex conjugated. Note that 



\a*)(a*\ = \a)(a\T. 
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Another example of a positive, but not completely positive map is the reduction map [76j |. This is 
on one system defined as 

A R (X) = Tr(X) ■ 1 - X. (13) 

Consequently a separable state has to fulfill \a®^ R {q) — QA®^- — Q > 0. This is the reduction criterion 
for separability. In Ref. [76| . it has been shown that the reduction map is a so-called decomposable 
map, that is, it can be written as 

K R = P l +P 2 o T, (14) 

where Pi,P2 are completely positive maps, and T is the transposition. This implies, that ILa <S) A fl 
can never detect entanglement, unless 1a ® T detects itOO Therefore, the reduction criterion is weaker 
than the PPT criterion. Nevertheless, the reduction criterion is interesting for several reasons: For 
two qubits, it also provides a necessary and sufficient criterion for separability 76], and it is further 
closely connected to the process of distillation (see Section 2.4) and certain entanglement witnesses 



(sec Section 7.2). Furthermore, one can extend the reduction map in the following way [72|, |77[: If U 
is a unitary matrix with U T = —U, then the map A(X) = Tr(X) ■ 1 — X — UX T W is a positive, but 
not completely positive map. Such a unitary transformation can only exist in even dimensions. But 
then, the map clearly improves the reduction criterion and can detect states that are not detected by 
the PPT criterion. 

For our purposes, the theory of positive maps will become important for the following reason. We 
will see that entanglement witnesses are good tools for the experimental detection of entanglement. 
Any entanglement witness, however, gives rise to a positive, but not completely positive map via 
the Choi-Jamiolkowski isomorphism, as will be explained in Section 2.5.3. Furthermore, the theory of 
positive maps will later enable us to construct nonlinear witnesses and to derive further entanglement 
tests (see Sections 7.3.2 and 7.4). 
(iii) The majorization criterion relates the eigenvalues of the total state with reduced states [78j ]. For a 
general state g one takes qa — Ttb(q) as the reduced state with respect to Alice and denotes by 
V = (pi,p2f-) the decreasingly ordered eigenvalues of g and by Q = (gi, 92, ■■■) the decreasingly 
ordered eigenvalues of gA- The majorization criterion states that if g is separable, then 
fc fc 

X><X> (15) 

8=1 t=l 

holds for all k. The same inequality holds, when gA is replaced by the reduced density matrix of the 
second system gs — TrA{g)- Physically, this criterion means that for a separable state g is more 
disordered than qa, and this criterion results also in inequalities for the entropies of g and gA Z^, 8C 



mini]. Jt 

was shown in Ref. (83|, however, that the majorization criterion follows from the reduction 
criterion hence any state which can be detected by the majorization criterion, can also be detected 
by the reduction criterion and, consequently, by the PPT criterion. Similarly, one can derive criteria 
for entropies from other positive, but not completely positive maps, leading to a generalization of the 
majorization criterion (84. l85l. l86j] . 
(iv) There are further algorithmic approaches to the separability problem. These methods formulate separa- 
bility in terms of (convex) optimization problems or semidefinite programs, in order to derive numerical 
algorithms for separability testing. 



Let us explain the method of the symmetric extensions [87ll88l]. since this is one of the most powerful 
separability criteria. The idea of this algorithm comes from the following observation. If a state g = 
^2kPk\ a k)(a-k\ <8 \bk)(bk\, is separable, we can find a symmetric extension to three parties (A, B and 
A\) as 



g (AB Al) = y pk \ ak ){ ak \ g, \ bk ) {bk \ g, | 0fc )(o fc |. (16) 



^(ABAi) = 



6 In fact, all positive maps with dimHs = 2 and dimWg < 3 are decomposable 1681. 17011 . from which it follows, that the PPT 
criterion is necessary and sufficient for low dimensions. 
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The state g( ABAl *> has now the following properties: (a) it is PPT with respect to each partition, (b) 
the reduced state on the first two parties is given by g and (c) it is symmetric under the exchange 
of the parties A,A\. Similarly, one can define for a separable symmetric extensions to more parties 
A,B,A u ...,A k . 

The key observation in Ref. [871 ] is that the question "Is there a symmetric extension with the 
properties (a) - (c)?" can be directly formulated as a feasibility problem in semidefinite programming. 
Semidefinite programs are a family of optimization problems that can not only be solved efficiently, but 
under weak conditions the global optimality of the found solution can also be proven (see Ref. [89( for 
a review) . Practically this means that the above mentioned feasibility problem can be directly tackled 
with standard numerical packages, and if no extension is found, the algorithm can also prove that no 
extension exists, consequently the state must be entangled. 

In fact, the method of symmetric extensions delivers a hierarchy of separability criteria, the criterion 
looking for the extension to k + 1 parties is always stronger than the criterion with the extension to 
k parties. It was then shown in Ref. (see also Ref. 88]) that this hierarchy is complete in the 
sense that any entangled state is detected in some step of the hierarchy. In addition, due to symmetry 
requirements, the number of parameters in the semidefinite program increases only polynomially in 
the number of extensions. Still, however, the numerical effort is considerable: while the first step of 
the hierarchy is, for system up to 5 x 5 feasible on a standard computeiQ, the higher order steps are 
difficult to perform. 

Beside the above mentioned algorithms, there are several other algorithmic approaches. First, the 
method of symmetric extensions was extended to the multipartite case [9lj. Then, in Ref. [921 ] it 
was shown that on can view the separability problem as an optimization problem with polynomial 
constraints, which can, in principle, be tackled with a hierarchy of semidefinite programs 
algorithms for separability tes ting have been proposed by making use of convex geometry 



97, 98 



93 
94 



Other 
96 



99L lOCj . see also Ref. 101 [ for an overview. In addition, it should be noted that semidefinite 
programs can also be used to tackle othe r pro blems in quantum infor mat ion theory, such as the 
existence of local hidd en v ariable models [1021 ] . distillability of states [1031 ] and the estimation of 
localizablc information 110411. 



Finally, these algorithms pose the question concerning the algorithmic complexity of the separability 
problem, that is, the question how much time is needed to decide whether a given state is entangled 
or not. It has been sho wn th at the separability problem, considered as a weak membership problem, 
is an NP hard problem [l05j . That is, deciding separability with an accuracy of e in an d x d system 
requ i res a n effort exponentially in d, if the inverse accuracy 1/e should increase exponentially in d. 



101 



105( | . As recently sho wn, t his result also holds if the if the inverse accuracy 1 je is required to 
increase only polynomially [106j . It should be noted, however, that this statement does not exclude 
the possibility of a simple solution for any fixed dimension. In addition, as this statement concerns 
only the computational complexity, it does not exclude possible physical insights into the separability 
problem. 

(v) Anothe r approach to the separability problem makes use of criteria based on covariance matrices 



(cms) film [Toi . [iM llldf . Here, one considers for some observables M k the covariance matrix 7 of 



these observables in some state, given by its entries 

7ij = (MiMj) - (Mi)(Mj). 



(17) 



Then, one tries to formulate separability criteria from this representation. The motivation for this 
approach comes from the theory of infinite-dimensional systems. There, an important family of states 
consists of Gaussian states, which are described by such a CM, using the canonical conjugate variables 
X and P (for reviews on this subject see Ref s. 14, 15, 1|| 17|). Moreover, for these Gaussian states 



the separability problem has been solved 111 . 112 . 113| 



With the help of these covariance matrices one can formulate a general covariance matrix criterion 



(CMC) [I09l . lll4| . For this, one chooses the observables in Eq. (17) to be the set {Gf <g> 1, l®Gf } with 



Numerical code for this first step is available at http://www.iqi.caltech.edU/documeiits/spedalieri/pptsetestl.m. 
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local orthogonal observables Gf and Gf. Then, 7 has a 2 x 2 block structure, and if g is separable, 



one can find matrices ka/b = 12kPkl(\' t Pk 
on the subsystems, such that 



A/B 



A 


C 




K A 






> 


C T 


B 




k b 



)) which are convex combinations of CMs of pure states 



(18) 



On the one hand, this criterion is very strong, as is detects many bound entangled states, and, with 
the help of filtering operations (see Sections 2.5.3 and 6.8), it can detect all entangled states for two 
qubits, just as the PPT criterion[f] Moreover, the CCNR criterion can be shown to be a corollary of 
the CMC. For instance, from the CMC it follows that for separable states 



||<7||l < ^[l-Tr(g A )][l-Tr{el% 



(19) 



where ||C||i is the trace norm of the mat rix Ca = (Gf (8 Gf) — (Gf)(Gf), which is a systematic 
improvement of the CCNR criterion [58L Ill4l Ill5| . In addition, the CMC can be related to other 
known separability cr iteria (lTi Evfy . 

It is not clear, how to use the CMC for entanglement dete ction in the multipartite case. The closely 
connected criterion of the local uncertainty relations (LURs) [l07| is, however, experimentally relevant 
and will be discussed in Section 7.3.1. 
(vi) A similar approach uses the expec t ation valu e matrix (EVM) of a quantum state to investigate its 
entanglement properties 118 . Ill9l ll20L 121 1 . For a set of operators {Mfe} one defines the positive 
semidcfinite EVM \ as 

(M}M 3 ). 



Xij 



(20) 



Typically, one chooses the observables to be of the type {M^} = {A^ ® B{\. Then it is straightforward 
to see that for a separable state g = J2kPk0k ® 0k one nas X(f?) = J^kP^xief) ® x(2fc)> hence 
x(g) is also separable and its separability can be investigated by criteria such as the PPT criterion. 
Especially, if Bj = B\ one has x Tb = x(q Tb )- This has two interesting consequences: First, a violation 
of the PPT criterion for x may be indicated by the n egat i vity of some subdeterminants of x, leading 
to nonlinear or variance- like criteria for separability 118l 122 1 . Second, if the PPT criterion for \ is 
violated, there is a vector \x) = (xi) such that (x\x Tb \x) < and hence a positive operator P = X^X 
with X = Y,i%i M i such that Tr(gP TB ) < while Tr(gP TB ) = (x\x Tb \x) > for all PPT states. 
The observable P Tb will later be called an entanglement witness (sec Section 2.5) and therefore the 
EVM allows to investigate for which states a violation of t he PPT criterion can be proved with a 
given, restricted set of observables {M&} = {Ak ® B[} [l2C| . Furthermore, the EVM can be used to 



investigate entanglement and separability for qubit-mode systems [123l Il24| . 
(vii) A different line of developing separability criteria uses linear contractions and permutations [60(. The 
main idea behind this approach is the following: given a density matrix g — Y]^ k i Qijki\i){j\ <8> \k)(l\ 
expanded in a product basis, all information concerning the state is encoded in the coefficients gijku 
which build a (oIa x oIa) x [cLb x ds) tensor. For pure states, it is normalized, ||ftjfci||i = 1, where 
|| V||i = Tr(VXXr) denotes the tr ace norm, and QijH is considered as an (d A ) x (d%) matrix. It 
is easy to see that for a product state, any permutation Tt(ijkl) of the indices i,j,k,l yields also a 
tensor with || g„(ijki) \\ 1 = 1, where, depending on the permutation, Q v Ujki) might be considered as 
a (dyids) x (d^ds) matrix. Since the trace norm obeys the triangle inequality, this means that for 



any separable state ^ I2 a P<*\\ 



(q) 



1, where g = Y^ a P a Q ^ turns out, that in 



the end for two parties only two permutations yield independent separability criteria, resulting in the 
conditions 



|f?(ij'Zfc)l|l 



< 1 and 



\Q(ikjl) 111 



< 1. 



(21) 



For general systems, however, there are NPT states which are not detected by the CMC and filtering 
entangled PPT states which not detected by CMC. 



moreover, there are 
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It can be directly seen that the first condition is nothing but the PPT criterion, while the second one 
corresponds to the CCNR criterion. Thus, this approach delivers a way to formulate these two criteria 
in a unified way . Further, it can be extended to the multipartite case, where many new separability 
criteria arise 62, 63|, see also Section 3.5. 



Other approaches to detect entanglement use special observables, they are called either Bell inequalities 
or entanglement witnesses. Bell inequalities will be discussed in Section 5 and entanglement witnesses in 
Section 2.5. Before doing that, we would like to discuss the phenomenon of bound entanglement. 

2.4. Bound entanglement 

For many tasks in quantum information theory, like teleportation or cryptography one ideally needs 
maximally entangled two-qubit states, i.e., singlet states. In the real world, however, noise is unavoidable, 
thus only mixed states are available. So one has to deal with the question, if and how one can create a singlet 
state out of some given mixed state. This leads to the problem of the so-called distillation of entanglement. 
It can be posed as follows: Assume that there is an arbitrary, but finite, number of copies of an entangled 
quantum state g distributed between Alice and Bob. Can they perform local operations on the states, assisted 
by classical communication, such that at the end they share a singlet state, 

...® g L ^ c kr) = -J=(|oi) - | io))? (22) 



k copies 

If this is the case, the state g is called distillable, otherwise we call it undistillable or bound entangled. 

It is not clear from the beginning that one can distill quantum states at all. However, the first distillation 
protocols were derived in Ref. [125J, [126J, showing that it is in principle possible. The question whether a 
given state is distillable is, of course, difficult to decide, since there is no restriction to a specific kind of local 
operations, as well as to the number of copies. The failure of a special protocol is not enough to conclude 
that a state is undistillable. However, some simple criteria are known: 

Theorem 7. If a bipartite state is PPT, then the state is undistillable. If a state violates the reduction 
criterion (e.g., due to a violation of the majorization criterion) the n the state is distillable. 

Proof. For the sufficient condition to be undistillable, see Ref. [l27[. For the sufficient condition to be 
distillable, see Refs. [zEUl]. □ 



By definition, bound entangled states require some entanglement for its creation, however, afterwards 
the entanglement cannot be distilled back. Naturally, such states are interesting, and the above Theorem 
implies that the search for PPT entangled states is of great interest in studying bound entanglement. The 
first example of a PPT ent angle d state was given in Ref. [Hl j] and in the meantime many other examples 
have been found 6^, 6(| 74 , 128| . 



In order to demonstrate how such states may be constructed, let us discuss a simple example from Ref. [65| . 
In a 3 x 3-system, one can consider the vectors 

|^ ) = i=|0)(|0)-|l)), |^) = _L(|0)-|1))|2), hM = ^=|2)(|l>-|2», 

IV>3> - - |2»|0>, |^4) = \(\0) + |1) + |2))(|0) + |1) + |2)). (23) 

These five product vectors form a so-called unextendible product basis: they are all pairwise orthogonal and 
there is not another product vector orthogonal to all of them, as can be seen by direct inspection. Therefore, 
the mixed state 

eBE = s (i-X>>M) (24) 

i=0 

has no product vector in its range, and must be entangled due to the range criterion. Since £>be = £>be ^ ^ s 
also PPT, and thus bound entangled. This state has also the minim al po ssible rank, since the minimal rank 
of a PPT entangled state in a x d^-system is max{d J 4, ds} + 1 129]. 
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Beside the mere existence, there are several results and questions, which make the phenomenon of bound 
entanglement interesting: 

(i) At first sight, entangled states that cannot be distilled seem to be not very useful, since they cannot 
be used for certain tasks like quantum key distribution via the Ekert protocol Q . It has been proved, 
however, that some bound entangled states can indeed be used for secure quantum key distribution 
130l 131 1, showing that entanglement distillation and secure key distillation are not equivalent. 



(ii) Bound entangled states require entanglement for their creation, but the entanglement can then not be 
distilled again. Similarly, one may ask, whether there are classical probability distributions between 
two or several parties t hat require secret communication for their creation, but they do not allow 
genera ting a secret key [l 32t | . It has been shown for the tripartite case, that such bound information 
exists |133| . Interestingly, the proof is based on a bound entangled tripartite quantum state, from 
which the classical probability distribution is derived. 

(iii) It has already been conjectured by Peres in 1999, that bo und entangled states admit a local hidden 



variable model and thus do not violate any Bell inequality 134j . We will discuss this open question in 
more detail in Section 5.3. 

(iv) Since PPT entangled states are bound entangled, the questi on arises, whether NPT states are always 
distillable. It has been conjectured that this is not the case 11351 1 13611 . E xamples of states have been 
found, which c anno t be distilled, if only few copies are available [135l Il36j ] . However, despite of s evera l 
partial results 137 1 the question remains open. For 2 x ^-systems, any NPT state is distillable [l35| . 

(v) The given example of a bound entangled state may seem artificial. An interesting question is whether 
bound entangled states occur under natural conditions, e.g., as t hermal states in spin models. For 
multipartite bound entanglement, this is indeed the case 138l 139j | . 

To conclude, all these results indicate that bound entangled states are an interesting object of study. Their 
construction, characterization and detection is a challenging and important task in entanglement theory. 



2.5. Entanglement witnesses 

Let us now come to a completely different type of separability criteria. The criteria from above have 
all something in common: at first sight they all assume that the density matrix is already known. They 
all require applying certain operations to a density matrix, to decide whether the state is entangled or 
notPl. There is, however, a necessary and sufficient entanglement c riterion in terms of directly measurable 

observables. These are the so called entanglement witnesses m m urn Em, winch we introduce now. 



Definition 8. An observable W is called an entanglement witness (or witness for short), if 
Tr(Wg s ) > for all separable g s , 

Tr(Wg e ) < for at least one entangled g e (25) 

holds. Thus, if one measures Tr(Wg) < one knows for sure that the state g is entangled. We call a state 
with Tr(Wg) < to be detected by W. 

The fact that entanglement witnesses are directly measurable quantities makes them a very useful tool for 
the analysis of entanglement in experiment. As we will see in the further course of this review, entanglement 
witnesses are one of the main methods to detect entanglement experimentally. 

For the further understanding, it is crucial to note that entanglement witnesses have a clear geometrical 
meaning. The expectation value of an observable depends linearly on the state. Thus, the set of states where 
Tr(Wg) = holds is a hyperplane in the set of all states, cutting this set into two parts. In the part with 
Tr{Wg) > lies the set of all separable states, the other part [with Tr{Wg) < 0] is the set of states detected 
by W. This scheme is shown in Figure 2. From this geometrical interpretation it follows that all entangled 
states can be detected by witnesses: 



9 Note that in principle one can decide whether a state violates the PPT criterion or not without knowing the state completely. 
For schemes to do that see Section 7.4 
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Fig. 2. Schematic picture of the set of all states and the set of separable states as nested convex sets and two witnesses, WW 
and W' 2 ). The red lines represent the hyperplanes, where Tr(Wg) = 0. The witness W' 1 ' is finer than the witness w( 2 ) [see 
also Eq. (32)] . 

Theorem 9 (Completeness of witnesses). For each entangled state g e there exists an entanglement 
witness detecting it. 

Proof. This theorem was proved in Ref. [Hoj ] . The idea for the proof relies on the fact that the set of 
separable states is convex and closed. Thus, for a point outside this set, there is a hyperplane separating 
this point from the convex set. For finite-dimensional spaces with a scalar product this is easy to see, but 
the statement holds als o for infinite-dimensional Banach spaces, where the proof requires a corollary of the 
Hahn-Banach-thcorem 

Although this theorem ensures that any entangled state can in principle be detected with an entanglement 
witness, the task remains to construct witnesses. This is not an easy problem, since solving this problem 
would also solve the separability problem. Indeed, a large part of this review is concerned with the proper 
construction and evaluation of witnesses for the multipartite case. 



2.5.1. Construction of witnesses 

Let us discuss in three simple examples, how witnesses for bipartite entanglement can be constructed. 
Typically, if a state violates some criterion for separability, an entanglement witness detecting the state can 
be written down. 

As the first example, let us take a state g e which is NPT. Then there is a negative eigenvalue A_ < of 
Qe A and a corresponding eigenvector \rj). Now 

W=\v)(v\ Ta (26) 
is a witn ess detecting g e : Due to the fact that for general matrices X, Y the relation Tr(XY TA ) = Tr(X TA Y) 



holdf|_^J, we have Tr(Wg e ) = Tr{\rf)(r]\ TA g e ) = Tr{\rf)(r]\g1 A ) = A_ < and Tr(Wg s ) = Tr(\n)(n\g^ A ) > 
for all separable g s , since they are PPT. 

More generally, if g e is entangled and detected by some positive map A, then 1(g) A(g e ) has again a negative 
eigenvalue A_ < for the eigenvector \r/). Then, W = I <g> A + (\rj)(rj\) is a witness detecting g e . Here, A + 
denotes the adjoint map, i.e., the map fulfilling Tr[A+(X)Y] = Tr[XA(Y)] for all A, Y. 

As a second example, let us consider the case when g e violates the CCNR criterion. Then, by defini tion, 
we have in the Schmidt decomposition in Eq. (10) J^k ^ k > 1- A witness is now given by Ref. \M,\iM 



W = l-J2 G k®Gk, (27) 
fe 

where the G h are the observables from the Schmidt decomposition. Clearly, for g e we have Tr(Wg e ) = 
1 — J2k Afe < 0. A general separable state g s — J^ki f-^iG^ ® Gf can be expressed in the basis of G^ B . 



1 This is obvious for the case X = A ® B and Y = C 55 D, for the general case it follows from linearity. 
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Then, Tr(Wg s ) = 1 — J^k Mkfe — 0, where the last estimate comes from the CCN R criterion and the fact, 
that a trace of a matrix is always smaller than the sum of its singular values |l43| . The witness in Eq. (27) 
can also be derived from the realignment map 



As a third example, witnesses can also be constructed from the consideration that states close to an 
entangled state must also be entangled. Therefore, one may try for a given entangled pure state \ip) to write 
down a projector-based witness like 

W = al - \t/>){ip\. (28) 

This witness can be interpreted as follows: The quantity Tr(g\tp)(if;\) — (ip\g\i/)) is the fidelity of the state 
in the mixed state g, and if this fidelity exceeds a critical value a, then the expectation value of the 
witness is negative and the state g must be entangled. 

The problem remains to calculate the smallest possible value for a, in order to ensure that W is still 
positive on all separable states. This can be done by taking 

a = . max Tr(g\^){^\) = max |(^|0)| 2 . (29) 

q is separable 10) — |a;Q9|o) 

Here, we used that a linear function takes its maximum on a convex set in one of the extremal points, and for 
the con vex s et of the separable states these extremal points are just the pure product states. Then, as shown 
in Ref. [l45j the maximum can be directly computed and is given by the square of the maximal Schmidt 
coefficient of Finally it is worth mentioning that the witness W in Eq. (28) does only detect states 
which are NPT therefore the construction is less general than the one in Eq. (26) [HJ It has the advantage, 
however, that it can be generalized to the multipartite case, as we will see later. 

Finally, let us mention other ways of constructing entanglement witnesses. Most of them are connected 
to special separability criteria. 

(i) As already mentioned, any positive, but not completely positive map gives rise to a construction of 
entanglement witnesses. As we will see soon, also a connection in the other direction holds. 

(ii) Many bound entangled states violate the range criterion in a way that for any product vector \a, b) G 
R(g) in the range of g we have already \a, b*) ^ R(g TB ). States with such an extremal violation are 
called edge states, and if an edge state g is PPT, then it lies on the boundary between the PPT and 
the NPT states. For such a state, one can write down a witness of the type 

W = oil - P - Q Tb (30) 

where P (Q) deno tes the projector onto the range of g (g TB ) and a is defined as a = max|0) = | aj 6 ) (<j>\ (P— 
Q Tb )\4>) < 2 [l4l| . This construction guarantees that W is a valid witness, and g is clea rly detect ed. 
The parameter a can either be estimated analytically , or calculated numerically [H, Il48l . Il49j . 

(iii) The semidefinite program in the criterion of the symmetric extensions gives as a byproduct an 
entanglement witness, if no symmetric extension was found. 

(iv) Entanglement witnesses can also be constructed from geometrical considerations: If g e is an entangled 
state, and a is the closes separable state to g e [in Hilbcrt- Schmidt norm, \\X\\ :— yjTr{X^ X)], then 

W = ^(a-g e + Tr[a( Qe -a)}-l) (31) 

with JV = \\g e — °1 is an entanglement witne ss, since the set of states with Tr(Wg) = is the 
hyperplane orthogonal to the line from g e to a [l50L Il5lj . The normalization is chosen in s uch a wa 



that its expectation value for g e equals (up to the sign) the distance to the separ able states I151U15 
Such geometric entanglement witnesses have been applied to different situations 153L 154 , 155l . 156| . 



(v) The criteria based on variances and covariance matrices also result in certain types of entanglement 
witnesses. These witnesses, however, are nonlinear, we will discuss them in Section 7.3. Furth er, special 
constructions of witnesses have derived from estimations of entanglement measures [l57j (see also 
Sect ion 7.2). Also, from a given witness detecting some state further witnesses can be constructed 



11 Th is follows from the fact that the witness W = 1/d — \ip) (ij}\, where is maximally entangled, does not detect NPT states 
Il46ll , combined with an argument as in Theorem 3 of Ref. |l47ll . 
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2.5.2. Optimization of witnesses 

Every entanglement witness can, by definition, detect some entangled states. However, as can directly be 
seen in Fi g. 2, some witnesses are better in this task than others. In this sense, one can optimize entanglement 
witnesses [141| . First, one calls a witness finer than another witness W^ 2 ** if detects all the states 
detected by and also some states in addition. This is the case, if 

yy( 2 ) = + P, (32) 

where P is a positive operator. Consequently, for any state Tr(gW^) < Tr(gW^) holds. 

Furthermore, a witness W is called optimal, if there is no other witness, which is finer than W. This 
implies that for any positive operator P the observable X = W — P is not a witness anymore. From this it 
can be easily seen that W is optimal, if and only if the product vectors \(f>i) = \aibi) with ((j>i\W\4>i) = span 
the whole space. A necessary condition for a witness to be optimal is that it "touches" the set of separable 
states (see Fig. 2), i.e., there must be a separable g with Tr(gW) = 0. This is, however, not a sufficient 
condition, witnesses which fulfill this condition are sometimes als o cal led weakly optimal. 

The optimization of a given witness can be done as follows [l4l| : If one finds a positive operator P, 
which fulfills (a,ibi\P\ai, bi) = for all product vectors with (a,ibi\W\ai, bi) — 0, then W = W — \P is an 
entanglement witness finer than W if 

A<A :=mf min [_Lw a -i=] (33) 

a eigenvalues \/ ± J a -v/ } J a 

with the notation X a = (a\X\a) G B(Hb) for \a) £ Ha and X £ {W,P}. This can be iterated, until the 
product vectors with (aibi\W\aibi) = span the complete space. However, the optimization in Eq. (33) is 
not straightforward, as it is equivalent to checking whether inf| a )| b ) (a, b\[W — XP] \a, b) > 0. In fact, such an 
optimization for general observables is equivalent to the separability problem [1591 ] , and it can be shown that 
the question whether a given state is separable or no t is e quivalent to the question whether some witness in 
a higher dimensional space is weakly optimal or not [l6C| . 



2.5.3. The Choi-Jamiolkowski isomorphism 

Now we explain a close connection between witnesses and positive maps. This connection is given by 
an isomorphism between quantum states on a space of matrices, and maps bet ween such spaces . In the 
literature, this isomorphism is known as the Choi-Jamiolkowski isomor phism [H EM EH E63J OH 

According to this, an operator E in B{TLb)®B(TLc) corresponds to a map e that maps states from B(Hb) 
to an element of B(Tlc)- This map is given by 

e{g)=Tr B {Eg T < 



lc). 

The inverse relation reads 

E= (I B ^ £ )(|0 + )(0+|), 

where Hb' — "H B and \4> + ) — 
the fact that \<j) + )((j) 



(34) 



(35) 



= £ fc (Gf ) T 



be expressed as 



is a non-normalized maximally entangled state on TLb' ® Ti.B- Using 
g> Gf (see, e.g., Ref. [TT^I) for any orthogonal basis G*f ^ B , this may 



(36) 



from which it is also easy to see that Eq. (34) is the inverse of Eq. (35). Physically, these relations can be 
interpreted as follows. If the map e can be realized, then the state E can be prepared by applying the map on 
one part ofamaximally entangled state. On the other hand (due to the relation Ttb'b[Xb' <8>Yb|</> + )(<^ + |] = 



12 To our knowledge, an isomorphism of this type was first investigated by J. de Pillis [16ll |. who used a map as in Eq. (36) 
witho ut t he transposition on the first subsystem and proved that for this map e is positive if E is positive semidefinite. In 
Ref. [l62l |. it was proved for the same map that e is positive, iff E is positive on all product states, i.e., it is a n ent anglement 
witn ess, and the explained version and the proof of main properties in Theorem 10 are given in Refs. |69l Il63| . See also 
Ref. H6l. 
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Trs' b[{Xb') t <£> Yb\4> + ) ((/> + \ Tb '] = TrB[X T Y]) one can view the map e as generated via teleportation with 
the state E. 

The Choi-Jamiolkowski isomorphism is useful for several purposes. For instance, it can be used to define a 
distance measure for maps as quantum oper ation s by looking at the distance of the corresponding states. This 
distance has then several natural properties [l65l | . It can also serve to investigate the entanglement generating 
capacity of positive maps on a bipartite system [l66j . For the theory of witnesses, the followi ng theore m 
concerning the relation between the map e and the state E is relevant (for a proof see Refs. [69l. [l63l[l67t ): 



Theorem 10 (Properties of the Choi-Jamiolkowski isomorphism). The Choi-Jamiolkowski iso- 
morphism has the following properties: 

(i) The map e is completely positive, iff E is a positive semidcfinitc operator. 

(ii) The map e is positive but not completely positive, iff E is an entanglement witness 

(iii) The map e is a decomposable map, iff E is a decomposable entanglement witness. 
Concerning the last point, note that an entanglemen t wi tness is called decomposable, if it can be written 



as W = Pi + Pf B , where Pi, Pi are positive operators 141 1. Clearly, such a witness cannot detect any PPT 
entangled state. 

Due to this isomorphism, any witness, which detects some PPT entangled state gives rise to a positive, 
but not completely positive map, which is non decomposable. In this way, entanglement witnesses can be 
used to construct new and interesting positive maps. 

The question arises, whether the map e(W) detects more states than the witness W. Indeed, one can 
directly see that if A : B(Hb) — > B(Ha') is a map detecting the sta te q ab, then it detects all the states 
(Fa <8> ~§-b)qab(F a ® Is), where Fa is an arbitrary invertible matrix [l67j |. On the other hand, if Wab is a 
witness with Tr (Wab Qab) < 0, then the adjoint map [e(>V)] + : B(Hb) —* B(Ha') detects the state, from 
which it follows that the map e + detects all states which are detected by the witness 

W = [(F a )- 1 ®1b]W[(f\)- 1 ®1 b ]. (37) 



This is a large r set of states, for examples see Refs. 168l 169| . The operator Fa corresponds to local filters 



170l . 171 . 172]. Such transformations do not change the fact whether a state is entangled or not, hence they 



may be used to simplify entanglement detection. For instance, it has been shown that any state g of full 
rank can be simplified via filters to (assuming cIa < ds) 

q ~ q = (F A <8> F b )q(F a ® F B y = — — + V ® Gf , (38) 

d A d B ^ 

where the G k are traceless orthogonal observables fl7ll[l7i Such normal forms under filtering operations 
have been used to improve given separability criteria [109l . lll4l ] and we will see in Section 6.8 how filters can 
be used to improve witnesses systematically. 

Finally, another interesting consequence from the Choi-Jamiolkowski isomorphism is that any witness can 
be written in the form Eq. (35). This form is similar to Eq. (26), but with a different positive map. This 
idea will be later used for the construction of nonlinear witnesses. 



3. Multipartite entanglement 

In this Section we discuss the structure of entanglement when more than two parties are involved. It 
will turn out that this structure is much richer than the structure of entanglement in the bipartite case. 
Especially, it will turn out that several inequivalent classes of entanglement exist. We first discuss different 
notions of entanglement and separability for the case of three qubits. Then, we discuss the case of general 
multipartite systems. After that, we introduce the main families of multipartite entangled states, such as 
GHZ states, Dicke states and graph states, which are interesting because of their entanglement properties 
or their usefulness for applications. Finally, we discuss the construction of entanglement witnesses for the 
multipartite case. 
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3.1. Entanglement of three qubits 



3.1.1. Pure states 

Let us first consider pure three-qubit states. There are two different types of separability: the fully separable 
states that can be written as 

\^ s )a\b\c = \a) A ®\(3) B ® \l)c, (39) 

and the biseparable states that can be written as a product state in the bipartite system. A biseparable state 
can be created, if two of the three qubits are grouped together to one party. There are three possibilities of 
grouping two qubits together, hence there are three classes of biseparable states. One example is 

\<p hS )A\BC = W) A ® \S)bC (40) 

The other possibilities read \4> hs )B\AC = \P)b <8> \S)ac and \4> hs )c\AB = \l)c <£> \5)ab- Here, \5) denotes a 
two-party state that might be entangled. Finally, a pure state is called genuine tr ipartite en tangled if it is 
neither fully separable nor biseparable. Examples of such states are the GHZ state jlzl ^7^ 

|GtfZ 3 ) = -^(|000> + |m», (41) 

and the so-called W state [Izl Oil] 

\W 3 ) = ^(|100) + |010) + [001). (42) 

From a physical point of view, the generation of fully separable or biseparable states does not require 
interact ion of all parties. Only for the creation of genuine tripartite entangled state all three parties have to 
interact 
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The genuine entangled three-qubit states can, however, be further divided into two inequivalent classes in 
the following way: Given two three-qubit states, \4>) and one can ask whether it is possible to transform a 
single copy of |</>) into with local operations and classical communication, without requiring that this can 
be done with certainty. The se operations are called stochastic local operations and classical communication 
(SLOCC). It turns out \l7^ that \4>) can be transformed into iff there exist invertible operators A, B,C, 
acting on the space of one qubit with 

\ip)=A®B®C\<t>). (43) 

Since the operators A, B, C are inver tible , this defines an equivalence relation with a clear physical meaning. 

Surprisingly, it was proved in Ref. |l76| that there are two different equivalence classes of genuine tripartite 
entangled states, which cannot be transformed into another by SLOCC. One class, the class of GHZ states 
is represented by the GHZ state \GHZ 3 ) defined in Eq. (41). The other class, the class of W states can be 
transformed via SLOCC into \W 3 ) given in Eq. (42). In this sense there are two different classes of tripartite 
entanglement. There are much more pure GHZ cl ass states than W class states: By local unitary operations 
one can transform any pure three-qubit state int d 14 1 

|V> = Ao|000) + Aie^lOO) + A 2 |101) + A 3 |110) + A 4 |lll), (44) 

where A^ > 0, Af — 1 an d G [0; tt], see Ref. [l82| . Thus, six real parameters are necessary to characterize 
the nonlocal properties of a pure state. For the W class states, however, 9 — A4 = holds, which shows that 
they are a set of measure zero in the set of all pure states. 

Physically, there are also differences between the two classes: On the one hand, the GHZ state is maximally 
entangled and a generalization of the Bell states of two qubits. For instance, for the most known Bell inequali- 
ties the violation is maximal for GHZ states [l83| . On the other hand, the entanglement of the W state is more 



13 The notion of genuine multipartite entanglement was (in the context of Bell inequalities) already anticipated in Ref. [17: 
see also Section 5.2.3. 

14 This is a generalization of the Schmidt decomposition to three qubits, see Refs. [i7ll . [l7^ . [l7i .[i79l . CT . [T8l for other results 
on this problem. 
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Fig. 3. Schematic picture of the structure of mixed states for three qubits. The convex set of all fully separable states (fs) is 
a subset of the set of all biseparable states (bs). The biseparable states are the convex combinations of the biseparable states 
with respect to fixed partitions, sketched by the three different leafs. Outside are the genuine tripartite entangled states, the 
W class and the GHZ class. 



robust against particle losses: If one particle is lost in the GHZ state, the state qab = Trc{\GH Z3) (GH Zs\) 
is separable, for the W state the resulting reduced density matrix qab = Trc{\W)(W\) is entangled. In- 
deed, it can be shown tha t the W s tate is the state with the maximal possible bipartite entanglement in the 
reduced two-qubit states [l76l . fl8^ . 



3.1.2. Mixed states 

The classification of mixed states according to Refs. [l85l Il86l | is similar to the definition of bipartite 
entanglement via convex combinations. We define a mixed state g is as fully separable if g can be written 
as a convex combination of fully separable pure states, i.e., if there are convex weights pi and fully separable 
states \<f)f) such that we can write 



Pi\<t>t)(<f>t\- (45) 



A state g hs that is not fully separable is called biseparable if it can be written as a convex combination of 
biseparable pure states: 



Q hs = J2Pi\$ S )($ S \- ( 46 ) 



Note that the biseparable states \$ s ) might be biseparable with respect to different partitions. We will see in 
Section 3.2.2 why it is reasonable to exclude mixtures of biseparable states with respect to different partitions 
from the genuine multipartite entangled states. Of course, one may define three classes of biseparable mixed 
states that are biseparable with respect to a fixed partition. 

Finally, g is fully entangled if it is neither biseparable nor fully separable. There are again two classes of 
fully entangled mixed states: A fully entangled state belongs to the W class if it can be written as a convex 
combination of W-type pure states 

e w = $>ICXCI, (47) 

i 

otherwise it belongs to the GHZ class. One can show that the W class forms a convex set inside the GHZ 
class. Also, in contrast to t he ca se of pure states, the set of mixed W class states is not of measure zero 



compared to the GHZ class |185j . 



Again, this classification can be represented in a schematic picture of nested convex sets (see Fig. 3). 
Here, it is important to note that this picture is, of course, only a schematic picture, which does not take 
all properties into account. For instance, it has been shown that there exist states that are biseparable with 
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resp e ct t o each fixed partition, however, they are not fully separable (for some examples see Refs. [65, 138 , 

The question to which class a given mixed state belongs, is, as the separability problem for the bipartite 
case, not easy to decide. Some conditions will be discussed in Section 3.5. Methods to distinguish between 
the mixed W class and GHZ class are, however, very rare. As we will see later, witnesses give the possibility 
to prove that a state belongs to the GHZ class. However, it is not clear how one can show that a state is 
tripartite entangled and belongs to the W class. This can not be done with witnesses, since they are designed 
to show that a state lies outside a convex set, they fail to prove that a state is inside a convex set. 



3.2. Multipartite entanglement distillation 



3.2.1. Multipartite distillability 

It is also possible to define distillability for three and more qubits. For instance, one may aim at the 
distillation (also called purification) of the GHZ state \GHZs), that is, one asks whether the transformation 

q®q®...®q^ ^? \GHZ 3 ) (48) 

k copies 

is possible. If this is possible, the state g is called multipartite distillable. The fact that the target state is 
a GHZ state, is not important for the question, whether or not g is distillable, since it is always possible to 
transform a GHZ state in a W state and vice versa, if many copies of the state are available. 

M any p rotocols for the distillation of different multipartite states have been proposed (for a review see 
Ref. t 188]). Also, multipartite bound entangled states exist. For instance, the above mentioned states that 
are separable with respect to each bipartition, but not fully separable, are multipartite bound entangled. 



3.2.2. Multipartite distillability vs. multipartite entanglement 

Let us compare the notion of multipartite distillability with the notion of genuine multipartite entan- 
glement and their importance for the analysis of experiments. First, one has to stress that both notions 
are inequivalent. There a re ge nuine multipartite entangled states that are not distillable. Examples of such 
states were given in Ref. [128J. 

To investigate the relation better, consider the mixed three-qubit state 

Q = ^(I^ + X^ + Ub ® |0)(0| c + \cj> + ){4> + \AC ® |0)(0| B + |0 + )(0 + |sc ® |0)(0U), (49) 

where = (|00) + \ll))/\/2 is a two-qubit Bell state, g is a mixture of biseparable states with respect 
to different partitions, and therefore biseparable. However, as can be easily checked, g is entangled with 
respect to each bipartition, and also any reduced two-qubit state is NPT and hence entangled and two-qubit 
distillable. So, if many copies of g are available, the three parties can first distill singlets between A — B 
and A — C. Then, Alice may prepare a GHZ state locally in her lab, and teleport two of its parties via 
the existing singlets to Bob and Charly. By this simple procedure one sees that g is multipartite distillable, 
showing that multipartite distillability and genuine multipartite entanglement are really inequivalent. 

The important point for our discussion in the above example is that the state g is extremely easy to 
prepare experimentally. First, one prepares a product state 1 000) and then one entangles randomly two of 
the three qubits. Equivalently, one can also prepare the state |</> + )(</> + |ab ® 1 0) (0 1 , and then one forgets 
which pair of qubits was the pair A, B. 

From this example, we can conclude three things. First, multipartite distillability is not a good criterion 
for claiming success of an experiment, since some multipartite distillable states are very easy to prepare. 
Second, since g is also NPT with respect to each bipartition, it is not sufficient to investigate all bipartitions, 
and the PPT criterion is a very limited tool to investigate multipartite entanglement. Finally, this example 



a This counteri ntuiti ve behavior can lead to surprising applications, like the possibility of distributing entanglement using 
separable states Il87ll . 
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justifies that we did not define biseparability with respect to a fixed partition, and considered also mixtures 
of biseparable states that are biseparable with respect to different partitions as not highly entangled. Due to 
these facts, we will mainly investigate how genuine multipartite entanglement can be verified in experiments. 



3.3. Entanglement classes for the general case 



Here we shortly discuss the classification of entanglement for general multipartite systems. This is a simple 
generalization of the three-qubit case: First, we distinguish different types of entanglement for pure states. 
Then we extend this definition to mixed states by considering convex combinations of pure states. 

Let us assume that a pure JV-partite state \ip) is given. We call this state fully separable if it is a product 
state of all parties, that is if 

N 

M = ®l&) (50) 

i=l 

holds. A mixed state is called fully separable if it can be written as a convex combination of pure fully 
separable states, that is, it can be written as 

Q = Y.p^k ) ®i ) ®-®i f) - ( 51 ) 

fe 

If a pure state is not fully separable, it contains some entanglement. Again, as in the three-qubit case, this 
does not have to be true TV-partite entanglement. Thus, we call a pure state m-separable, with 1 < m < N, 
if there exist a splitting of the N parties into m parts Pi, P m such that 

m 

= <S> I'M Pi ( 52 ) 

i=l 

holds. There are (m N )/m\ possible partitions of the N parties into m parts. Note that an m-separable state 
still may contain some entanglement. Again, we call mixed states m-separable, if they can be written as 
convex combinations of pure m-separable states, which might belong to differen t partitions. O ne can refine 
these definitions, by considering separability with respect to fixed bipartitions 186l 189l 190l |. Finally, we 



call a state truly TV-partite entangled when it is neither fully separable, nor m-separable, for any m > 1. 

The above classification distinguished different classes, according to the question: "How many partitions 
are separable?" . Alternatively, one m ay ask the ques tion "How many particles are entangled?" A classification 



in this sense was proposed in Refs. [19ll . Il92l . Il93j ]. According to that, a pure state \ip) is said to contain 



only m-party entanglement, if can be written as 



A" 



IV>)=(g)l&) (53) 

i=l 

where the \<f>i) are states of maximally m qubits. Consequently, K > N/m has to hold. If is not of this 
form, it contains (at least) m + 1-party entanglement. Again, this definition can be extended to mixed states 
via convex combinations. A mixed state that is fc-producible requires only the generation of fc-party pure 
entangled states and mixing for its production. Consequently, a mixed state g contains fc-party entanglement, 
iff the density matrix cannot be obtained by mixing pure states that are (k— 1) -producible. This classification 
of separability classes does not play a significant role in the analysis of experiments, since there typically 
only few particles are considered, and the aim is to prove genuine multipartite entanglement, which coincides 
with the statement that g contains TV-partite entanglement. H owever, fo r spin models with a large number 
of particles, this notion can be a reasonable property to study [HI [nil (see also Section 8.3). 

Another characterization is possible based on the number of uncntangled particles. An TV-qubit pure 
product state contains n unentangled particles if it can be written as = l \4>k) <8> l^n+i.-jv) where 
\4>k) are single-qubit states, while |^„ + i.jv) is a (TV — rt)-qubit state vector. Consequently, a mixed state 
contains at most n unentangled spins, if it can be constructed mixing pure states with n or less unentangled 
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spins. There are entanglement conditions that can provide upper bounds on the number of unentangled 
spins (see Sec. 8.1.2). If there are no unentangled spins, the state is not necessarily genuine multi-partite 
entangled: It can be a tensor product of two-qubit entangled states. However, creating very many (e.g., 10 5 ) 
two-qubit singlets can be a goal of quantum control in itself. 

Finally, it is worth mentioning that a classification of the pure truly TV-partite entangled states like in 
the three-qubit case is not straightforward. At least a classification via equivalence classes under SLOCC is 
not simple, sinc e it has b een shown that already for four qubits there are infinitely many equivalence classes 
under SLOCC 194 , 1 1 951 ] . This continuous set can then be separated into nine or eight di ffere nt classes. 
Other app roaches try to identify different classes via a generalization of the Schmidt rank [l9tj ] or tensor 



ranks [1971 ] . but there it is difficult to determine in which class a pure state lies. 



3.4. Families of multipartite entangled states 



In order to further investigate multipartite entangled states, it is useful to identify interesting families of 
multipartite pure states. Here, one might be interested in maximally entangled states, or states that are 
useful for certain applications. For the three-qubit case, we already know the GHZ state and the W state 
as such examples. There are, however, several families of interesting multi-qubit states with an arbitrary 
number of qubits: 



3.4.1. GHZ states 

The GHZ state for N qubits is defined as 

1 



\GHZ_ 



V2 



(10) 



1 



(54) 



and is a simple generalization of the three-qubit GHZ state. It is a superposition of two maximally dis- 
tinct states, and therefore also sometimes called a "Schrodinger cat state". GHZ states have been studied 
intensively for a long time and it turned out that they are useful fo r many applications. These include 
entanglement enhanced s pectr oscopy and quantum me trolo gy 198 . 199| . quantum secret sharing [ 200l| . open 
destination teleportation 201 1, quantum computation [2021 ] and cryptographic protocols 203. 204lirFurt her- 
more, it has been shown that GHZ states lead to an extremal violation of local realism 174, Hlioi. In 



fact, it can be s hown that for many Bell inequalities, the unique state leading to a maximal violation is the 
GHZ state (l8^ . 

These properties make the GHZ states interesting from many perspectives. In some sense they may be 
considered as the maximally entangled multi-qubit states. Therefore, many experiments aim at th e generat ion 
of GHZ states, and GHZ states wi th up to ten qubits have already been prepared using ions @, 12071 1208| or 

photons d Ea uni, m, nm . 



3.4.2. Dicke states and W states 

These states have been investigated already in 1954 by R.H. Dicke, while studying light emission from a 
cloud of atoms [212l |. Dicke found that, if the atoms are in certain highly entangled states, then the intensity 
of radiation is much larger than if the atoms were emitting light independently, starting all from an excited 
state. In particular, the intensity scales quadratically with the number of atoms, while for independent atoms 
it scales linearly. 

Dicke states, in general, are simultaneous eigenstates of the collective angular momentum operators J z — 
i/ 2 Efc^ fe) and J 2 Of] Typically, for a certain eigenvalue for J z and J 2 there are several eigenstates and 
additional parameters are needed to distinguish them from each other. Since very often many-particle systems 
start from a symmetric initial state and they evolve under symmetric dynamics, perhaps the most important 
examples for quantum information arc the symmetric Dicke states. Such states are all eigenstates of J 2 with 
the maximal eigenvalue y(t an< ^ they can ^ e characterized uniquely by their eigenvalues j z for J z . 



We will discuss these operators and their properties later in Section 8.1 in more detail. 
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Beside j 2 , it is also common to use for their characterization k = N/2 — j z . With this, an (N, k) symmetric 
Dicke state on N qubits is defined as an equal superposition of k excitations, i.e., 



\D 



k.N) 



^Pj\\l) 9k ®W> N - h \ (55) 



where Y^j denotes the sum over all possible permutations of the qubits. 

The W state corresponds to the Dicke state l-Di.,3) and is therefore an example of a Dicke state. Another 
important Dicke state is the four-qubit state \D 2A )'= (|0011) + |0101) + |0110) + |1001) + |1010) + |1100))/\/6. 
In general, the state \D n ) is the state for which the superradiance is the strongest. Also, since its overlap 



2 ,N> 



with biseparable states is the minimum \ for large N |213| . for an experimental investigation of genuine 
multipartite entanglement only a fidelity of h is required. Finally , it s hould be noted that the entanglement 



in the Dicke state \D%,a) is relativ ely robust aga i nst d ecoherence [214j . Dicke states and W states have been 
prepared in many experiments [ci l208l I215I 216l 217 1. We will discuss two of them in Section 6.8 and 8.2. 



3.4.3. Graph states and cluster states 

Let us come now to the family of graph states 218 . 219l . 220l |. These states are in a mathematical way 



defined as follows: First one considers a graph, i.e., a set of N vertices and some edges connecting them. 
Examples of graphs are shown in Fig. 4. For each vertex i the neighborhood Af(i) denotes the vertices that 
are connected with i. Then one can associate to each vertex i a correlation operator or stabilizing operator 
9i by 

9i = Xi (g) Zj. (56) 

jeA/"(») 

Here Xi,Yi,Zi denote the Pauli matrices <7 x ,a y ,a Zl acting on the i-th qubit; we will use this notation often 
when working in the graph state formalism. For instance, for the three vertex graph No. 2 in Fig. 4(a), the 
stabilizing operators are 171 = Z1X2I3, <?2 = X\Z 2 Z^ and 53 = Z1I2X3. The graph state \G) associated with 
the graph G is the then the iV-qubit state fulfilling 

9i \G) = \G), iori=l,...,N. (57) 

Clearly, \G) is then also an eigenstate of all possible products of the gk, the commutative group of this 
operators is called the stabilizer |22l| . Since the gi form a maxi mal set of commuting observables, the graph 
state \G) is uniquely determined by these eigenvalue equations P^l 

The physical meaning of the graph is that it describes the perfect correlations in the state \G), since 
(gi) = (Xj (^jgjvm Zj) = 1- At the same time, the graph denotes a possible interaction history leading to 
\G) : the state \G) can be produced from the product state [(|0) + | l))/v / 2] <8>JV by an Ising t ype interaction 
Uki = exp{— it(tk — Zi)(lk — Zi)} acting for the time t — it /A between the connected qubitsl 18 1 Finally, the 
graph state can also directly be written in the stabilizer elements as 

N 

\G){G\ = H±Y*. (58) 
fe=i 

There are several interesting properties of graph states. First, it should be noted that the GHZ states 
belong to this class. If we consider the state in Eq. (54), it can be easily checked that this is an eigenstate of 

gx = X 1 X 2 ...X N ; g k = Z x l 2 ...Z k ...l N ; for k = 2, N. (59) 

Up to some local rotations, these are nothing but the stabilizing operators of the star graphs [Nos. 2, 3, and 
5 in Fig. 4(a)]. Therefore, graph states can be seen as an extension of the GHZ states. 



17 For a given graph, one can consider a whole family of 2 N states that are eigenstates of the gi with eigenvalues ±1. This 
orthogonal basis is the so-called graph state basis, we will use it in Section 6.6 in more detail. 

18 If the interaction time depends on k and I, the resulting state is a so-called weighted graph state, see Ref. \'22'S i for discussions 
and applications. 
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Fig. 4. (a) Several types of graphs for up to five qubits. The shown graphs represent all inequivalent graph states up to five 
qubits [220j | ■ (b) The graph for a five qubit linear cluster state, (c) The graph for a 2D cluster state (for the case of 15 qubits). 
This type of graph states is a universal resource for measurement based quantum computation. 



Second, it is important to note that different graphs may lead to graph states, which can be transformed 
into each other by local unitary transformations, implying that the entanglement properties of the states 
are th e same. I n general, it is not known yet, under which conditions two graphs yield equivalent graph 
states 2231 . |224| . However, it is known that some graph transformations (called local complementation) do 
not chang e the graph state, and for a small number of qubits (N < 8) a full classification has been achieved 



220LI22 



This leads directly to the question, what other graph states exists beside the GHZ state. For 2 qubits, 
there is only one graph [no. 1 in Fig. 4(a)], corresponding to a Bell state. For three qubits, there is again only 
one independent graph state, corresponding to the three-qubit GHZ state (No. 2). For four qubits, there are 
already two independent graph states: the four-qubit GHZ state (No. 3), and the four-qubit cluster state, 
given by 



\CU) = i(|0000) 



10011) + |1100) - |1111», 



(60) 



which corresponds (up to a local change of the basis) to the graph no. 4. For five qubits there are already 
four different graph states (nos. 5-8), and with increasing number of qubits, the number of inequivalent 
graph states increases rapidly. 

Third, it should be stressed that graph states are of eminent importance for applications and protocols 
of quantum information processing. First, a special graph state, the 2D cluster state [see Fig. 4(c)] can be 
used as a resource for measurement based qu antu m computation @, |226| . In addition, graph states occur a s 
codewords of quantum error correcting c odes 227 1 , violate certain Bell inequalities maximally 147l [228. 229] , 
and are al so ro bust against decoherence [2181 ] . Finally, graph states can be used for the simulation of anyonic 
statistics 230] and can be related to the computation of the partition function of classical spin models 231 1 . 

Due to this importance , it is not surprising that graph states have been prepared in many experiments 
(9I. I23I . I23I . I234L I235L l236| . We will discuss in Section 6.7 an experiment for the generation of the four-qubit 
cluster state. 



3.4.4. Multi-qubit singlet states 

Another interesting family of multi-qubit states are multi-qubit singlet states. These are pure states 
that are invariant under a simultaneous unitary rotation on all qubits, i.e., we have ((7)® Ar | , i/>) = e l ^\ip) (or, 
equivalently [IV'XV'li (U)® N ] = 0). Such states exist only for the case that N is even. For two qubits, the only 
state of this form is the two-qubit singlet state ■ For four qubits, these states form a two-dimensional 
subspace of the total Hilbert space, spanned by the vectors Y$>\) = <g> \tp~) and 
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|\£r 2 ) - -^[10011) + |1100) - |(|01) + HO)) ® (|01) + |10))] . (61) 



This is the so-called four-qubit si nglet state [237j . In general, there are d N = N\/[(N/2)\(N/2 + 1)1} linearly 
independent states in this family [2381 . |239| . 

The importance of these states comes from the fact that they are by definition robust against certain 
types of noise, where all qubits undergo the same (random or unknown) unitary transformation. They lie 
in a so-called decoherence free subspace 24Cll 24 lj . This can be used to encode quantum information, e.g., 



one qubit can be encoded in superpositions of the four-qubit states and | \&2) - 

Experimentally, it is remarkable that the state |\&2) can simply be prepared by a seco nd order process in 



parametric down-conversion [237j | (see also Section 6.5). This has been done in Ref. [242j |. also pr otocols lik e 
secret sharing or the encoding of quantum information as described above have been performed [243l . I244] . 

The expression singlet is also used in another context, denoting states with zero total angular momentum. 
That is, these are the states that are the eigenstates of the total angular momentum components J x , J y 
and J z with eigenvalue zero. Such states appear very often in condensed matter physics, for example, the 
ground state of the anti-ferromagnetic Heisenberg chain is such a multi-qubit singlet state. We will discuss 
entanglement detection in the vicinity of such states in Section 8.1. 



The connection between the two definitions can be seen as follows ! 19 1 One can rewrite a simultaneous 
multilateral unitary transformation as 

U® U <&...<& U = exp (ia ■ J), (62) 

where J = (J x , J y , J z ) and a is a three-element real vector. Hence, we can see that if a pure state is an 
eigenstate of expect./) for any a, then it is an eigenstate of Jk for k = x,y, z, and vice versa. Thus, the two 
definitions coincide for pure states. 

3.4.5. Werner states 

In the general case, iV-qudit mixed states invariant under U ® U <£> ... <£> U are called Werner states [Hj]. 
They do not necessarily have zero total angular momentum variance. Such multi-qudit states naturally arise 
as the result of the twirling operation 

Twirl[ e ] := / U® N g{U^)® N fi(dU), (63) 



where the integral is over the unitary group of d x d m atric es U(d), and \x is the Haar measure. Multi- 
qudit twirling can efficiently be realized experim ental l y | 245l| . The resulting Werner states are the linear 



combination of relatively few basis matrices [44], 12461 . |247| . For example, for N = 2 there are two such 
matrices for any dimension d, while for N = 3 there are 5 for d = 2 and 6 for any d > 2. Thus, such states 
can be characterized with few and simpler measurements than the measurements needed for full tomography. 
For example, for the two-qubit case the single free parameter of a two-qubit Werner state can be obtained 
by measuring (a^a^). Twirling has also been suggested for simplifying process tomography [2481 ]. 

3.4.6. Maximally entangled states for comb monotones 

As already mentioned, it is a natural question to consider the maximally entangled states for N qubits. 
This, however, depends on the entanglement measure chosen: As we will discuss later, there are several 
inequivalent measures of multipartite entanglement, and the maximally entangled states differ for different 
measures. 

In Ref. (Hi, l250j a class of entanglement measures has been introduced, which are based on anti-linear 



operators and combs. There are three different measures of this type for four qubits, resulting in three 
different maximally entangled states. The first one is the GHZ state, the second one the cluster state, and 
the third one is the state 

| X ) = -pfv^lllll) + |0001) + |0010) + |0100) + |1000)). (64) 
v6 



! We thank Sofyan Iblisdir for discussions about this point. 
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This state has also been mentioned in Ref. 25 1| , where it has been noticed that it is the symmetric four-qubit 
state that maximizes certain bipartite entanglement proper ties. Also, \x) is probably the most entangled 
symmetric state for the geo metri c measure of entanglement [252| (see also Section 4.2.6). 

As it was shown in Ref. |250l | one can generalize this state to five and six qubits, where it is also a 
maximally entangled state for some comb measure. It is not yet clear, which physical properties follow from 
this fact. Also, the generation of this state is still an experimental challenge. 

Finally, it should be added that there are further type s of states, which deserve attention. For example, 
there are states for optimal cloning via teleportation 253], or the family of multipartite states which can be 
maximally entangled with local auxiliary systems 254j |. 



3.5. Separability criteria for the multipartite case 

As for the bipartite case, one may ask how to decide whether a given mixed multipartite state is entangled 
or not. This problem is not so well understood as for the bipartite case, however, several criteria are known: 

(i) Permutation criteria. As discussed already for the bipartite case, one can view the PPT as well as 
the CCNR criterion as norm conditions on a density matrix where some indices are permuted [see 
Eq. (21)]. This approach can be generalized to more parties 0, [62|, [63| . One first writes the density 
matrix in some product basis, g = J2i 1J1: ...,i N ,j N fti,ji,-.,iw,iJvKi)(ii| ® ■■■ ® |*JV>(jJv|i and then for 
separable states we have 

\\Qw(ii,jl,...,iN,jl,)\\l < 1 (65) 

where 7r(...) is an arbitrary perm utat ion of the indices. This approach generalizes the CCNR and PPT 



criterion to the multipartite ease l 20 1 Clearly, not all permutations ir result in different criteria. As we 
have already seen, for two parties there are only two inequivalent permutations, resulting in the PPT 
and CCNR criterion. Further, it has been shown in Rcfs. m [63j that for the case of three parties 



there are six independent permutations, and for four parties there are 22 independent permutation 
criteria. All these criteria can only rule out full separability, 
(ii) Quadratic Bell-type inequalities. It was already demonstrated in Ref. [2551 ] that quadratic Bell-type 
inequalities can be a useful tool for the investigation of multipartite entanglement. Recently, a general 
method to derive quadratic i nequaliti es for the investigation of the different classes of multipartite 
entanglement has been found (l90l . |256| . This method can be sketched as follows: for two qubits and i = 
0,1 one considers the families of observables^^ = {Ai®A 2 ^Bi®B 2 ) /2, B { p = (Bi(g>A 2 ±Ai(g>B 2 )/2, 
C\ 2) = (Ci ® 1 2 ± li O C 2 )/2, and l\ 2) = (l x ® 1 2 ± C x ® C 2 )/2, where A k ,B k , C k denote orthogonal 

(2) (2) 

spin observables, i.e., Pauli matrices in arbitrary directions. The two-qubit observables A\ 
fulfill similar relations as the Pauli matrices, and for separable states it can be shown that 

max{(4 2 >) 2 + (Bf>) 2 } < min{(4 2 >) 2 - (Cf ) 2 } (66) 



holds. This is in fact a necessary and sufficient condition for two-qubit separability 257]. This con- 



struction can then be iterated: One defines four sets of three-qubit observables A? , as A^ 3 ' — 

(2) (2) 

(Ai ® AY =p B\ eg) Br k ')/2 etc., and it is possible to prove similar relations as above for the different 
classes of separability. These conditions can be translated into conditions on the density matrix ele- 
ments by fixing Ai, Bi and GV For instance, for three qubits it can be shown that the matrix elements 
of biseparable states fulfill 

1 018 1 < y/OZZQri + V 033066 + \/044£>55- (67) 



These and the other conditions improve some already known conditions 258]. In addition, they are 
relatively easy to check in experiments. 



'For a different generalization of the CCNR criterion see Ref. 1611 
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(iii) Algorithmic approaches. Also some of the computational criteria can be extended to the multipartite 
scenario. The method of the symmetric extensions can be formulated for the multipartite case (9lj . 
and it can be shown again that the hierarchy is complete. However, this method is practically not very 
feasible. Also, the method of Brandao and Vianna [98| can be applied to show the entanglement of 
some multipartite bound entangled states. 

(iv) Further criteria. There is a series of other a pproa ches that try to extend bipartite criteria to the 
multipa rtite cas e. These include positive maps [259j and conditions on the generalization of the Bloch 

ety of entanglement witness like criteria for full separability, these 



vector 



will be discussed in Section 3.6.3. Finally, it should be noted that any bipartite separability criterion 
can also be used to rule out full separability, since a fully separable state is also separable with respect 
to each partition. 

Finally, it should be added that also criteria for the distillability of a multipartite quantum state are 
known. First, in a similar argument as in Section (3.2) one can conclude that a mixed state is distillable, 
if enough reduced two-qubit states are NPT, allowing the distillation of singlet pairs between them. For 
some special families of states (i.e., stat es that ar e diagonal in a graph state basis or the GHZ basis) more 
elaborate criteria have been worked out 189l 262 1. 



3.6. Entanglement witnesses for multipartite entanglement 

Let us now discuss entanglement witnesses for multipartite entanglement. Clearly, witnesses can be used 
to distinguish the different classes of multipartite entanglement. This is due to the fact that all the classes 
defined above are defined via convex combinations, thus they are convex sets. Contrary to the bipartite case, 
one has several types of witnesses, for the different types of multipartite entanglement. 

3.6.1. Three qubits 

First, there are GHZ class witnesses Wghz that allow to detect mixed states belonging to the GHZ class. 
Thus, their expectation value is positive on all fully separable, biseparable, and W-type states: 

Tr(WcHzg) < =>• g is in the GHZ class. 

Tr(W G Hzg) > => g is not detected. (68) 
An example for such a witness is 

Wghz = |l - \GHZ 3 )(GHZ 3 \. (69) 

This witness is constructed as in Eqs. (28, 29), and the constant 3/4 is the maximal overlap between the 
state \GHZ 3 ) and the pure W class states [185[. In other words, this witness expresses the fact that if the 
fidelity of the GHZ state is larger than 3/4, then the state belongs to the GHZ class. 

In a similar fashion, other witnesses can be defined: Witnesses for genuine tripartite entanglement, denoted 
by W3, have a positive expectation value on all biseparable states thus a negative expectation value indicates 
the presence of true tripartite entanglement. Witnesses for biseparable entanglement W e have a positive 
expectation value on all fully separable states, thus a negative expectation value is still a signature of 
entanglement, which might be only biseparable entanglement. 

Such witnesses can be constructed in a similar manner as above: one makes the ansatz 

W = at- (70) 

where a is the maximal overlap between \ifj) and the biseparable (or fully separable, respectively) states. 
Note that it suffices to consider pure biseparable states, as they are the extremal points of the convex set of 
mixed biseparable states. The calculation of a depends on the type of states, one wishes to exclude. All the 
types of witnesses can be cast into a schematic picture as shown in Figure 5. 

The computat ion of the maximal overlap between a given state \ip) and the pure biseparable states is 
straightforward [l45j : Any pure biseparable state is separable with respect to some fixed bipartition, and 
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Fig. 5. Schematic picture of the set of mixed three-qubit states with three witnesses: The witness Wghz detects GHZ-type 
states, the witness W3 detects three-qubit entanglement, and the witness W e rules out full separability. 

for this case we know already that this overlap is given by the maximal squared Schmidt coefficient [see 
Eq. (28)]. So one can take for the construction of witnesses for genuine tripartite entanglement 

a= max { max {[ s fc(W] 2 }} 

bipartitions bp Schmidt coefficients s k (bp) 

as the maximal squared Schmidt coefficient over all bipartitions. The computation of other overlaps, e.g., 
the overlap with fully separable states is, in general, not so simple. However, for many cases, re sults on this 
problem are known, especially in the context of the so-called geometric measure of entanglement [263l |. which 
will be explained in Section 4.2. To give an example, a witness for genuine multipartite entanglement 
around the state \Wg) would be W3 = 2/3 • 1 — |W3)(Ws| and a witness excluding full separability would be 
W 3 = A/9-l~\W 3 )(W 3 \. 



3.6.2. The general case 

For the case of more than three qubits similar methods can directly be applied. Especially for the detection 
of genuine multipartite entanglement the method from above provides a standard method for the construction 
of witnesses. In addition, the required o verla ps are for many classes of states known. 

For instance, it has been shown in Ref. p64] that the maximal overlap between graph states and biseparable 
states is always 1/2, hence 

W Gjv = I - \G N )(G N \ (72) 

is for any graph state \Gn) a suitable witness. Especially, 

W GH z N = I - \GHZ N )(GHZ N \ and W C l n = | - \CL N ){CL N \ (73) 

are suitable witnesses for GHZ states and cluster states, r espe ctively. For some symmetric Dicke states, the 
maximal overlap can also be computed straightforwardly [213| . e.g., a witness for the \Wn) state reads 

N — 1 

Ww N = 1 - \W N )(W N \. (74) 

and for the \Dn ) state it is 

2 < N 

W* = 5^i-I^><^> (75) 



21 However, such witnesses are not so relevant for this review, as the detection of genuine multipartite entanglement has more 
experimental relevance. 
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Other witnesses for the states in Eqs. (61, 64) read 



W* a = |l-|* a )<*2| and W x = l -l-\ X ){x\- (76) 

Note that these are only the basic constructions for witnesses for multipartite states, more sophisticated 
constructions and their experimental implementation will be discussed in great detail in Section 6. 

3.6.3. Witnesses for full separability 

Let us finally mention some other ways of constructing witnesses (or witness like inequalities) for multi- 
partite systems which result in conditions excluding full separability. 

First, t here are witnesses for multipartite entanglement that are derived in a similar way as Bell ineq ualities 
255l 265]. Then, one can use geo metri cal arguments to construct entanglement witnesses 266l 267 1 or use 



properties of correlation functions [2681 ]. Also, similarly as demonstrated above, one can construct witnesses, 
if the expectation value for a give n observable c an be bounded for product states. To compute such bounds, 
linear programming can be used 



Finally, other ways of constructing witnesses will be discussed in the other parts of this review. Bell 
inequalities will be discussed in Section 5. In Section 6.6 we explain the construction of witnesses using 
stabilizer theory. In Section 8.1 we will discuss separability criteria based on collective angular momentum 
observables (spin squeezing inequalities). In Section 8.3 we will review the possibility to use the Hamiltonian 
of a spin model as an entanglement witness. 

4. Entanglement measures 

In this Section, we explain the basic notions of entanglement measures. We restrict our attention to the 
basic definitions and requirements of entanglement measures that are needed for the methods presented 
later to analyze entanglement in experiments. For a deeper discussion, we can recommend several excellent 
general reviews on entanglement measures 



4.1. General properties of entanglement measures 

4.1.1. Requirements for entanglement measures 

By definition, entanglement measures (or entanglement monotones) should quantify the amount of en- 
tanglement in a given state. As such, a general entanglement measure E{g) should have several desired 
properties [2721 ] . which are listed below. It should be noted, however, that not all of the following properties 
are fulfilled by all discussed entanglement quantifiers. 

(i) First, a natural requirement is that the entanglement measure E{g) vanishes, if g is separable. 

(ii) Second, an entanglement measure should be invariant under a local change of the basis. This means 
that it should be invariant under local unitary transformations, 

E{g) = E(U A ® U B gU\ ® U%). (77) 

(iii) As entanglement can not be created by LOCC, it is reasonable to require that E{g) does not increase 
under such transformations. That is, if A LOCC is a positive map that can be implemented by LOCC, 
then 

E[A LOCC (g)]<E(g). (78) 

Often, this condition is replaced by a different requirement, namely that E{g) should not increase on 
average under LOCC. That is, if a LOCC transformation maps g to states gk with probabilities pu, 
then 

Yp k E(g k ) < E(g). (79) 
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This condition is stronger than Eq. (78), but many entanglement measures also fulfill this stronger 
condition. The monotonicity under LOCC in Eq. (78), implies invariance under local unitary transfor- 
mations. 

(iv) A further property that is often demanded and that is fulfilled by most entanglement measures is 
convexity. That is, one requires that entanglement decreases under mixing of two or more states, 

k k 

This inequality expresses the fact that if one starts from an ensemble Qk of states, and looses information 
about the single instance Qk, then the entanglement decreases. 

Not all measures fulfill this property, and it has been argued that the convexity condition may be 
relax ed b y requiring only that E(g) should not increase if locally distinguishable states are mixed (see 
Ref. [273l | for a discussion). 

(v) Further questions arise, if more than two copies of a state are available. For instance, if Alice and Bob 
share n copies of the same state g, it may be reasonable to require additivity, that is 

E(g® n ) = nE{ e ). (81) 

Even stronger, one may require full additivity. This means, that if Alice and Bob share two different 
states, Q\ and g 2 , then 

E{g l ®g 2 ) = E{g 1 ) + E{g 2 ). (82) 
Such additivity requirements are not fulfilled for some measures, or they are difficult to prove [l9j . 

4.1.2. Convex roof constructions 

There a several strategies to define an entanglement monotone. A first choice is to take the usefulness 
of a state for a certain task as a measure of entanglement. For instance, one may define for a state g the 
optimal distillation rate Ed(q) as a measure of entanglement, as we will discuss in Section 4.2.1. Clearly, such 
quantities are difficult to compute, especially for mixed states, as their computation involves a complicated 
optimization over all distillation protocols. Only for special cases, exact results are known (see Refs. 18, ]j| 
for details). 

Another strategy for defining entanglement measures uses the so-called convex roof construction. For that, 
one first defines the measure E(\i/j)) for pure states, and then defines for mixed states 

E(g) = inf V ft £(W), (83) 

Pk,\<Pk) k 

where the infimum is taken over all possible decompositions of g, i.e., over all pu and \4>k) with g = 
X^fcPfcl^fcX^fcl- I n other words, E(g) is defined as the largest convex function smaller than E(\tp)). 

The advantage of the convex roof construction lies in the fact that the resulting entanglement measure 
has by construction some desirable properties, as E(g) is convex. Moreover, from the properties of E(\ijj)) 
for pure states one can often directly see, whether the measure for mixe d states fulfill s the conditions for 
entanglement monotones, e.g., whether it is non-increasing under LOCC 274 . 275 . 276j |. 



Clearly, the optimization in Eq. (83) is also not straightforward to compute. Again, only for sp ecial case s 
results are known. There are however, some general recipes to give lower bounds on the convex roof [2771 . 278] . 
The main idea is to take an easily computable convex function F(g) and to derive a lower bound in terms of 
F for pure states. Then, as the convex roof is the largest convex function smaller than E for pure states, the 
same lower bound also holds for mixed states. For instance, the trace norm of the partially transposed state 
F{q) = \\q Ta \\i is such a convex function, and from this the lower bound on the concurrence [see Eq. (89)] 
in a d x d-system 



can be established. The resulting lower bounds are often surprisingly good [2771 1278| ]. In the Sections 6.9 



and 7.2 we will investigate how entanglement measures can be evaluated or estimated in experiments. 
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4.2. Examples of entanglement measures 



Now we discuss some examples of entanglement measures. We start with measures for the bipartite case, 
and then explain some measures for the multipartite case. 

4.2.1. Entanglement cost and entanglement of distillation 

The entanglement cost Ec(g) is de fined as the minimal rate of singlets that have to be used to create 
many copies of the state g via LOCC [2791 ] . More formally, it is given by 

E c (g) = inf lim — . (85) 

LOCC iiout — >oo n ou t 

Here, mfLOCC denotes the min imization over all LOCC protocols that map rii n input singlets onto n out 
output copies of the state 



Conversely, the entanglement of distillation Er>{g) is defined as the optimal singlet distillation rate from 
many copies of g, 

E D (g)= sup lim (86) 

where the LOCC protocols map now n- m input copies of g onto n out output singlets. 

With these definitions, Ejj(q) < Ec{g) has to hold. Since there are bound entangled states that require 
entanglement for their generation {Eq > 0), but from which no entanglement c an b e distilled (Ejj = 0), it 
can happen that Eo{g) ^ Ec{g)- For pure states, however, it was shown in Ref. [2791 ] that the entanglement 
cost equals the entanglement of distillation and they are given by the von Neumann entropy of the reduced 
state qa, 

E c (m = = S(g A ) = -Tr[g A log 2 (gA)]. (87) 

Therefore, pure states can be transformed reversibly into singlet states. 

4.2.2. Entanglement of formation 

The entanglement of formation is defined as the convex roof of the von Neumann entropy, 

E F (g)= inf (88) 

Pk,\<Pk) fc 



where (gA)k is Alice's reduced state of the state \ipk), and the optimization is defined as in Eq. (83) 279]. 
Physically, the entanglement of formation may be interpreted as a minimal number of singlets that are 
required to build a single copy of the state. 

A central problem in the study of entanglement is whether or not the entanglement of formation is fully 
additive. This question is not solved yet, but its solution would have wide ranging consequences also on 
other topics, like the additivity of classical capacity of quantum channels (see Ref. [191 ] for a discussion) . 

4.2.3. Concurrence 

A very popular measure for the quantification of bipartite quantum correlations is the concurrence [28i 



2811 ] . This quantity can be defined for pure states as 



C{m = ^2[l -Tr(g A )], (89) 

where gA the reduced state of l^) for Alice. For mixed states this definition is extended via the convex roof 
construction. 



22 To be more precise, one should also include in this definition that the LOCC protocols do not generate prefect copies of the 
state q, but states arbitrary close to it 11911 . 
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The popularity of the concurrence stems from the fact that for two qubits the convex roof can analytically 
be computed [46|. Namely, one has 

C(g) = max{0, A x - A 2 - A 3 - A 4 }, (90) 



where the Ai are the decreasingly ordered eigenvalues of the matrix X = y^/oipy ® < J y)Q*{ (J y ® <J y)\[Q an d 
the complex conjugation g* is taken in the standard basis. Moreover, for two qubits the entanglement of 
formation can be expressed in terms of the concurrence as 



B*)-K( i±£Em y (9i) 

where h(j>) = — plog(p) — (1 — p) log(l — p) is the binary entropy function. For other dimensions, however, 
such a relation does not hold and the physical interpretation of the concurrence is not so clear. Moreover, 
the concurrence is not additive. 

4.2.4. Negativity 

The negativity is just given as the violation of the PPT criterion [52], [H^| 

m = Jl^i, (92) 

where ||...||i denotes the trace norm (i.e., the sum of all singular values) of the partially transposed state. 
Two of the main advantages are that negativity is very easy to compute and it is convex. In order to make 
this quantity additive, one can consider the logarithmic negativity E?f(g) — log 2 \\g TB ||i- This gives an upper 
bound on the ent angl ement of distillation, E^(g) > Ejj(g) [53j ]. The logarithmic negativity is, however, not 
convex anymore [273J . By construction, the negativity fails to recognize entanglement in PPT states. 

4.2.5. Distance measures 

Another class of entanglement quantifiers uses the distance to the separable states as a measure for 
entanglement. As there are several possible notions of a distance, several ent angl ement parameters arise. 
One of the first distance measures is the relative entropy of entanglement [282| , 

E R (g) = MS(g\\a), (93) 

a 

where the infimum is taken over all separable states, and the distance S(f>||cr) is the relative entropy, S(g\\(j) = 
Tr[glog(g) - £>log(cr)]. 

Another class of distance-like measures that are relevant in our future discussion are the so-called robust- 
ness measures. Roughly speaking, these measures quantify, how much noise must be added, in order to make 
the state separable. For instance, the robustness of entanglement R(g) is defined as the minimal s, such that 
the state 

X=^—{Q + so) (94) 
1 + s 

is separable [283| . Here, a denotes an arbitrary separable state. One can s tudy vari ations of this measure, 
by allowing a to be an arbitrary state (the generalized ro bust ness R g (g) [284 285]) or by fixing it to be 
the maximally mixed state (the random robustness R g (g) [283| ). The interesting point for our purpose is, 
that the negative mean value of an entanglement witness can be directly turned into a lower bound on the 
robustness (see Section 6.9). 

4.2.6. The geometric measure of entanglement 

A further interesting en tanglement measure for multipartite systems is the geometric measure of entan- 
glement Eq [2 HI286I |287| . This also quantifies the distance to the separable states. The geometric measure 



is defined for pure states as 

E G {\i,)) = l- sup |<</#}| 2 , (95) 
|*)=|o>|&)|c>... 
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i.e., as one minus the maximal squared overlap with pure fully separable states, and for mixed states via 
the convex roof construction. In many cases, the geometric measure is similarly defined as Eq S (\iP)) = 
—2 log 2 (supi ^ _i n u f,\ | ^ | (01?/;) |). This is then a lower bound on the relative entropy of \ip ), however, it is 



still not additive 



and its convex roof is not an entanglement monotone anymore 



The geometric measure is a multipartite entanglement measure, as it is not only a summation over bipartite 
entanglement properties. Despite of its abstract definition, it has turned out that Eq can be used to quantify 
the distinguishability of multipartite states by local means }290| | . Furtherm ore, the optimization i nvolved in 



263l l29lll or certain graph states [29H29I 



Eq. (95) can be performed in many cases, e.g. for W-type states 



Also the convex roof can be computed in important cases [2 14 . |263j. We will see in Section 6.9 how the 



geometric measure can be estimated from the mean values of entanglement witnesses. 



4.2.7. The three-tangle and comb monot ones 

The three-tangle r, introduced in Ref. [294j . is an entanglement measure for three-qubit states. For states 
in the form of Eq. (44) it is given by 



-(M)=4A*Ai 



(96) 



Interestingly, an arbitrary pure three-qubit state fulfills the monogamy relation C^, BC (\ip}) — C\ b (qab) 



C\ c (qac) 



A\BC\ 

where C^ BC (\ip)) = 2y/det(gA) is the concurrence between A and the other two qubits, 



and C ab (qab) is the concurrence between A and B 294L 295l |. For mixed states, the tangle is defined via 
the convex roof construction, which can, in important cases, be perfor med analy tically 2961 . 1297 1. For more 
than three qubits, the tangle can be generalized via combs and filters 249l 250| . 



5. Bell inequalities 

In this Section, we review Bell inequalities as the oldest tool to detect entanglement. Originally, Bell 
inequalities were designed to rule out local hidden variable (LHV) models. We will therefore first explain 
LHV theories and then present different Bell inequalities for various situations. While Bell inequalities were 
violated in many experiments, several loopholes make it still possible that in principle LHV models could 
describe the measurement results. Wc will discuss the detection efficiency loophole and the locality loophole. 
Finally, we will discuss a two-photon experiment, where the locality loophole has been closed. 



5.1. Bipartite systems 

The first Bell inequality was published by John Bell in 1964 in Ref. 0|. The goal of this paper is to capture 
the Einstcin-Podolsky- Rosen paradox in a quantitative way. The basic idea behind Bell inequalities is the 
following. Let us consider a bipartite system, and perform measurements on both parties. If we assume that 
measurement results existed locally at the parties before the measurement then it is possible to obtain bounds 
on certain quantities composed from two-body correlation terms. In quantum mechanics, it is possible to 
design experiments in which a higher value is measured than these bounds. This, on the one hand, shows 
that quantum physics violates local realism: The measurement is not simply reading out a pre-existing local 
quantity and its results cannot be described by a LHV model. On the other hand, as we will see later, if 
measurements on a quantum state violate a Bell inequality, it implies that the state is entangled. 

5.1.1. Local hidden variable models and CHSH-inequality 

Let us consider a simple example of a bipartite system. We assume that Alice can measure two quantities 
at her party, called Ai and A2, while Bob can also measure two quantities at his party, called B\ and B2, 
as shown in Fig. 6. The results of these experiments are ai, 02 and b\, 62 and we assume that these results 
can take the values +1 or — 1. Now let us perform simultaneous measurements of these four quantities 
M times. Expectation values can be obtained simply by averaging the measurement results, (AiBj) = 
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Fig. 6. The basic setup for the CHSH inequality: in a bipartite system Alice can choose to measure two quantities, denoted by 
Ai and A2, each giving the results a; 6 {—1, +!}• Similarly, Bob can measure B\ and B2 on the right party. 



jf X)fe=i CLi(k)bj(k). If the probabilities of the outcomes are given by expressions like p(a^,b 2 
+1, 62 = — 1) this can be rewritten as 



(A i B j )=p(af,bf)-p(ai,bf)-p(ai 



>bj )+p{a l ,bj ). 



p{ai = 
(97) 



Now one can ask whether these probabilities and expectation values can be described by a local hidden 
variable (LHV) model. Latter assumes that the probabilities of the measurement results are independent 
of whether they are indeed measured or not (the assumption of reality) and that Alice's probabilities do 
not depend on the choice of the observables by Bob and vice versa (locality assumption) . So one makes the 
ansatz 

*>(<*?, ft?) = Jd\p(\)Ax(a?)Bx(b$), (98) 

where a, (3 denote the possible outcomes ±1 and A is the hidden variable, occurring with a probability 
p(A). The response function A\ (af)B\(b^) factorizes due to the locality assumption, as for a fixed A Alice's 
probabilities should not depend on Bob's choice of observables. Furthermore, one can always assume the 
LHV model to be deterministic (i.e., for a fixed A the A\(a™) and B\(bj) take only the values or 1) as a 
nondeterministic model corresponds to a deterministic one, where the hidden variable A is not known (29l.[30j. 
Finally, it should be noted that if the probabilities are of the form as in Eq. (98), then the expectation values 
can be written as (AiBj) = JdX p(A)21a(^)55a(.B,) with %\(Ai) = A\(af) - A\(a~) etc. 

From the fact that the probabilities are described by a LHV model, bounds on the correlations foll ow. 
These are the Bell inequalities. On the level of probabilities, the Clauser-Horne (CH) inequality [298| read j 23 1 

p{a^,b^) + p(a+ } b 2 ) + p(a 2 ,b+) - p{a% ,b%) > 0. (99) 

Eq. (99) reflects the fact that the last event (a^~ : b 2 ) can only occur if at least one of the three other 
possibilities also occurs, and follows directly from the structure of the response function. 

Using the normalization of the probabilities, one can deriv e a si milar inequality for the mean values, which 
is the Clauser-Horne-Shimony-Holt (CHSH) inequality 299l . 300( | . It reads 



{A x Bx) + (A 2 Bi) + (AxB 2 ) - (A 2 B 2 ) < 2. (100) 

Eq. (100) can also be directly understood by considering all possible measurement results 01,02,61,62, as 
for all possible of them ai&i + a\b 2 + 02^1 — 0262 < 2. 

Naturally, inequality Eq. (100) is satisfied if there is a hidden variable model for the probabilities of the 
measurement results. Such an inequality is called Bell inequality. Note at this point that for defining Bell 
inequalities quantum physics is not used at all, only simple basic assumptions lying behind classical physics 
were used. 



'This is the simplest possible form from Ref. [2911 . the original CH inequality is equivalent, but contains six terms. 
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Let us connect these ideas to quantum physics. Assuming that Alice and Bob both have a quantum 
mechanical spin-1/2 particle one can identify Ai and Bj with measurement results of the Pauli spin operators 
X = o~ x and Y — o~ v , respectively. Then, the left hand side of Eq. (100) defines a qu antum mechanical 
observable B = A\ ® B\ + A 2 ® B\ + A\ ® B 2 — A 2 <E> £?2, the so-called Bell operator [30 lj ] that can be 
measured. One finds that there are quantum states that violate the inequality and give 2^/2 for the left- 
hand side. Taking the observables A\ = —cr x ,A 2 = —<r y ,Bi = (a x + <r y )/\/2 and B 2 — (<r x — <T j/)/"V / 2 the 
quantum state with the highest violation is the cigenstate corresponding to the largest eigenvalue of this 
Bell operator. It is given by 

|* CHS H> = ^(|01) - |10», (101) 

which is just the two-qubit singlet state. It can be shown that quantum mechanics does not allow for 
violations larger than 2s/2, a fact which is known as the Tsirelson bound [302| . Furthermore, the question 
whether a give n two-qubit state violates the CHSH inequality for some optimized settings can directly be 
answered ;303j: If one considers the 3x3 matrix r = (<jj ® <jj) with i,j S {x, y, z} then the state violates the 
CHSH inequality iff the sum of the two largest eigenvalues of the matrix r T r is larger than one. For general 
states one can show that any entangled state violates the CHSH ineq ualit y, after local filterings and adding 
an auxiliary system that does not violate the CHSH inequality itself [3041 ] . In general, the problem whether 
a given state violat es violates a Bell inequality for some choice of settings can be tackled using semidefinitc 
programming [305l |. 

How is a violation of the Bell inequality possible? The reason is that the simple assumptions we made 
for the Bell inequality are not valid for quantum mechanics. Indeed, the violation of Bell inequalities in 
experiment shows that either locality or realism must be abandoned, however it is not clear which one of 
them. Interestingly, one can derive Bell-like inequalities that rely only on realism and on some ass ump tions 
about the correlation f uncti ons (AiBj) [3061 ] and are also found to be violated in experiments [3071 ] (for 
generalizations see Ref. [3081 ]). This might be interpreted as if the assumption of realism must be dropped, 
but it should be noted that there are realistic (but non-local) theories such as Bohmian mech anics that 
can describe these experiments. Another interesting non-local realistic model is described in Ref. 309]. It is 
based on simulating the correlations arising from quantum mechanics without communication between the 
parties, using Popescu-Rohrlich non-local boxes [310j . 

5.1.2. Violation of the CHSH inequality and entanglement 

The violation of a Bell inequality implies non- locality, i.e., that the measurement probabilities cannot be 
described by a LHV model. In quantum mechanics, the violation implies that the state is entangled. This 
can be seen as follows. By construction, a separable state in a bipartite system can be written in the form 
g = J2kPkQk ® Bk- Hence, one can write 

(A^) = ^pkTriAiQ^TriBjQ*), (102) 

which is, if we define 2tfc(A;) = Tr(AiQ^) and %5k(Bj) = Tr(Bjg^), directly a LHV model. This means 
that if some measurements on a state violate a Bell inequality then it cannot be separable, i.e., it must be 
entangled. Note, however, that the converse is not true: there are entangled states for which a local realistic 
model for all von Neumann measurements can be written down explicitly (see also below) [44 1 Therefore, 



they do not violate any Bell ine qual ity. The same statement holds even if one allows positive-operator- 
valued measure men t s (POVM s) [31 lj . For pure states, however any bipartite entangled state violates a 
Bell inequ ality [3121 13131 13141 ]. Due to that, Bell inequalities may be viewed as non-optimal entanglement 
witnesses [3151 ] . (Concerning the comparison of Bell inequalities and witnesses see Sec. 6.4.) 

Finally, before considering further Bell inequalities, let us mention ways to characterize them from the 
point of view of their usefulness in experiments. We have seen that the CHSH inequality is maximally 



24 The question whether a state allows for a LHV model or not can also be tackled by looking for symmetric extensions of the 
state using semidefinite programming 102], 



35 



violated by the state Eq. (101). However, in a real experiment one can never prepare the state Eq. (101) 
perfectly, it is always mixed with noise 

ftioisysingletb) := p\ *CHSh) (^CHSH | + (1 -p) J, (103) 



which is a two-qubit Werner state 44 1. The CHSH inequality is violated if p > l/v2 « 0.71 and, using 



the PPT criterion, the state is entangled for p > 1/3. Moreover, it is kn own that this state allows for a 
LHV model (for all von Neumann measureme nts) for p < 0.6595 44, 316j and it violates a Bell inequality 
(different from the CHSH one) for p > 0.7056 [317|f^1 

The bound for the CHSH inequality characterizes the minimum visibility that the CHSH inequality re- 
quires. In general, the visibility is the ratio between the maximum of the Bell operator B for quantum states 
and the maximum for LHV models 

;= max^H (104) 
maxLHV \ {d)lhv\ 

Clearly, a large visibility is an advantage from the experimental point of view. 



5.2. Multipartite systems 



5.2.1. Multipartite LHV models 

In the previous subsection we discussed Bell inequalities for bipartite systems. The ideas can straightfor- 
wardly be generalized to multipartite systems. Let us consider an TV-partite system, in which each party can 
measure some observables. Then the measurement results fit a local hidden variable model if for any choice 
of measured variables we have 

A.ll.CJh... = JdXpiX^xiA^xiBjXxiCk^xiDt)..., (105) 

where — 1 < %\(Ai) < +1 etc. If a multipartite Bell inequality is violated then there is no local hidden 
variable model that could explain the measurement results. 

Similarly to the definition of biseparability one can also ask, whether i t is possible to explain our mea- 
surement results by bipartite models in which several of the parties join 173l 318l 319} . For example, for 
the three-partite case this would mean that any expectation values could be explained by 



(A i B j C k )=p 1 Jd\p(\)Vl x (A l )X x (B J C k )+p 2 ydAp'(A) ( 8 A (S,)2) A (A l C fe )+P3 J d\p"(\)£ x (C k )3x(AiB j ), 

(106) 

where pi + P2 + P3 = 1- If there is not such a description for the measurement results, we say that there is 
genuine three-partite non-locality in the system. 



5.2.2. Mermin and Ardehali inequalities 

In the previous subsection we discussed the CHSH inequality for bipartite systems. It needed the mea- 
surement of two quantities for each party and these measurements had two outcomes. Simila r two -setting 
inequalities exist also for multipartite systems. The most famous one is the Mermin inequality [206| . For N 
qubits it is given by 

Y J {M 2 X 3 X i X^--X N )-Y J {yiY 2 X^X A X^--X N )+Y J ^y2Y 3 Y i X^--X N ) + < L M ormin, (107) 

7T 7T TV 

where and Yi represent the Pauli matrices a Xl a y on the z-th qubil l 26 1 . represents the sum of all 
possible permutations of the qubits that give distinct terms and ^Mermin is the maximum for local states. 
It is i M crmin = ^ N/2 for even N and L Mcrmin = 2 ( ^ N - 1 ^ 2 for odd N. 

25 See also Problem 19 on http://www.imaph.tu-bs.de/qi/problems/. 

26 Of course, in the sense of LHV models, they can be replaced by arbitrary observables with two outcomes ±1. 
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The Mermin inequality is maximally violated by the GHZ state defined in Eq. (54 j 27 1 . The Bell operator 
is the sum of 2( A '~ 1 ) terms that all give +1 for the GHZ state. Hence the maximum of the Bell operator 
is 2' JV ~ 1 ) for quantum states. Thus, the visibility is VMermin = 2 JV//2 ~ 1 for even N and 2 ( - N ~ 1 ^ 2 for odd N. 
Note that the visibility increases exponentially with the number of qubits. It can be proved that for odd 
number of parties the Mermin inequality has the largest visibility possible among two-setting inequalities 
with full correlation terms. Latter means that the Bell operator is the sum of correlation terms that contain 
a variable for each spin. For example, for N — 4, X\Y 2 X 3 Y 4 is a full correlation term, however, X1Y2I3X4 
is not. The I3 should indicate that on the third particle no measurement has been done. 

The Mermin inequality consists of stabilizing operators of the GHZ state (see Sec. 3.4.3). This implies that 
there is a quantum state that gives +1 for al l ter ms, while LHV models have to result in —1 for at least one 
term. This is a form of the GHZ argument [l T4j | . which leads to an obvi ous contradiction to local re alism. 
Similarly, one can derive Mermin-like inequal ities for other graph states 147 . 228 . 229l l32d( . 1 3 2 lk l322| . 



Let us now discuss the Ardehali inequality [3231 ] . For an even number of parties this inequality is superior 
to the Mermin inequality. It is defined as 

((4 +) - A[- ] )( - J2 W 3 X 4 X 5 --Xn + Y, Y 2 Y 3 X 4 X 5 ■■Xn-Y, Y 2 Y 3 Y 4 Y 5 X e ■ -X N + ... - ...)) 

TV TV TV 

+ ((A[ +) + 4 _) )(S y 2 ^3^4X5 --Xn-Y 12*3*4*6 --X N + Y X 2 Y 3 Y 4 Y 5 Y 6 X 7 ■ -X N - ... + ...)) 

TV TV TV 

< -^Ardehali, (108) 

where A± are operators corresponding to measuring the first spin along directions corresponding to the 
quantum operators A± = (=F*i — *i)/\/2- Note that on the other spins X or Y is measured as in the case 
of the Mermin inequality. The constant £ Ardehali appearing in Eq. (108) is the maximum for l ocal states. It 



is 2 N I 2 for even N and 2 {N+1 ^ 2 for odd N. It is again maximally violated by the GHZ stated The Bell 
operator is the sum of 2 N terms that all give l/y/2 for this state. Thus, the visibility VMermin = 2 N ' 2 for 
even N and 2^ N ~ 1 ^ 2 for odd N. 

An interesting feature of the Ardehali inequality is that the Bell operator is essentially the same as for the 
Mermin inequality, but the rewriting with A^ leads to a higher violation of local realism. This shows that 
the violation of local realism of GHZ states with an even number of q ubits can increase, if non-stabilizer 
observables are considered. This also holds for some other graph states 320l 



Note that so far all the inequalities needed the measurement of two operators for each qubit. Also, they 
were the sum of full correlations terms. Among such inequalities, it can be sh own that for any N the Mermin- 
Ardehali construction, also discovered independently by Klyshko 325, 326] provides Bell inequalities with 



an optimal violation. It turns out that the full set of such Bell inequalities can be generated and written 
down concisely in the form of a single nonlinear inequality [3271 l328j . Sur prisingly, there are multi-qubit 
pure entangled states that do not violate any of these Bell inequalities [329j . 

The situation is more complicated with inequalities that are not sums of full correla tion terms. Then, it has 
been shown that such inequalities can detect any pure entangled multi-qubit state |330| . Also, inequalities 
of this type c an be constructed such that they are maximally violated by cluster states and graph states 
fl47l [228LI320I ] (see Sec. 3.4.3 for the definition of these states). In particular, for the four-qubit cluster state 
this inequality looks like [2281 13311 j 



(X^XsZi) + (Z 1 Y 2 Y 3 Z 4 ) + (X 1 1 2 Y 3 Y 4 ) - {Z 1 Y 2 X 3 Y i ) < 2. (109) 

Note that on all of the qubits two operators are measured except for the second qubit for which only Y 2 is 
measured. One can show that for a large class of graph states, e.g., for linear cluster states, i t is p ossible to 
construct two-setting Bell inequalities that have a visibility increasing exponentially with N [320] . 



27 Note that the original Mermin inequality is slightly different and is maximally violated by the state (|000...) + i|lll...))/%/2. 

28 Note that the original Ardehali inequality has a different definition for ' and is maximally violated by the state ( |000. . .} — 
|lll...»/^2. 

29 This is the inequality for the cluster state in the basis as defined via Eqs. (56,57), for the cluster state as in Eq. (60) it reads 
(Z1I2X3X4) + (X 1 Y 2 Y 3 X 4 ) - (Zilal^n) + (X 1 Y 2 X 3 Y i ) < 2. 
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5.2.3. Bell inequalities detecting genuine multipartite non-locality 

So far all the Bell inequalities presented could be used to confirm that a state is nonlocal, that is they 
ruled out local hidden variable models as in Eq. (105), but did not give any information on whether or not 
the state possesses genuine multipartite non-locality [see Eq. (106)]. However, some Bell inequalities exist 
that can rule out such hybrid LHV models. Let us consider the three-particle Svetlichny inequality as an 
example |l73| . It reads 

(A^id) + (A 1 B 2 C 1 ) + (A 2 B 1 C 1 ) - (A 2 B 2 C 2 ) + (A 2 B 2 d) + (A 2 B 1 C 2 ) + (A X B 2 C 2 ) - (A&d) < 4, (110) 

where the Ai,Bj and C'k are two outcome measurements (e.g., Pauli matrices) on each particle. This Bell 
inequality is a sum of two three-particle Mermin inequalities. However, one can directly see that for hybrid 
LHV models (e.g., for an 1 1 23 case, where the particles 2 and 3 are combined, the Svetlichny inequality can 
be written as a sum of two CHSH inequalities) the same bound holds. For GHZ states, the Bell operator 
can take values of Ay/2 |319l . |332| . however, not all genu ine m ultipartite entangled states lead to a violation 
of the Svetlichny inequality [256j . even if they are pure [333| . 

A generalization of this inequality for an arbitrary number of parties is presented in Ref. |318| . Ref. (3 19( | 
gives bounds for Bell inequalities both for genuine multi-qubit non-locality and genuine multi-qubit entan- 
glement. 



5.2.4. Further Bell inequalities 

There are Bell inequalities that use more than two operators for each particle (i.e. , more than two settings) 
or need the measurement of observables with more than two outcomes. Ref. [334j considers bipartite sys- 
tems with qudits having dimension d > 2, and studies two-setting Bell inequalities with observables with d 
outcomes. It fi nds t hat with the increase of d, it is possible to find Bell inequalities with increasing maximal 



violation. Ref. [3351 ] presents a family of Bell inequalities, called the Collins-Gisin-Linden-Massar-Popescu 
(CGLMP) inequalities for bipartite quantum systems of arbitrarily high dimensionality, which are strongly 
resistant to noise. They a re fo r systems with two measurement settings per site, with more than two mea- 



surement outcomes. Ref. 336| presents Bell inequalities stronger than the CHSH inequality for isotropic 
states of two three-state particles. 

Systematic study of such inequalities were also presented for small number of outcomes. Such an inequality 
can be characterized by four numbers ijmn, where Alice measures one of i operators that have m outcomes. 
Similarly, Bob measures one of j operators that have n outcomes. In the bipartite case, for two t wo-ou tcome 
observables the CHSH inequality is the only Bell inequality apart from simple transformations [3371 ]. That 
is, it is the only tight Bell inequality, whe re tight means that the inequality is a facet of the polytope of 
correlations allowed by LHV models [3381 ]. Non-tight inequalities also exist, but they do not detect more 
states than the tight inequalities. 

For the 3322 case, th at is, when thr ee two-outcome observables are measured at each party, then there is 
a single new inequality [339l l340l . l34lj . It reads 

(A 1 l)-{A2l) + {lB 1 )-{lB2)-{A 1 B 1 )+{A 1 B2)+{A2B 1 )-{A 2 B2) + {^3) + {^B 3 )+{A 3 B 1 )+{A 3 B2) < 4. 

(Ill) 

This is a re levan t inequality in the sense that it detects two-qubit state s tha t cannot be detected by the CHSH 
inequality [340| . For the 4322 case there are three new inequalities |340l |. For systems with measurements 
with more outcom es th e num erica l search is exceedingly difficult. A set of 26 inequalities for the 4422 case are 
presented in Ref. [342j |. Ref. [343| determines the quantum violation of 241 tight bipartite Bell inequalities 
with up to five two-outcome measurement settings per party, for up to eight dimensional complex and six 
dimensional real Hilbert spaces, using semidefinite programming. The violation of bipartite inequalities with 
arbitrary number of two-outcome observables is bounded by the Grothendieck constant [3161 ] . 

For th ree p arties and for two-outcome observables, inequalities with arbitrarily large violation can be 



designed [3441 ] . Ref. [3451 ] presents three-setting two-outcome inequalities for GHZ states that have a higher 



violation than the Mermin inequalit y. A set of multipartite thre e-setting Bell inequalities were studied, with 



two-outcome observables in Refs. [346l 13471 ] . Finally, Ref. [3391 ] presents the full list of Bell inequalities for 



three parties, two measurements and two outcomes. 
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5.3. Consequences of a Bell inequality violation 



We have seen already that for the bipartite case a violation of a Bell inequality implies entanglement. For 
the multipartite case similarly non-locality implies entanglement, however, only ruling out LHV models as 
in Eq. (106) can prove genuine multipartite entanglement. As in the bipartite case, entanglement does not 
imply non-locality. For example, there are three-qubit state s with genuine multipartite entanglement that 
are local if von Neumann measurements are considered |348l | . 

Interestingly, while Bell inequalities are primarily designed to rule out LHV models (as described in 
Sec. 5.2.1) they can sometimes be used to detect genuine multipartit e entangl ement. For example, the 
three-qubit Mermin inequality used in the photonic experiments of Ref. 210l 349j is 



(X 1 X 2 X 3 ) - (Y 1 Y 2 X 3 ) - (X 1 Y 2 Y 3 ) - (ftXiYa) < 2. (112) 

As a Bell inequality, violation of this inequality does not imply three-qubit entanglement. In quantum 
mechanics, however, it can be proved that, if Xk and Yk refer to t he measur ement of Pauli spin matrices 



at qubit k, then for biseparable states the maximum is also two [191L 13501 ] . T hus, the violation of this 
inequality implies genuine multi-qubit entanglement and in the experiment of Ref. [210l ] g enuin e multi-qubit 
entanglement was detected ! 30 1 In another experiment creating a four-qubit GHZ state 211 j, a bound for 



multipartite entanglement for the Ardehali inequality was used to detect four-qubit entanglement. In general , 
the detection of multi-qubit entanglement with many-qubit Bell inequalities was studied in Re fs. lim [352^ . 



These papers presented bounds for various forms of multipartite entanglement. Refs. 255 . 353| present 
nonlinear Bell inequalities for the same aim. 

Besides, there are some other interesting connections between the violation of a Bell inequality and the 
properties of a state: 

(i) Quantum communication. It has been observed that in the BB84 protocol the error thresholds where a 
Bell inequality can be viola ted and where the mutual information between Alice and Bob and Alice and 



Eve a re th e same, coincide [354J . 13551 . 13561 ] . Similar results have been found in the multipartite scenario 
357 . |358{ . This suggests a connection between Bell inequality violation and security of quantum 
cryptography. Moreover, the Ekert protocol is based on distributing entangled particles [J]. Its security 
is directly related to the violation of a Bell inequality. Note that Bell inequality violation can imply 
security even in the device independent case when no assum ption is made on the quantum system we 
are using for communication (e.g., dimension of the system) |359| . 
(ii) Communication complexity. The violation of a Bell inequality means that measurement results on a 
quantum state cannot be replaced by shared classical randomness. Because of that, using the quantum 
state as a resource has an advantage in quantum communication. Indeed, for any state violating a 
Bell inequality one can write down a communication complexi ty problem , for which the quantum state 
improves the solution compared with classical resources [360l . I361L 362]. 



(iii) Distillability. Quantum states violating two-setting Bell inequalities are bipartite distillable 3631 . 1364 
This fact, however, relies on the spectral decomposition of the Bell operator, and not on the nonex- 
istence of a LHV model. In the multipartite case, this means distillability at least with respect to 
some bipartitioning of the qubits, which arises if two groups of parties join. On the other hand, if the 
distillation protocol has to be local wi th r espect to each party, there exist multipartite non-distillable 
states which violate a Bell inequality |365| . 

This directly leads to the conjecture of Peres, stating that bipartite states with a positive partial 
transposition allow a LHV description. While no proof or counterexample has been found, in several 
special cases it is shown that such states do not violate large classes of Bell inequali ties, e.g., for the 
Mermin-Klyshko inequalities [3661 ] and for a wide range of multi-setting inequalities |367l ] . 



This has also been shown in Ref. l35lll based on a different argument. 
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5.4. Loopholes 



While Bell inequalities have been violated many times in experiments, several loopholes make it still 
possible in principle that the experimental results can be explained with a LHV model. There are mainly 
two loopholes: 

(i) The fair sampling assumption or detector efficiency loophole. In many setups based on photons, photo- 
detectors detect photons only in part of the experiments, while those experiments in which there are 
no photons detected (no "clicks") go unnoticed. We assume, naturally, that the unnoticed experiments 
have the same statistical properties as the noticed ones. However, a so far unknown mechanism could 
cheat us in this respect. It could happen, that all the experiments tog ether wou ld not violate a Bell 
inequality, however, the subset for which our detectors click, violate it [368l l369j . 

The question, for which detector efficiency the loophole becomes closed depends on the Bell inequal- 
ity, and is, in general, not straight forward to answer. For several important experimental situations, 
results have been obtained in Refs. [373, [373, HzE HzE . The detection loophole was eliminated in 
experiments with trapped ions with very large detection efficienc y in R ef. |375| . In a recent experiment, 
also eliminating the detection loophole, the ions were lm apart |376| . 

(ii) The locality loophole. Ideally, before each correlation measurement the circuits connected to the de- 
tectors should decide locally which operator to measure. This has to be done such that the detectors 
cannot communicate their choice if we assume that the speed of this communication cannot be larger 
than the speed of light. In practice this means that the event of deciding the measurement direction 
for one of the detectors must be space-like separated from event of re ading out the measurement on 
the other detector. This requirement is very challenging. In Ref. |377l |. an experiment is desc ribed in 
which the directions for detection are determined after the photons left the source. In Ref. [3781 ] . a 
photonic experiment is described that eliminates the locality loophole, however, in this realization the 
detection loophole is still present. We will discuss it in the following section. 

Another assumption in the derivation of Bell inequalities is that Alice's and Bob's choices of the mea- 
surements do not depend on the h idden variables. This corresponding loophole is sometimes called the 
freedom-of-choice- loophole l38dl ] and has been addressed experimentally in Ref. [380j |. 



5.5. Experimental violation of a Bell inequality 

In this Section we des cribe as a first experiment the violation of Bell inequalities in a two-photon experi- 
ment by G. Weihs et al. [3781 ]. 

5.5.1. Generation of entangled photons 

First, we review briefly the fundamentals of the experimental techniques, since there are numerous similar 
many-photon experiments. An excellent and exhaustive review of this topic can be found in Ref. [12]. As is 
well known, a photon can be used to store a single qubit information since it can be in a superposition of 
the horizontally and vertically polarized states, denoted by \H) and \V), respectively, as 

\*) = a\H)+0\V), (113) 

where |a| 2 + |/3| 2 = 1. Other basis states can be also chosen. The left and right circularly polarized states 

can be expressed as \L) := -75G-H) + i\V)) and \R) := ^(\H) - i\V)). 

The basic linear optical elements for the manipulation of a single photon are as follows: 
(i) Wave plates are devices that delay one of the polarization components with respect to the other, and 
from the state Eq. (113) they lead to 

I*) = a\H) +e**j0|V), (114) 

where 4> describes the delay. For half-wave plates (HWP) and quarter-wave plates (QWP) the delays 
are 4> — n = 180° and <f> = tt/2 = 90°, respectively. 
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(ii) A beam splitter (BS) is a device that has two spatially separated input modes. The state of the two 
spatially separated output modes can be obtained through linear transformations. That is, let us denote 
the two spatially separated modes {\T) and \B)} referring to Top and Bottom. Then a 50 — 50% beam 
splitter that does not introduce any additional phase shift makes the transformation 



(iii) Polarizing beam splitters. If photons carry also a polarization degree of freedom, then a polarizing beam 
splitter (PBS) can be used to direct photons with H polarization into one direction, while photons with 
V polarization into the other direction. Using of polarizing beam splitters is typical before detecting 
the state of the photons. That is, after the polarizing beam splitter two detectors detect the photons 
of the two different polarizations. 
Next we discuss, how to realize multi-photon experiments. The difficulty of using photons, rather than 
atoms or ions, for creating entangled states is that they do not interact with each other. Thus, special tech- 
niques are needed to create entangled states. The most often used method is using spontaneous parametric 
down-conversion and post-selection. 

(i) Spontaneous parametric down- conversion (SPDC) is a process in w hich a photon with frequency v 
is converted to two photons with lower frequencies v\ and vi |38lj . Due to conservation of energy, 
v = V\ + 1>2 . In a typical experiment, by filtering out the unwanted frequencies one sets the frequencies 
of the outgoing photons equal, that is, v\ — = |- Due to conservation of momentum, the sum of the 
momenta of the outgoing photons equal the momentum of the incoming photon. The down-conversion 
is collinear, if the outgoing photons proceed in the same direction as the incoming one, while it is 
non-collinear if this is not the case. Non-collinear down-conversion allows the spatial separation of the 
two photons right after the down-conversion. In case of type-I down-conversion the two photons have 
the same polarization, while for type-II down-conversion they have perpendicular polarizations. 

(ii) Conditional detection or post-selection. Since the SPDC happens with some small probability, e.g., 
when a laser beam is directed to a /3-Barium-Borate (BBO) crystal, and the photon detectors do not 
have perfect efficiency, only the events where two photons are detected are relevant. In experiments 
with more than two photons, the photons are usually generated by higher order emissions of SPDC 
(see Sec. 6.5) or by interaction of entangled photon pairs via beam splitters. Then, it happens that the 
N photons are not equally distributed over the N spatial modes, i.e., with some probability one mode 
may contain more than one photon and some modes may contain none. In this case, the observation 
of a multipartite state is usually based on conditional detection (or post-selection) : only if each spatial 
mode contains one photon, i.e., all detectors register one photon, the desired state is observed (for an 
example see Sec. 6.5). 

5.5.2. The experiment 



Now we are in the position to review the experiment of Ref . [3781 ] . In the experiment they created a two- 
qubit singlet state [Eq. (101)] and used a CHSH Bell's inequality [Eq. (100)] to prove the non-locality of the 
state. The key point of the experiment is that the two detection events are space-like separated. Thus, what 
is measured at one of the qubits cannot influence the results of the measurement of the other qubits, even 
if we allow for a physical mechanism, so far unknown, that could make communication possible between 
the two detectors with the speed of light. Thus, this experiment closed the locality loophole discussed in 
Sec. 5.4. Note, however, that it did not close the detection loophole since the detection efficiency for photonic 
experiments is relatively low (in the present case it was 5%). 

In the experi ment , type-II parametric down-conversion was used to generate an entangled photon pair in 
state Eq. (101) [38 1| . Then the two photons were led to the two observers stations (see Fig. 7) with optical 
fibers. In order to achieve space- like separated detector events ("clicks"), the observer stations were 400m 
away from each other. 



\T)^^(\T) + \B)), 
\B)^M\T)-\B)). 



(115) 
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Fig. 7. One of the two observer stations. All alignments 
and adjustments were pure local operations that did not 
rely on a common source or on com munication between 
the observers. Figure from Ref. |378ll . 
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Fig. 8. Four out of sixteen coincidence rates between var- 
ious detection channels as functions of bias voltage (ana- 
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tor with switch having been in position "1" and Bob's 
"— " detector with switch position "0" . The difference in 
height can be explained b y dif ferent efficiencies of the de- 
tectors. Figure from Ref. |378|| . 



It was decided locally, which one of the two operators to measure at each qubit. This was done using 
a quantum random generator, rather than a pseudo-random number generator. Latter would be based on 
a deterministic process and in principle one party could predict what the other measures knowing the 
algorithm of the generator. Moreover, the timing of the generation of the random decision for the detectors 
was such that, again, communication that is not faster than light, could not bring information on the choice 
of what to measure at one detector to the photon at the other detector. This means that the generation of a 
random bit, the choice between the two possible measurements and even detection of the result was realized 
locally at the two parties, as can be seen in Fig. 7. The operator to be measured was set by an electro-optical 
modulator. The result was recorded in a logic circuit together with a time tag and the information on which 
one of the two operators were measured. 

If the two polarization analyzers are set into directions a and /3, respectively, then the coincidence rates 
show a sinusoidal dependence on a — (3, which can be seen in Fig. 8. The solid lines are the theoretical curves, 
while the dots correspond to the measurement values. The experimental results fit very well the theoretical 
calculations. The mean value of the Bell operator was 2.73 ± 0.02 and the visibility is 97%. For that 14700 
coincidence events were collected in 10s. Note that the coincidence rates in multi-qubit experiments based 
on down-conversion and post-selection are much lower, and the experimental time is much longer. Such 
experiments are presented in Sec. 6.5 and Sec. 8.2. 

6. Entanglement witnesses 

In this Section, we will explain how entanglement witnesses can be used for the experimental detection of 
entanglement. The theoretical concept of entanglement witnesses has already been introduced in the Sec- 
tions 2.5 and 3.6, and some major construction methods were shown. Now we will discuss how witnesses can 
be implemented in experiments. Furthermore, we will describe advanced methods to construct witnesses for 
special experimental situations and will also discuss several recent experiments from the point of view of 
entanglement detection. Finally, we will discuss how entanglement witnesses can be used to derive quanti- 
tative statements about the entanglement in a quantum state and will explain applications of witnesses in 
quantum key distribution. 
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6.1. Witnesses for two qubits 



Let us start our discussion with a simple two-qubit example 148l . 382]. From an experimental point of 
view, the investigation of two qubits seems trivial, as most experiments nowadays aim at the generation of 
multiparticle entangled states. However, the investigation of the two-qubit case will show already the main 
features of entanglement detection via witnesses in a simple setting. 



6.1.1. Construction of the witness 

Let us consider a setup that is intended to produce a particular pure state \ip){ip\, but due to the im- 
perfections some noise is added. As a simple approximation, one may consider a mixed state g of the form 



g(p) :=p\il>){il>\ + (l-p)~, 



(116) 



where 1 /4 is the maximally mixed state of two qubits. The aim is to give a scheme to decide whether g(p) is 
entangled or not. Clearly, as the form of the state in Eq. (116) is only a rough approximation, the conclusion 
that the state is entangled should not depend on any assumption on the form of the state. However, it is 
reasonable to choose the detection method such that it detects states of the form Eq. (116) even for small 
p. The optimal entanglement witness for an entangled g(p) is easy to construct from the results of Section 
2.5. First, one has to compute the eigenvector corresponding to the negative eigenvalue of g(p) Tn , then the 
witness is given by the partially transposed projector onto this eigenvector. 

If the Schmidt decomposition of \ip) is \tp) — a\01) + 6 1 10) the eigenvalues of g(p) TB are given by {(1 — 
p)/4 + pa 2 : (1 ~p)/4:+pb 2 - 1 (1 -p)/A+pab] (1 - p)/A - pab}. Therefore g(p) is entangled iff p > 1/(1 + Aab). 
The eigenvector corresponding to the minimal eigenvalue A_ is given by 



|0_) = i(|OO)-|ll)), 
and thus the witness W can be written in several ways as 



(117) 



W = \4>-){<t>- 
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(118) 



with \ip+) — ( 1 1) + |10))/a/2. Note that this witness does neither depend on p, nor on the Schmidt coeffi- 
cients a, b. Due to the special assumption about the state, it detects g(p) iff it is entangled, since we have 
Tr[\^)(^\ T -g(p)} = Tr[\^)(^\g(p) T -} = A_. 



6.1.2. Local decomposition of the witness 

In the form of Eq. (118) the witness W cannot easily be implemented in an experiment. The reason 
is that W is up to a constant a projector onto one of the Bell states and such a nonlocal obs ervable 
is not straightforward to measure in practice (for a theoretical proposal to do that see Ref. |383| ). The 
observables that can easily be measured in any experiment are local observables like (a z ® er 2 ) or projectors 
like (|0) (0| <£> |0)(0|). Therefore, for the experimen tal im plementation it is necessary to decompose the witness 
into operators that can be measured locally 

muni- Thus, we need a decomposition into projectors onto 

product vectors of the form 



W = £)c i |e i )<c,| »!/,)(/, 



(119) 



Such a decomposition can be measured locally: Alice and Bob measure the expectation values of the |e<) (e^ | (g> 
\fi)(fi\ and add their results with the weights Cj. One can construct such a decomposition in many ways, 
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but it is reasonable to do it in a way that corresponds to the smallest possible experimental effort for Alice 
and Bob. 

The experimental "cost" that Alice and Bob have to pay when measuring the witness W via such a decom- 
position is the number of measurements they have to perform. As already mentioned, local measurements 
are easily possible, and by measurements we always mean von Neumann (or projective) measurements. One 
measurement on Alice's side in this sense consists of a choice of one orthonormal basis for Alice's Hilbert 
space. For qubits, this corresponds to a polarization measurement in a certain direction. For a particle with 
spin s one may interpret this as the choice of a direction for a Stern-Gerlach-like apparatus. With this setup, 
Alice can measure the probabilities for all the outcomes. 

A similar argument holds for Bob, thus, if one considers a d x d system a term of the form 
d 

M=J2 c u \a k )(a k \®\bi){bi\ (120) 

k,l=l 

with (a s \at) = (b s \b t ) = 5 st can be measured with one collective setting of measurement devices of Alice 
and Bob. Alice and Bob can discriminate between the states \akh), measure the probabilities of these states 
and add their results with the weights c k i using one collective setting and some classical communication. 
See also Fig. 12 for a four-qubit example of such measurement results. 
It is therefore reasonable to find a decomposition of the form 

m m d 

w^^ = E E 4*14X41 ®W)( b i\ (121) 

i=l i=l k,l=l 

with (a l s \a\) — (b\\b\) — S st and minimal number of measurements, i.e., a minimal m. Note that a decom- 
position of the form (121) is more general than a decomposition into a sum of tensor products of operators 
W = YITLi Ji^i <8 £>i. The difference is that a measurement as in Eq. (120) consist of a choice of observables 
(say, a z (g> <t z for definiteness), but then all expectation values like (a z <g> <r z ), (a z <g> 1) and (1 <g> a z ) can 
be de term ined from the same data. These terms would occur as different Ai % Bi in a decomposition like 
above Hj] 

The optimal decomposition of the witness can be obtained as follows: W is a partially transposed projector 
\tfj}(tp\ TB , and we consider the Schmidt decomposition of — a|00) + /3|11). The actual W is the special 
case a = l/\/2 = — /?, but one can directly deal with the most general We define as usual the spin 
directions by \z + ) = |0), \z~) = |1), l^*) ~ ^(lO) ^ |1); Iz/*) = 7^(1^) i *|1) an d obtain the decomposition 

\ip)(tp\ TB =a 2 \z+z + ){z + z+\ + f3 2 \z- Z -)(z-z-\ +a/3 (\x+x+) (x + x+ 1 + \x~x~) (x~x~ \- 

-\y + y-)(y + y-\ - \y-y + )(y-y + \) (122) 

= - [l <E) 1 + a z <E> g z + (a 2 - (3 2 ){a z ® 1 + 1 <g> a z ) + 2aj3(a x ®a x + a y ® a y )] . (123) 

This decomposition into six product vectors requires only a measurement of three settings: Alice and Bob 
only have to measure the coincidence probabilities of a x <g> <j x , a y <g> a y , and <r z (g) a z to measure \ip)(ip\ TB . 
We will discuss an implementation of this witness in Section 6.2. 

The question a rises , wh ether this decomposition is really optimal. Indeed, this is the case. It can be 
shown (see Refs. 148, 382]) that the witness W can not be decomposed in two measurements, hence the 



decomposition is optimal. 

Finally, it should be noted that for the case that the target state is given by \tp) — (|01) + |10))/\/2 = 
\tp + ), the witness W also allows for a determination of the fidelity of the state. This quantity is defined as 
Ftp = {^\qW) an d quantifies to which extent the desired state was produced. Due to Eq. (118) we have 

F^ = \- (W) (124) 



31 One can also ask for a decomposition of W that minimizes the number of produ ct vectors in Eq. (119) 1471 . although this is 
not optimal from an experimental point of view. This has been discussed in Refs. ll4Stl383l . 
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and the decomposition of W in three local measurements just means that the fidelity of can be 

determined with three measurements. This connection between mean values of witnesses and the fidelity of 
the target state will also occur in the multipartite setting. 



6.1.3. Discussion 

We have seen that with the help of a local decomposition entanglement witnesses can be used for the 
experimental detection of entanglement. For two qubits, this required only three local measurements, which 
is less than quantum state tomography that would require 3 2 = 9 measurements, namely the measurement 
of all correlations (<Ji <g> <jj). In general, two criteria determine the efficiency of an entanglement detection 
scheme: 

(i) Robustness to noise: The robustness to noise is given by the minimal p in Eq. (116) that has to be 
reached for the state to be detected. Equivalently, one can define a p no isc = 1 — p and ask for the 
maximal p n oise that can be tolerated. This robustness quantifies the noise that can be allowed in an 
experiment, and if some amount of noise is expected one can estimate, whether there is a chance to 
confirm the success of the experiment at all. 

Clearly, in an experimental setting the noise does not have to be white and the state may not be 
of the form as in Eq. (116). However, this is often a reasonable approximation, hence the robustness 
to noise is a good indicator to compare different entanglement witnesses or different entanglement 
detection methods. 

(ii) Number of measurement settings: Another crucial point is the number of local measurement settings 
that are required for an implementation. To use as few measurements as possible is important from 
several points of view: First, in multi-photon experiments the count rates (i.e., the number of suc- 
cessfully generated states per time unit) is often low. Therefore, the measurement of a single setting 
with a given accuracy needs a certain time, and only a restricted set of measurements can be done 
while keeping the setup stable. But even if measurements can be repeated fast, and a high rate of 
states is available (as it is in ion trap experiments), state tomography would require an exponentially 
increasing effort, making it practically impossible for ten or more qubits Finally, for experiments 
using hyper-entangled photon states, the switching of one measurement setting to the other requires 
the adjustment of an interferometer, which needs a considerable amount of time and work 234, 235] 
(see also Section 6.7). 

In the example above, three measurement settings suffice, which is a modest improvement com- 
pared with full state tomography. However, the difference becomes large, if the number of particles 
increases. For tomography of an A-qubit state, 3^ measurements are necessary, which amounts to an 
exponentially growing effort. For the measurement of the witness, typically only a few measurements 
are required. As we will see in Section 6.6, for important states, even for arbitrarily many qubits, only 
two measurements are required for the detection of entanglement and an estimation of the fidelity. 



6.1.4. Further results for bipartite systems 

Let us finally mention further results on the decomposition of witnesses for bipartite systems. As a first 
question, one may ask what happens if the witness is of the form W = (i/j\ Tb , where 

M=X}Ai|it) ( 125 ) 

i 

is a state on a d x d-system. It has been shown in Ref. [l48| that if \ip) has a Schmidt rank of k, then at 
least k + 1 measurements are necessary. Also, if k is even (odd), an explicit decomposition into 2k (2k — 1) 
measurements has been found. For the case that \ip) is a projector onto a maximally entangled state in a 
d x ri -system and d is prime, it has been shown that W can indeed be decomposed into d + 1 measurements 
[3841 ] . hence for this special case the problem of the optimal decomposition is solved. 

A different interesting problem, which is somehow reverse to the original problem is the following: If we 
assume that some correlation measurements like o~ x ® a Xl a x ® a z ,a z ® o~ x and o~ x ® o~ x have been done, 
which states can then be detected with that data? This question is important in the context of quantum 
key distribution, and we will discuss it in Section 6.10. 
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Fig. 9. Layout of polarization entangled photo n sou rce. 
For details see the text. Figure taken from Ref. [385ll . 
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Fig. 10. Experimental results of the entanglement wit- 
ness for Werner states: Theoretical curve (solid), experi- 
mental results (dots), and experimental best fit (dotted). 
The horizontal dashed line indicates the transition be- 
twee n sep arable and entangled states. Figure taken from 
Ref. [3831. 



6.2. Implementation of a two-qubit witness 



In this Section we describe the experi men t by M. Barbieri et al. 385], which was one of the first im- 
plementations of entanglement witnesses ! 32 1 In the experiment they created two-qubit Werner states (see 
Sec. 3.4.5) of photons using parametric down-conversion and linear optics (see Sec. 5.5). 

The experimental setup can be seen in Fig. 9. The incoming photon with wavelength Ap and y-polarization 
arrives from the top from the laser. Then, after a mirror, it propagates to the — k p direction towards the 
BBO crystal. Here with some probability it is converted down (type I) into two photons with A := Ap/2 and 
H polarization. The two photons must have total momentum — k p due to momentum conservation. Thus, 
if one of them is propagating in a direction slightly up, the other has to propagate into a direction slightly 
down. Both of these photons are reflected on the mirror M, go back through the crystal and a lens. Another 
beam of photon pairs arises from the photons that went through the crystal without down-conversion, got 
reflected on the mirror and then converted down on the BBO crystal when propagating from the left to 
right. These photons also arrive to the same lens. Due to the wave plate between the mirror M and the BBO 
crystal, the photon pair from the first down-conversion arrived in \HH) state to the lens, while the other 
arrives in the \VV) state. As the contributions are indistinguishable, finally the setup generates photons in 
the state 



HH) 



i<t> 



\VV)). 



(126) 



There is a A/2 wave plate inserted before detector A, which leads to a l^smgiot) = ^Tjd-^^O — I'^O) state. 
It is also possible to "spoil" the entangled state Eq. (126) and create a mixture 



g(p) = p|*singlot)(* S inglot| + (1 



">4' 



(127) 



This can be done using the coated glass plate G and the wave plate before the detector A. This partially spoils 
the indistinguishability and leads to the Werner state. For large values of p the state created is entangled, 
while for small values it is separable. The entanglement of the state is detected by the following witness 
operator [see Eq. (123)] 



W 



I ^singlet) (^singlet | 



(128) 



32 For a similar experiment see Ref. [38 
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that can efficiently be measured locally. The experimental results are shown in Fig. f 0. Note that the witness 
detects the states also for p < 0.65, where no Bell inequality can detect the entanglement of that state. 

6.3. Witnesses for multi-qubit systems 

Let us now discuss witnesses for multipartite systems. Here, the general strategy is quite similar to the 
two-qubit case: The experimenter aims to prepare a state that is affected by some noise. The problem is to 
find a witness for this scenario that is robust against noise, and can be measured with a moderate effort. 
In this section we will mainly discuss the "standard" witnesses for genuine multipartite entanglement as 
described in Eq. (70). Different witness constructions will be discussed in the following sections. 

6.3.1. Three- qubit GHZ state 

Let us discuss first the three-qubit GHZ state \GHZ 3 ) = (|000) + |IlI))/%/2 as a simple example. The 
experimentally created state may be of the form 

Q = (l- P noisc)\GHZ 3 )(GHZ 3 \ + Pnoisc l (129) 

o 

As we already know from Eq. (73) the standard witness would read W — 1/2 — \GHZz){GHZ 3 \ and this 
detects the state for p no i S c < 4/7 ~ 0.57. 

For this witness, a decomposition into local measurements can be found as follows: first one decomposes 
VV into tensor products of Pauli matrice s, th en, one optimizes this decomposition in order to minimize the 



measurement settings. So one can write [387 1 



W = - (3 ■ 1 <g> 1 <g> 1 — 1 <£> a z ® a z — a z <g> 1 <g> a z — a z (g a z (g 1 — a x <£> a x (g a x + a x <£> a y <£> a y + 



x + o- y <g a x 63 o-y). 

l-l(g) o z <g a z - a z (g) 1 (g) a z -a z ®u z <g 1 - 2 • a? 3 



This decomposition requires only the measurement of the four correlations cf 3 , erf 3 and (a x ± a y /V2)® 3 . 
These are simple measurements on the Bloch sphere that can be implemented experimentally. From the 
data, also the fidelity of the GHZ state can be determined, as in Eq. (124). It can be further shown that this 
decompos ition is optimal, i.e., it is not possible, to determine (W) with three local measurement settings 



387, 38 



6.3.2. General states 

Let us discuss now general states. Using the presented methods one can calculate the witness and determine 
its robustness to noise straightforwardly for an arbitrary state \ip). To obtain a good local decomposition 
requ ires often some effort, especially proving that a given decomposition is optimal, is often very difficult 
.388]. However, for many interesting states results on this problem are known, and we summarize them in 
Table 1. Interestingly, one can also make some general statement about the re quired m easurement settings: 
For any pure state there exists a witness that requires 2N — 1 measurements 389L 390j , but the robustness 



to noise may be small. Furthermore, there exist observables, for which the local decomposition requires 
2 • 3 N ^ 1 /(N + 1) local measurements, which means that a local measurement of these observables re quire s 
nearly the same effort as state tomography, but specific examples of such observables are not known [3881 ] . 

6.4. Comparing witnesses with Bell inequalities 

Let us now compare the entanglement detection via Bell inequalities and entanglement witnesses in some 
more detail. This will also shed light on the underlying assumptions of different entanglement verification 
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Table 1 



Results on local decompositions of different entanglement witnesses for different states. 



# of qubits 


state 


witness 


maximal p no ise 


local measurements 


references 


remarks 


3 


\GHZ 3 ) 


±1- \GHZ 3 )(GHZ 3 \ 


4/7 


4 (optimal) 


[3871 


a 


3 


\W 3 ) 


|l-|W 3 >(W 3 | 


8/21 


5 (optimal) 


[3871 


b 


4 


\CU) 


\l-\CU)(CU\ 


8/15 


9 (optimal) 


[388, 3911 


a 


4 


1*2} 


fl-|*a)(*a| 


4/15 


15 


[145, 3901 


c 


4 


|£>2,4> 


§1 - |D 2 , 4 >(D 2 ,4| 


16/45 


21 


[213] 


d 


N 


\GHZ N ) 


§1 - \GHZ N )(GHZ N \ 


1/2. [1/(1-1/2")] 


N + 1 


[390] 


a 


N 


\W N ) 


^ll-\W N )(W N \ 


X/N- [1/(1 -1/2")] 


2N - 1 


[8, 390] 


b 


N 


\Gn) 


±l-\G N )(G N \ 


1/2- [1/(1 -1/2")] 


depends on the graph 


[261] 


a 


N 






1/2 • (N - 2)/[(N - 1)(1 - 1/2")] 


not known 


[213] 





a Witnesses that tolerate less noise but require less settings exist. See Section 6.6.1. 

b Witnesses that tolerate more noise with the same measurements exist . See Sections 6.8.2 and 8.2. 

c Witnesses that tolerate more noise and require less settings exist [l90l l390j | . 

d For witnesses that tolerate less noise with less settings see Section 8.2, for witnesses which tolerate more noise see Ref. Il9fj|| . 



methods (for a detailed treatment on this see Ref. [2lT]). 

When comparing different entanglement verification schemes, two main questions are of interest. First, 
one may ask for the power of the method in an experiment, i.e., to which extent it gives useful and relevant 
information about the state. Here, also the question for the experimental effort becomes relevant. A second 
main question concerns the underlying assumption for the scheme, and to which extent these assumptions 
are realistic in a given experimental situation. 

6.4. f. Usefulness 

Concerning the first question, we have seen in Section 5.2.3 that a violation of most of the Bell inequalities 
for multipartite systems does only exclude full separability of the state. But, as seen in the same Section, 
there are some Bell inequalities that can be used for proving genuine multipartite entanglement. However, 
up to now such constructions work only for few states (like the GHZ states) and there is no construction 
known for other important states, such as cluster states. Furthermore, the existing Bell inequalities for the 
detection of genuine multipartite non-locality require a high fidelity of the target state and an exponentially 
growing number of measurement settings. In contrast to that, we will see in Section 6.6.1 that one can 
write down witnesses for the GHZ state that tolerate more noise and need only two measurement settings, 
independently from the number of qubits. 

6.4.2. Underlying assumptions 

A further interesting question concerns the underlying assumptions of the different entanglement detection 
methods. As the loopholes are in most experiments not closed, Bell inequalities require the assumption that 
the parties act like as if they were space-like separated and that the detector efficiencies do not play a role. 
Apart from that, however, no more assumptions are needed. Especially, the measurement directions do not 
have to be fixed, a misalignment of the measurement directions does not invalidate the conclusion that the 
state was entangled. 

For entangle ment witnesses, more assumptions are needed. First, loopholes can also arise for entangle- 
ment witnesses |392| . Then, the validity of quantum mechanics is assumed. Most importantly, entanglement 
witnesses require that the experimentalist knows what he is doing: The measurements have to be m ade in 
the proper direction and the measurement apparatuses must work as they are intended to work 1 33 



For some schemes to deal with non-ideal detectors see Refs. [3931 . l394ll . 
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A further question is whether the use o f entanglement witnesses requires an assumption on the dimension 
of the underlying Hilbert space [395l l396| . In many experiments one works with multi-level systems like ions 
in a trap, however, for the theoretical description one considers it as a two level system and disregards the 
levels that are not occupied. The question arises, whether this might lead to false detection of entanglement. 
To discuss this in a simple setting, consider two qubits and the witness from Eq. (122). A measurement of this 
witness consists of the determination of the outcome probabilities for the eigenvectors \x ){x |, \y ){y |, 
and |2 ± )(z ± | for Alice and Bob. These probabilities are defined on an arbitrary dimensional Hilbert space, 
and if these probabilities violate the witness inequality, the corresponding reduced two-qubit state must be 
entangled. This implies entanglement of the total state l 34 1 

Finally, it should be noted that due to the finite number of copies available in any experiment, Bell 
inequalities and entanglement witnesses can only detect entanglement with a finite stat istic al significance. 
This is quantified by the standard deviation of the experimentally measured mean value 1 35 



6.5. Implementation of a multi-qubit witness 

Now we discuss an experiment by M. Bourennane et al. [l45j ] . which was the first experiment, where a 
multi-qubit witness has been implemented. In this experiment, the four-qubit singlet state [Eq. (61)] 
has been generat ed u sing the second order emission of parametric down-conversion and postselection as 
proposed in Ref. [2371 ] . 

One possibility to generate entangle d sta tes of polarized photons in a laboratory uses type-II spontaneous 
parametric down-conversion (SPDC) [38ll |. In this process, an ultraviolet (UV) laser pumps a nonlinear 
/3-Barium-Boratc (BBO) crystal (see Fig. 11). Multiple emission events during one pump pulse lead to the 
emission of the state 

\ip) ~ exp[-ia(aJ )ff 6j iV + «4 )V &J iH )] 1°) ( 131 ) 

distributed onto two modes, ao and bo- Here, aj H denotes a creation operator, creating a horizontally 
polarized photon in mode ao; b v creates a vertically polarized photon in mode bo and |0) denotes the 
vacuum. In the second order of this process only the quartic terms are relevant, leading to ~ (a H 6 J v + 
a' Q v bl H ) 2 \0). This results in the state 

IV>> ~ mj2V) bo + \2V) ao \2H) bo +2\lH,lV)JlH,lV) bo , (132) 

where \2H) denotes a state of two horizontally polarized photons in mode ao, and \1H, lV) bo denotes one 
horizontally and one vertically polarized photon in mode bo, etc. 

Two beam-splitters are used to distribute the photons in the four modes a,b,c and d. Then postselection 
is used and one can directly calculate that fourfold coincidences in all the four modes a,b,c, and d can only 
occur when the state 

W ~ \H) a \H) b \V) c \V) d + \V) a \V) b \H) c \H) d + \{\V) a \H) b + \H) a \V) b ) (\V) c \H) d + \H) c \V) d ) (133) 

is produced. This state is, up to a sign which can be removed by a local change of the basis the four-qubit 
singlet state | \I/ 2) ■ In the experiment, this state was created with a count rate of around five four-photon 
events per minute. 

The witness for this state, is given by Eq. (76) as W* 2 = |l — I^K^I- This witness is equivalent to 
a measurement of the fidelity and requires fifteen measurement settings. It should be noted, that in the 
meant ime witn esses have been found that require significantly less measurements and tolerate even more 
noise 190l 39C| . The measurement results are given in Fig. 12. The fact that the coincidence probability 



34 To be precise, if gtot is the total state on a d X d system, violation of the inequality implies entanglement of the (not 
normalized) state g r cd = V^? <g> Vg ^ Qtot 

' P A ) ® with ^ t2) = 1 0) <0| + This implies entanglement of Qxox, as the 

projection is local. 

35 For an interesting proposal to determine the probability of entanglement via Bayesian updating see Sec. 6.9.2 
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Fig. 11. Setup for the generation of the four-photon singlet Fig. 12. Fourfold coincidence probabilities for the fifteen 

state 1*2 >• See the text for details. The figure is taken from measurement settings, needed for the fidelity m easu rement 

Ref. |242l . °f the |*2)-state. The figure is taken from Ref. |l45|l . 

results are similar for different settings is a signature of the U <S> U (£> U <g> {/-symmetry of the singlet state. 
From these data, one obtains 

(W* 2 ) = -0.151 ± 0.01 and F<p 2 = 0.901 ± 0.01 (134) 

which clearly proves the four-partite enta nglement. The state l ^) can then be further used to encode a 



qubit in a decoherence f ree s ubspace [244j or for secret sharing [243J. With a similar setup, three-qubit W 
states can be generated [21fj| . 

6.6. Stabilizer witnesses 

In this Section we will explain a method for constructing witnesses that are very convenient for experi- 
ments, since they can be implemented with very few, often only two measurements. 

6.6.1. The main idea and GHZ states 

The main idea of the stabilizer witnesses is the following: If an existing witness W requires many measure- 
ment settings, one can look for witnesses that might be less powerful, but easier to measure. The relationship 
between the old witness that is difficult to measure and the new one is demonstrated in Fig. 13. 

How can such a new witness be constructed? Let us first consider a general observable W. If one can find 
an a > such that 

W-aW>0, (135) 



then Tr(gW) < implies Tr(gW) < 0[££] Therefore any state W detects is also detected by W. Moreover, 
if we know that W detects states with genuine multipartite entanglement, so does W. The question remains, 
how W can be found. It turns out that stabilizer theory is a well suited tool for this task, but it should be 
noted that it works also for non-stabilizer states. 

Let us shortly recall some facts about it, see also Section 3.4.3. Stabilizing operators are locally measurable 
operators that have the state |^) as an eigenstate with eigenvalue +1, 

|*> =£*!*}. (136) 



36 The condition in Eq. (135) is indeed equivalent to the "finer" relation in the sense of Section 2.5.2, i.e., VV is finer as W iff a 
relation as in Eq. (135) holds. 



50 



original witness — ► 



<W " — witness with less 

slab 

measurements 




Fig. 13. Schematic view of the construction of stabilizer wit- 
nesses: Instead of measuring the original witness W, (which 
might be difficult) one considers the witness W ata ,h- This 
witness detects fewer states, however, it is often much eas- 
ier to measure than the original witness. The witness W sta b 
can be constructed from the stabilizing observables for the 
target state. 
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Fig. 14. (a) Measurement settings needed for detecting gen- 
uine multi-qubit entanglement close to GHZ states with the 
Mermin inequality [sec Eqs. (107, 141)]. For each qubit the 
measured spin component is indicated, (b) Settings needed 
for the entanglement witnesses detecting entangled states 
close to GHZ states [Eq. (139)] and (c) cluster states [Eq. 
(143)]. 



It is common to express Eq. (136) saying that operator Sk stabilizes the state |\E'). If operators SV- and 
Si stabilize then so does their product SkSi. It is easy to s ee that these stabilizing operators form a 



commutative group. This group was named stabilizer in Ref. 221). If |^) is an N— qubit quantum state then 
the group of stabilizing operators has 2 N elements. The group has N generators gk from which all elements 



are possible products of the Sk- Finally, the projector to can be written as 

N 



1 2 N 1 _L 



fc=i fe=i 

Let us show a simple example and consider the witness operator Eq. (73) detecting genuine multipartite 
entanglement around a GHZ state. The elements of the stabilizer for the GHZ state are products of the 
generating operators 

N 

g (GHZ N ) . = -Q Xfe; g( GHZ N ) . = Zk iZk for k = 2 ^ ^ ^ (138) 
fe=l 



Then, one considers the following observable 
VV GHZn := 31 - 2 



(GHZ N ) N (GHZ N ) 

9l ±J_ , TT Ek ±J. 

O ^11 O 



k=2 



(139) 



Knowing that Wghz n = 1/2 • 1 — \GHZn)(GHZn\ is a valid witness, one can verify Eq. (135) for a = 2 : 
The left hand side of Eq. (135) is then diagonal in the GHZ-basis, and it suffices to prove that the diagonal 



elements are positive, which can be directly checked [26413971 ]. Therefore, Wqhz n is a valid witness for the 



detection__of genuine multipartite entanglement around the GHZ state. For the special case of three qubits, 
it reads W GHZ3 = §1 - X X X 2 X 3 - \{Z X Z 2 1 3 + Z X 1 2 Z 5 + \ X Z 2 Z Z ). 

The first remarkable thing is that WchZn requires only two measurement settings: The first term in the 
square bracket can be measured with a {X, X, X, X} setting, the second with a {Z, Z, Z, Z} setting. 
This is shown in Fig. 14(b). Second, the witness Eq. (139) is quite robust against noise. It detects the GHZ 
state as entangled if the fraction of the white noise fulfills 

Pnoisc < 3 _22-iV ( 14 °) 
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Note that for large N the right hand side of Eq. (140) converges to 1/3 J 37 1 Therefore, GHZ states with 
a fixed amount of noise can always be detected with two measurement settings, independently from the 
number of qubits. This shows that entanglement detection does not become necessarily more complicated 
when increasing the number of qubits. 

If more than two settings can be measured then criteria with a higher noise tolerance can be obtained. 
At this point it is interesting to show a connection between stabilizer witnesses and the Mermin inequality 
Eq. (107). Let us consider the following witness 

WMcrmin := 2 N ~ 2 1 - £ X X X 2 X^X A ■ ■ ■ X N + ^ ^2X3X4 • • • X N - ^ Y 1 Y 2 Y i Y i X 5 ■ ■ • .V v • 

7T 7T 7T 

(141) 

where ^2 indicates that every term represents the sum of all of its permutations. This witness is just 
constructed from the Bell operator of the Mermin inequality given in Eq. (107). We used the bound for 
genuine multipartite entanglement for this inequality presented in Refs. 353l 1397 1. Thus, WMermin detects 
genuine multipartite entanglement. Note that the witness is also a stabilizer witness since all the terms in 
Eq. (141) are stabilizing operators of the GHZ state. This witness detects the noisy GHZ state as entangled 

if Pnoisc < 1/2. 

6.6.2. Linear cluster states and other graph states 

Similar methods can be applied to create witnesses with few measurements for cluster states and in 
general, graph states, introduced in Sec. 3.4.3. The generators for the group of stabilizing operators of the 
cluster state are 

gf N) := X X Z 2) gf N) := Z k _ x X k Z k+l for k = 2, 3, JV - 1, 9 ^ n) := Z N _ X X N . (142) 

For detecting genuine multipartite entanglement we have already presented a witness Eq. (139). Based on 
arguments similar to creating witnesses for the GHZ state, one can construct a witness that also detects 
genuine multipartite entanglement, however, it is easier to measure 



n L^+n^ 1 



■ even k odd k 

It detects the cluster state as entangled if 



(143) 



f(4-4/2-)- 1 for even N, 

Pnoise < <^ N + 1 JV _ 1 (144) 

I [4-2(1/2 a +1/2— )]-* for odd N. 

For large N the limit is p n oisc = \- The measurement settings needed are shown in Fig. 14(c). 

The ideas presented for GHZ and cluster states can be straightforwardly generalized for graph states. 
This requires the notion of the colorability of a graph. A graph is called fc-colorable, if there exist fc disjoint 
subsets Mi, Mfe of the set of vertices, such that for any j there are no edges between any pair of vertices 
in Mj. For instance, the star graph (see Fig. 4 in Section 3.4.3) and the cluster graph are two-colorable, 
while a three- vertex ring graph is three-colorable. 

Using this notation, for a fc-colorable graph state |Gjv) the following witness detects genuine multi-partite 
entanglement. 



W Gjv := 31 



±_ , (G N ) _ 

" , ( 145 ) 

j=i ' ieMj 



This witness requires the measurement of k settings, since each of the terms in the sum can be measured 
with a single setting, namely with a setting measuring {Xi}i^Mj an d {Zi}igM r The robustness to noise 



7 It can be shown, that the witness in Eq. (139) is optimal within a certain class of witnesses with two measurements 1397 
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depends on k. It should be noted, however, that the notion of £;-colorability is a property of the graph 
(and not the state) and if a graph state can be represe nted (up to local unitaries) by different graphs, the 
colorability properties of these graphs may be different [22d| . 

Such witnesses have been used in m any experimen ts. That i s, the two-setting witness for cluster states 



Eq. (143) was used first in Ref. [2321 ] . later in Refs. 234L 235 1 . In another experiment, six photon graph 



states were created in a photonic experiment and genuine multipartite entanglement was detected with a 
witness using six measurement settings Q . The advantage of using six, rather than only two settings is the 
larger robustness to noise of the witness. We will discuss in Section 6.7 an experiment for the generation of 
the four-qubit cluster state in detail. 



6.6.3. Non-stabilizer states 

Similar ideas can also be used for the detection close to states that are not stabilizer states. For such states 
it is not possible to find 2 N operators such that they stabilize the state and also they are the tensor-products 
of single-qubit operators. An example for such a state is the three-qubit W-state IW3). If nonlocal operators 
are also considered, however, the following operators, which are the sum of several operator products, 
stabilize this state: 



g (w 3 ) .-l{Zx + 2Y 1 Y 2 Z 3 + 2X 1 Z 2 X 3 ), 
g (w 3 ) . = 1 ^ + 2 Yl Z 2 Y 3 + 2Z X X 2 X 3 ) . 



9^ 



= \{Z 2 + 2Z 1 Y 2 Y 3 + 2X 1 X 2 Z 3 ), 



(146) 



There is only a single state, the W-state, that gives +1 for all the three operators, as in the case of stabi- 
lizer states. These operators commute with each other. Multiplying these operators gives further operators 
that stabilize the W state thus these operators form an eight-element group. The g^ 3 ^ operators are the 
generators of the group. However, now the choice of th ese eight stabilizing operators is not unique. A simple 
method for finding such opera tors is given in Ref. [397J. Using these operators, and the strategy from above, 
one can find now the witness 13971 



(147) 



% = jl + 2Z X Z 2 Z 3 - \ ^(2A fc A z + 2Y k Y l - Z^), 



which requires only three measurement settings. For another application of similar ideas see Ref. [395 



6.6.4. Fidelity estimation 

Detecting genuine multipartite entanglement is closely related to measuring how close the experimental 
state is to a given multipartite quantum state. This can be seen noting that the projector-based entanglement 
witnesses Wqhz n and Wc N of Eq. (73) were simply based on measuring the projector, that is, the fidelity 
F = (V-'l^cxpIV') with respect to a highly entangled quantum state Thus, methods similar to the ones 
for detecting genuine multipa rtite enta nglement with few measurements also work for estimating the fidelity 
with few measurements 391 . 397 . 3991 ] . In particular, one can consider the relation 



\GHZ N )(GHZ N \ > 



(GHZ N ) , N (GHZ N ) 

91 9 +1 +n gfc o 



fe=2 



1. 



(148) 



which is essentially equivalent to Eq. (135) for the stabilizer witness for the GHZ state. This relation implies 
that measurement of the right hand side of Eq. (148) can be used to give a lower bound on the fidelity 
of the GHZ state. Equivalently, the stabilizer witness allows to give a lower bound on the fidelity via 
F G hz n > (1 - (W GH z N ))/2. 

There are other ways to estimate the fidelity, similarly to Eq. (148) [390|,[39l|. To see the main idea, let us 
consider the example of a four qubit cluster state \CL 4 ) = (|0000) + |1100) + |0011) - |llll))/2, for which 
the usual fidelity-based witness is given by W = 1/2 — \CL4)(CL,4\. Then, the observable 
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VV = 5I- |CL 4 )(CL 4 | + |OL 4 )(Cl/ 4 | 

= i[i - |oo)(oo| ® (|oo)(ii| + |ii)(oo|) + (|oo)(ii| + |ii)(oo|) ® |ii)(ii|)] (149) 

with |CZ4) = (— 10000) + 1 1100) + |0011) + |llll))/2 is, due to Eq. (135), also an entanglement witness and 
allows to estimate the fidelity via F > 1/2 — (VV). As the off-diagonal ter ms o f VV in the computational 



basis are of a simple form, VV can be measured with four measurements only [39 1| , whi ch is significantly less 
than the nine measurements for W. Simi lar ideas can be applied to other states |390j and have been used 
in several six-photon experiments 0, |236| . 

6.6.5. Witnesses for full separability 

Up to now, we considered stabilizer witnesses that detect only genuine multi-qubit entanglement. Typically, 
the noise in an experiment must be relatively low if we want to detect this form of entanglement. On the 
other hand, if we want to detect entanglement in general, including partial or biseparable entanglement, 
then typically in an experiment much higher noise is allowed. 

The basic idea behind detecting entanglement with stabilizing operators is that a product state cannot 



give +1 for two stabilizing operators that do not commute locally [397[. Based on this idea, and after simple 
algebra, a witness detecting entanglement around an iV-qubit GHZ state can be obtained as 

W^ HZ ^ := 1 - g[ GHZN) - g£ HZ »\ (150) 

where m — 2, 3, N. The proof is based on the Cauchy-Schwarz inequality. Using this and (Xi) 2 + (Zi) 2 < 1, 
we obtain for pure product states 

{g (GHZ N)) + {g (GHZ N)) = {Xl){X2 )...(X N ) + ( Z m - X ){Z m ) 

< \{X m - X )\ ■ \{X m )\ + |(Z m _i)| ■ \(Z m ) \ < v/(A^-i) 2 + (Z m _i)V(*m) 2 + (Z m ) 2 < 1. (151) 

By convexity, the bound is also valid for mixed separable states. This proof can straightforwardly be gener- 
alized for arbitrary two locally non-commuting elements of the stabilizer of any graph state. Witnesses can 
be constructed with more than two elements of the stabilizer as 

VV|f HZn) := 1 - g[ GHZN) ~ g^ HZN) - g[ GHZN) gi GHZN) (152) 

for m = 2,3, ...,N that rule out full separability. The proof is similar to the proof of Eq. (150). These 
witnesses tolerate noise if j? n oise < \ and |, respectively. As can be seen, more terms make it possible 
to have higher noise tolerance. One can even construct a witness from the Bell operator of the Mermin 
inequality 

^Mermin ■= ^ 1 ~ M 2 X S X A • ■ • X N + ^ ^2X3X4 • ■ • X N - ^ 1^2*3*4 • • • X N + ... - .... 

TV TV TT 

(153) 

This differs from Eq. (141) in having a smaller constant term corresponding to the relaxed conditions for 
biseparable entanglement comp ared to genuin e multipartite entanglement. The constant comes from the 
calculations presented in Refs. (l90l . I2651 [267} . While the number of terms increases with the number of 



qubits for W^[ crmin , the noise tolerance is also increasing. 

For the cluster state similar witnesses are [393, EoqI W {Gn) := 1 - 5 [ Cjv) - g^p for k = 1, 2, N - 1, 
and, with a better noise tolerance W {Gn) := 1 - gf N) - g {G + f - gf^g^ for k = 1,2, ...,N - 1. Both 
witnesses involve only the qubits of a quadruplet and tolerate noise if p n oise < \ and |, respectively. It i s 



also possible to obtain a series of witnesses for which the noise tolerance approaches 1 for increasing N 397 1 . 



6.6.6. Mixed stabilizer states 

So far, entanglement conditions were presented that detect entanglement in the vicinity of pure multipar- 
tite states. It is also possible to use witn esses made out of stabilizing operators for detecting entanglement 



in the vicinity of mixed stabilizer states 397]. The general idea is considering witnesses of the from 
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M 

Wmixed := const. - ^ 9k, 
fc=l 



(154) 



where gk are some of the generators of the stabilizer for some graph state. If M = N then there is only a 
single quantum state giving the minimum for (W m ixcd)- However, if M < N then the minimum is degenerate 
and the states giving the minimum make up a subspace of dimension (N — M +1). Then, if one can determine 
the constant in Eq. (154) such that (yV m i XC d) is positive on product states, yet it is still negative for some 
state, then W m i XC d is a valid witness. 

As an example, let us consider the biseparable (but not fully separable) mixed quantum state 

qs ■= + irxri), r> := ^(ioo> + |n» ® io>, m := ^=aoo) - in» ® ix>, (155) 

which occurs if the fourth qubit of a cluster state is traced out. The states l^) are stabilizer states. Any 
mixture of states is stabilized by the common elements of > . The generators of the stabilizer of these 
mixed states is given by S^ 3 ) = {Z1Z2, X1X2Z3}. Based on that, the following witness can be constructed 



W ea := 1 - Z X Z 2 - X X X % X Z . (156) 

The constant term in Eq. (156) can be obtained via methods similar to the ones used for obtaining a bound 
for product states for the sum of two stabilizing operators used for Eqs. (150). 



6.7. Experimental generation of cluster states 

Now we describe an experiment from K. C hen et al. [2351 ]. where a four-qubit cluster state has been 



observed and a stabilizer witnes s has been used l 38 ! This experiment used hyper-entanglement of photons for 
the creation of the cluster state [40 lj . This means that several degrees of freedom of each photon are used to 
carry the entanglement. In the given experiment, each photon represented two qubits, one polarization qubit 
and one spatial qubit. This technique is becoming m ore and more i mportant in multi-photon experiments 



as it allows to ex tend the Hilbert space significantly [402, 1403 . 1404] and also complete Bell state analysis 



becomes possible (405]. Quite recently, this technique has been used to realize a ten-qubit GHZ state with 
five photons [Io| . 

The setup of the experiment is shown in Fig. 15. A UV laser pulse passes two times through two combined 
BBO crystals. This leads either to emission of \HaHb) + |VaVb) in the upward direction (modes La, Lb) 
or to emission of \HaHb) — \VaVb) in the downwards direction, i.e., modes Ra, Rb- Here, H and V denote 
horizontal and vertical polarization. If the photons overlap perfectly at the beamsplitters BS\ and BS2, the 
state will be a superposition 

\1>) = ^[(\H A H B ) + \V A V B ))\L A L B ) +e*+{\H A H B ) - \V A V B ))\R A R B )]. (157) 

By adjusting the concave mirror, = can be obtained. Then the state is equivalent to the cluster state 
= (|0000) + 10011) + |1100) - |11H», if the replacements \H/V) h-> |0/1) and \L/R) ^ |0/1) are done. 
Since only two photons are required, the count rate in the present experiment is around 1.2 • 10 4 per second, 
which is about four orders o f magnitude higher than in experiments that produced the same cluster state 



using four photons 232, 12331 ] . 

In order to measure the state, the elements as in Fig. 15(b) have been used. A measurement without 
any beamsplitter as in (i) first reads out whether the photon is in L/R and also the polarization of the 
photon. Hence it can be used to measure Z sva ±X vo \ at one photon. Using a beamsplitter as in (ii) leads to a 
Hadamard transformation on the spatial modes, such that A spat Z po i can be measured. The problem is that 
such a setup is sensitive to small path length changes, hence it requires a significant amount of adjustment 



3 At the same time, a similar experiment has been performed by Vallone et al. 234] . 
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Fig. 15. (a) Setup for the generation of four-qubit cluster states using two photons and hyper-entanglement, (b) Required setups 
for the measurement of the polarization qubit and the spatial qubit of one photon. See the text for further details. Figure is 
taken from Ref. Ijjl . 



and is only stable for relatively short times ! 39 1 Therefore, despite of the high counting rate, an analysis of 
the state should require only few measurements. 
One can use the witness 

W = 21 - ipfiXa 3 ^4 + XxXxZzU + lll 2 ^3^4 + I1Z2X3X4 + ^ll 2 ^3^4 + Z 1 Z 2 1 3 U) (158) 

which is just the witness from Eq. (143) adapted to the form of the cluster state used here. This witness 
can be measured with the settings X1X2Z3Z4 and Z1Z2X3X4. This corresponds to the setups in Fig. 15b 
(i) and (ii). The experimental data give 

1 - (W) 

(W) = -0.766 ± 0.004 and F > ^— '- = 0.883 ± 0.002, (159) 

clearly confirming the four-qubit entanglement. Furthermore, this cluster state has then been used to demon- 
strate two-qubit gates. 

The nature of multipartite entanglement in experiments with hyper- entangled states deserves some more 
comments. First, it should be stressed that due to the hyper- entanglement not all qubits are spatially 
separated, hence the entanglement refers to different degrees of freedom, and not different particles. Second, 
in multiphoton experiments as in Ref. [l(| the process of postselection has a different effect now. Ideally, 
in any experiment, one aims to generate some entangled state g which can then be distributed and used 
it as a resource for some task to be determined later. In usual postselection experiments, a state gi n i is 
prepared, which contains different numbers of photons. Postselection results in a "state" g ps , which, however, 
is destroyed and cannot be stored. However, as usual postselection is a local projection, entanglement of 
Q ps implies entan glem ent of which may therefore be considered as the (in principle) distributable and 



storable resource [2lJ . When hyper- entanglement is used, the postselection acts collectively on two qubits at 
the same time, and there is no direct connection between the multi-qubit entanglement of g ps and the one 
of gini. For many applications, however, this does not matter. 

Finally, it is interesting that if one wants to detect entanglement only b etwe en the different degrees of 
freedom of hyper- entangled particles, stabilizer witnesses can again be used |406l | . 



Only recently, using a Sagnac-like interferometer, stable interferometers for this problem have been introduced jl0| . 
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Fig. 16. (a) Level scheme of the 40 Ca+ ion. The D 5 / 2 and 

the Sx/2 levels are used as the basis states of the qubit. Fig- 17- Absolute values of the reconstructed density matrix 

(b) Schematic view of the motional levels describing the bus for the four-qubit W state, 

mode, (c) Combined level scheme of the motion al qu bit and 
the internal qubit. The figure is taken from Ref. 140711 . 



6.8. Eight-qubit W states in ion traps 

In this Section we explain an experiment in which eight qubits have been entangled using trapped ions 
Q. Apart from being the most advanced experiment to generate multi-particle entanglement in ion traps, 
it has also some interesting aspects from a theoretical point of view. First, a witness different from the 
projector-based one was used, second, the usage of local filters enabled the detection of entanglement even 
in a very noisy environment. 



6.8.1. Description of the experiment 

Let us first shortly describe the setup. Detailed overviews on quantum information processing with ion 
traps and their physical principles can be fo und in Refs. 36|, 37, 38, 3t| and for a more detailed description 
of the present experiment see Refs. EJEqzI. 

In the experiment strings of up to eight 40 Ca + ions were trapped in a linear Paul trap. The D 5 / 2 (= |0)) 
and the S1/2 (= |1)) levels of the ion are used as the basis states of the internal qubit [see Fig. 16(a)]. The 
ions can be individually addressed via laser pulses on the |0) <-*• |1) transition. 

For the interaction between the qubits, the vibrational excitations of the ion string are used. These 
excitations arise from a harmonic oscillator potential [see Fig. 16 (b)], and the two lowest levels |0) m and |l) m 
are referred to as the motional qubit. Via driving the |0) m |l) «-> |l) m |0) transition ( see Fig. 16 (c)), each ion 

. rm . ri"ii . 



can individually interact with the motional qubit, allowing finally two-qubit gates 408, 
The idea to prepare the iV-qubit W states 



\W N ) = — =(|00...01) + |00...10) + ... + |10...00)) (160) 
v N 

can be sketched as follows. First, the state |0) |10...0) is created, where all but the first internal qubit are 
in the |0) state. Then, the motional qubit is partly excited and distributed over all the ions. That is, the 
state is changed according to the sequence 
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Table 2 

Fidelities and mean values for the witnesses of the 7V-qubit W states for different values of N. 



N 


4 


5 


6 


7 


8 


F 


0.846 ±0.011 


0.759 ± 0.007 


0.788 ± 0.005 


0.763 ± 0.003 


0.722 ±0.001 


<W> 


-0.460 ± 0.031 


-0.202 ± 0.027 


-0.271 ± 0.031 


-0.071 ± 0.032 


-0.029 ± 0.008 



|0>„J10...0) -> -!=|0> J100...0) + ^_I|1)J000...0) 



-^=|0) m |100...0) + -L|0)J010...0) + — =?|l) f j000...0) 



7W ia >- 



1100. 



|010...0) 



1000. ..1) 



(161) 



Details on the required pulse sequences are given in Refs. [|, 407 1. This scheme allows to generate W states 



/2 



level is t = 1.16 s. It should be noted 



in a time of 500 - 1000 (is, for comparison, the life time of the Dn 
that W states in ion traps can also be prepared in different ways 413 [ 

In order to characterize the experimentally generated state g cxp , state tomography has been performed. 
For that, all possible correlations of the type (Tj®^ ®...®crfe have been measured. Repeating each correlation 
measurement at least a hundred times requires 100 • 3 N repetitions of the state generation, for N = 8 qubits 
this amounts to a total measurement time of 10 hours. 

Due to the statistical nature of the measurement, the raw data from the correlation measurement s do 
not result directly in a density matrix that is positive semidefinite. Therefore, usin g the results of Ref. 414 1 
a maximum- likelihood approximation to the data has been performed. In Ref. [4141 ] . it was shown that 
the maximum- likelihood estimation is a fixed point of an iterative map, where the iteration depends on 
the measurement data. This fixed point was numerically determined by iterating the map, starting from a 
maximally mixed state until convergence was reached. Due to the exponential growing size of the density 
matrix this iteration requires a considerable effort, for the eight qubit case it requires more than one day. As 
a result, valid density matrices have been obtained, for the case N = 4 the absolute values of it are plotted 
in Fig. 
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From these matrices, the fidelity of the W states can be calculated, see Table 2. In addition, 
given the estimated density matrix one can generate with a Monte Carlo simulation further data sets, which 
were used to calculate the error bars. 



6.8.2. Analysis of the states with witnesses and local filters 

A first tool for the analysis of the states is the fidelity-based witness VV = ^r"l — |Wjv)(Wjv|- However, 
as can be seen from the fidelities in Table 2, this witness does not detect any entanglement for N > 5. To 
improve this witness, one starts with the following observation: The projector onto the W state is a matrix 
that is acting on the subspace with one excitation (= |1)) only. Therefore, if one calculates the overlap with 
a biseparable state \<p) = \a)\b) it suffices to assume that \a) and \b) have maximally one excitation, as any 
contribution with two ore more |1) in \a) or \b) has a vanishing overlap with the state |Wi\r)- From this 
consideration, one can directly see that 

N - 1 

W = —fi-V<2-\W N ){W N \, (162) 

where V<2 is the projector onto the space with maximally two excitations, is a valid witness P^H 

For larger N, the witness in Eq. (162) is already a significant improvement, as the difference between V<i 
and 1 increases exponentially. However, one can make the more general ansatz 

W = aP + (3Vi + iV2- \W N ){W N \, (163) 



The numerical values of the density matrices are given in the online material of Ref. [8j . 
For example, for iV = 3 we have "P< 2 = 1 - | 111) <111| . 
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where the Vi are now projectors onto the subspaces with exactly i excitations. It remains to determine the 
allowed values for a, fj and 7 in order to guarantee that W is indeed non-negative on all biseparable states 
|0) = |a)|b). From the symmetry of the W state, one can conclude that it suffices to prove non- negativity for 
I a) of the form \a) = a |00...00) + ai(|00...01) + ... + |10...00)) and similarly for |6) [|. Therefore, the resulting 
conditions on a, /3 and 7 are given by simple optimization problems of a polynomial in four variables with 
two normalization constraints, which can numerically simply be solved for any N. 

For given experimental data, one can choose th en t he tuples (a, f3, 7) that yield the most negative value. 
In practice, the further condition a — 7 was fixed ! 42 1 In order to make the values comparable, all witnesses 
have been normalized such that Tr[W(l/2 N )] — 1, i.e., the mean value for the maximally mixed state equals 
one. The witness values in Table 2 for 4 < N < 7 have been obtained in this way. 

For the case N = 8, however, this construction still does not yield a witness with a negative expectation 
value. In order to prove still the presence of entanglement, one can use the help of local filtering operations 
1Z3, Hill, Ell - Such a local filter is a map 



g h-> g = JV ■ {Ti <g> F 2 ® ... ® -FsM-T 7 ! ® T\ ® ■■■ ® T\) (164) 

where the Ti are arbitrary invertible 2x2 matrices and JV denotes the normalization (see also Section 2.5.3). 
As the Ti are invertible, the state g is entangled (or biseparable) if and only if g is entangled (or biseparable) 
and it suffices to detect the entanglement in the filtered state. Equivalently, one may also consider the class 
of witnesses 

W = <g> T\ ® ... ® Tl)yV{T x ®T 2 ® - ® F s ) (165) 

and try to find filters such that the mean value (W) gets negative. Here, one also has to renormalize the 
filtered witness with a factor JV [e.g., by fixing Tr(VV) = Tr(W)] in order to make the witnesses' values 
comparable. 

The optimization over all filters Ti can be performed by first optimizing T\. then JF 2 , etc., until a negative 
mean value can be found. With this technique, the mean value for the eight-qubit witness could be decreased 
significantly, and finally a negative value has been found. 

It should be noted that the filter technique is not restricted to the case where complete information about 
the state via state tomography is available. As shown in Refs. 1(| 39C| some witnesses that require only 



few measurements can be improved b y sp ecial filters that leave the required measurement settings invariant. 



Finally, it should be added that Ref. [4151 ] gave an interesting scheme to implement witnesses as in Eqs. (74, 
162) in linear optics experiments. 



6.9. Estimation of entanglement measures 



In the previous Sections we introduced methods for the detection of entanglement in experiments. More 
generally, one can ask how entanglement in experiments can be quantified. This l eads to the problem, how 
to quantify entanglement, if only incomplete information about the state is given [416j . 

In this Section, we will explain methods to quantify entanglement using entanglement witnesses. It should 
be noted that there are other proposals for the determination or estimation of entanglement measures in 
experiments, they are described in Sections 7.2 and 7.4. We will first describe a ge neral method t o derive a 
lower bound on a generic entanglement measure from the mean value of a witness 417 . 418l 419j . Then we 
will discuss results for specific witnesses and entanglement measures. 



6.9.1. A general method using the Legendre transform 

Let us consider the following situation: In an experiment, an entanglement witness W has been measured, 
and the mean value (W) = Tr(gW) = w has been found. The task is to derive from this single value a 



42 In Ref. Q the witness was found in a different way, leading to a different notation. However, it is equivalent to the construction 
in Eq. (163) with = 7. For the experimental data of Ref. a possible additional variation over 0,7 would not yield a big 
improvement. 
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Fig. 18. Geometrical interpretation of the Legendre transform. The dotted line is a general affine lower bound, and the dashed 
line is an optimized bound with the minimal /i . By varying the slope, one finds an affine lower bound, which is tight for a given 
wi (solid line). This figure is taken from Ref. [41911 . 

quantitative statement about the entanglement present in the quantum state. That is, one aims to derive 
statements like 

Tr(gW) = w => E(g) > f(w), (166) 

where E(g) denotes some convex entanglement measure or another convex figure of merit, which one is 
interested in (e.g. the fidelity). Clearly, one aims to derive an optimal bound f{w) and an estimate is 
optimal, if there is a state go with Tr(goW) — w and E(go) = f(w). 

In order to derive such lower bounds, one can use the Legendre transform |42dl . l42lj | of E for the witness 
W, defined via the maximization 

E(W) = sup{Tr(Wg) - E(g)}. (167) 

Q 

As this is defined as the maximum over all g, we have for any fixed g that E(W) > Tr(Wg) — E(g), hence 
E(g) > Tr(Wg) — E(W). The point of this rewriting is that the first term on the right hand side is the 
given measurement data, while the second term can be computed. Therefore, a measurable bound on E(g) 
has been obtained. 

More generally, one can use the fact that Tr(gW) = w is equivalent to knowing Tr(gXW) = Xw for any 
A. Therefore, optimizing over all A yields 

E{g) > swp{XTr(Wg) - E(XW)}. (168) 

A 

This is already the optimal bound in Eq. (166). This follows from the geometrical interpretation, of the 
Legendre transform (see Fig. 18): For a given convex function f(w) an affine lower bound g{w) — Xw — /i 
with a given fixed slope A is optimal, if ^ opt = sup tu [Au' — f{w)], which is exactly the Legendre transform. 
Such a lower bound is optimal for one wo. To obtain the optimal bound for any given w\ one optimizes over 
the slope A, corresponding to Eq. (168). Since / is convex, we obtain for each w the tight linear bound. 

Note that it is not needed that W is an entanglement witness, we may consider an arbitrary observable 
instead. Also one may consider a set of observables W = {Wi, ...,W rl } at the same time, by introducing a 
vector A = {Ai , A„} and replacing AW by J2k ^kWk- 

In any case, the main problem in this scheme lies in the calculation of the Legendre transform in Eq. (167). 
The difficulty of this task clearly depends on the witness W and on the measure E(g) chosen. For this 
problem, the following results have been obtained: 

(i) For entanglement measures defined via the convex roof extension [see Eq. (83)] the optimization can 
be simplified to 

E(W) = sup{(^|W|V) - E(\1>))}, (169) 



i.e., it suffices to optimize over pure states only [41 
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(ii) For the entanglement of formation [see Eq. ( 88)1 one can design a simple iterative algorithm that 
can perform the calculation numerically [417j |. For spec ial witnesses of the type W = \ip) {4>\ Tb the 
Legendre transform can be determined analytically [418j. C oncerning the concurrence, one can find a 
similar algorithm as for the entanglement of formation [419(. 

(iii) For the geometric measure of entanglement, one can also construct an iterative algorithm for general 
witnesses. F urthe r one can compute the Legendre transform for a witness like W = at — \^)(tp\ 
analytically 4171 ] . From this one can show that if an experimentally generated state has a fidelity F 
with the target state | ?/>), and the geometric measure of \tp) is E, then for F > 1 — E the geometric 



measure is bounded bj 



■a 



E G >^{l-F-E + 2EF + Ul^[(-in2F-ir-l]}. 



(170) 



For some other witnesses one can analytically estimate the Legendre transform 418l 419| |. 



6.9.2. Other relationships between witnesses and entanglement measures 

In order to investigate the relationship between witnesses and entanglement measures from a different 
viewpoint, one can ask whether the negative expectation value can directly serve for an entanglement 
quantification. Therefore, one considers the quantity [422j 

E(q) = max{0,max[-Tr(e>V)]} (171) 

where the optimizations is restricted over C, which is a certain compact subset of all witnesses. This E(g) is 
by definition already convex. 

Then, by choosing C appropriately, one can obtain different E(g) that can be used for the entanglement 
quantification: If C is the set of witnesses with W < 1 then E(g) is the generalized robustness (s ee Section 
4.2.5) and if C is the set of witnesses with Tr(W) = d, then E(g) is the random robustness [422L l423j . 
Also the distance to the set of separable states in Hilbert-Schmidt norm, the negativity and the concur- 
rence can be expressed in a similar way |l5l[ |422[ | . All these connections allow for a simple estimati on o f 



entanglement quantifiers in experiments, as any appropriately normalized witness gives a lower bound 424 1 . 
Experimentally, this has been used in Ref. [425j to probe the boundary of the set of separable states. 

Let us add that witnesses as in Eqs. (26, 27) can also be used to estimate the violation of the PPT or 
CCNR criterio n, which can further be used to estimate the concurrence or the entanglement of formation 
as in Eq. (84) |277l . |278| . Also other special witnesses can be us ed: One can relate the violation of a CHSH 
inequality for two qubits t o the entanglement of formation [4261 ] , use witnesses for the reduction criterion to 



estimate the con curr ence 1571 ] (see Section 7.2) and can also use the extended reduction map to estimate 



the concurre nce 427 1 . For special cases one can also estimate the negativity from inco mple te tomographical 
information [4281 . l429l | . A comparison of some of these methods can be found in Ref. |419l | . 

Finally, a completely different approach tries to estimate entanglement from incomplete information via 
Bayesian updating. If one assumes some a priori probability distribution p{g) for the quantum states and 
measures an observable with the result d [with a priori probability p{d) = J p(g)p(d\g)], one can update the 
probability distribution via the Bayesian rule, 

, I j\ p(d\g)p(g) . n . 

p{Qld) = p(d) ■ (172) 

This can then be itera ted. Consequently, one can obtain for the negativity a pro babil ity distribution p[N(g)] 
from incomplete data [430l | . Similar methods can be used for state tomography [43 ll ] . 



43 Clearly, for F < 1 — E the fidelity is compatible with a separable state, hence the only lower bound is the trivial bound 
E > 
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6.10. Entanglement witnesses in quantum key distribution 



Let us finally discuss another app l ication of entan glement witnesses, namely the security analysis of quan- 



tum key distribution (QKD) [432], 14331 |434j, |435[. There are many different protocols for quantum key 
distribution. They can be divided into two classes: First, in entanglement based (EB) schemes an untrusted 
source distributes an entangled state g, and Alice and Bob make measurements on it. From these measure- 
ment data they obtain the secret key, using classical post-processing of the data (e.g., key distillation or 
advantage distillation). An example is the Ekert protocol, where the source ideally distributes two-qubit 
singlet states \tp) = (|01) - |10»/V2 [|. 



In prepare & measure (P&M) schemes like the BB84 protocol [4361 14371 |438{ . Alice prepares randomly 
states from a certain set and sends them to Bob, who makes measurements. Then, using classical 

communication they obtain the secr et ke y. At first sight, these two types of protocols do not seem to be 
closely related, but, as shown in Ref. [4391 ] one can always view a P&M scheme as an EB scheme: The state 
preparation on Alice's side can be interpreted as if a bipartite entangled state 

i 

is prepared in Alice's lab, then one part is sent to Bob, and Alice measures in the computational basis \i) a- 
This is practically equivalent to an EB scheme, however, the possibilities of an eavesdropper are slightly 
more limited: In the usual EB scheme, Eve may have access to the source and can hence affect the complete 
g, while in the P&M scheme Eve can only access the part sent to Bob, i.e., the reduced state qa of Alice is 
fixed. 

What is the role of entanglement in these protocols? In the EB protocols, it is natural to ask whether the 
presence of e ntan glement is necessary for the possibility to create a secret key from the measurement data. 
Indeed, Ref. |440| showed the following result: 

Theorem. If the measurement data obtained in an EB scheme are compatible with a separable state, 
then no secret key can be extracted from that data and the QKD scheme is not secure. The same holds for 
P&M schemes, if one considers the corresponding EB scheme. 

Concerning the contrary implication, it has been shown in Ref. [44 lj that a state is entangled iff there 
exist measurements which result in probabilities which contain secret correlations. Hence, if Alice and Bob 
can verify entanglement, their correlations contain secrecy, however, as there might be bound information, 
they might not be able to distill a secret key. 

Note that the theorem requires that the entanglement can be proven with the measured data. Since many 
protocols Alice and Bob do not perform state tomography of the distributed state, this leads to the question, 
how one can prove the entanglement of a state, if only the expectation values of some observables are given. 

It is not difficult to see that if mean values (At), (Ak) of some observables are given, then entanglement 
of the underlying state can be proven, if and only if there is a witness that is a linear combination of the Ai 
that detects it. This can be seen from the fact that in the space of the mean values of the Ai the separable 
states form again a convex set and this set is characterized by linear functions of the (A%), (Ak). Therefore, 
the task remains to characterize the optimal witnesses for a given set of observables. 

As an example let us consider the Ekert protocol. There, Alice and Bob measure randomly the polarization 
in the x- and z-direction. As they can communicate via classical communication, one can assume that the 
mean values of 

Ai=a x ®t: A 2 = 1 ® a x ; A 3 = cr x <g> a x ; A 4 = a x <g> a z ; 

A 5 =a z ®l; A 6 = l®a z ; A 7 = a z <g> a z ; A 8 = a z ® a x ; (174) 



are known, while other correlations like a y ® a x are unknown ! 44 1 Using the decomposition of a witness into 
Pauli matrices, one can first prove that all entanglement witnesses W that can be measured with these 
observables fulfill W = W T = W Ta = W Tb . 



44 For the P&M protocol, the reduced density matrix of Alice is fixed, as Eve cannot access it. Theref ore, one can assume that 
also (<j y !g 1} is known. This does, however, not affect the optimal witnesses obtained later. See Ref. |442|| for a discussion. 
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Then, it has been shown in Refs. US SH that the optimal witnesses for entanglement detection in this 
class are given by 

W=|V)(^| + |VX^r B , (175) 

where is an entangled state with real coefficients. Therefore, it suffices to check these witnesses, and if 
they do not detect any entanglement, the QKD implementation is provably not secure. 

Similar ideas have then be co nside red for other QKD protocols, ofte n using th e EVM matrix [see Eq. (20)]. 



This includes the B92 protocol 442] , QKD with continuous variables [124J, 14431 ] and entanglement detection 
using Stokes parameters [l23j . 

7. Further methods of entanglement detection 

In this Section, we explain other methods for entanglement detection and quantification beside witnesses 
and Bell inequalities. We first discuss methods that have already been implemented in experiments, then we 
consider theoretical proposals that are still waiting for an implementation. 

7.1. Entanglement criteria using entropies 

Here, we describe an entanglement detection method that requires the measurement of some nonlinear 
properties of a quantum state. To start, let us recall the majorization criterion for separability [see Eq. (15)]. 
It states that for a separable bipartite state the global state is more disordered than the local (reduced) 
states. If we measure the disorder by the linear (or Tsallis) entropy S(g) = 1 — Tr(g 2 ), this implies that for 
separable states 

Tr{g\) > Trie 1 ) and Tr(g 2 B ) > Tr(g 2 ) (176) 
has to hold (80l.[8ll.l444|F F 1 For two qubits, it can be shown that this inequalities detect strictly more states 



than the CHSH inequality [444| . however, there are some entangled states that are not detected by it. 



From an experimental point of view, it is no t obv ious how Eq. (176) can be tested, as it is a nonlinear 



4<i 



function of the quantum s tate . However, in Ref. 4451 ] an elegant measurement scheme has been implemented 



that we will describe now 

Let us first discuss the determination of Tr(g 2 ). This can be written on two copies of the state as Tr(g 2 ) = 
Tr(g<& gV), where V = J^jj Ki)0'*| is the swap operator between the two copies. Using the fact that V = 
Ps — Pa = 1 — 2Pa, where Ps (Pa) is the projector onto the symmetric (antisymmetric) subspace of the 
two-copy system, we have Tr(g 2 ) = Tr[g ® g(Ps — Pa)]- Decomposing Ps and Pa into local projectors 
onto the symmetric (or antisymmetric) subspaces of Alice's or Bob's local two-copy system, one can write 
it further as 

Tr{g 2 ) = Tr([g ® g]X) with X = P^®P b s - P%®P b s - Ps® p a + PX® P A, (177) 

where Pg denotes the projector onto the symmetric space of Alice's two copies, etc., <£> denotes the tensor 
product between the two copies, and (H> denotes the tensor product between Alice and Bob. Clearly, Tr(g A ) = 
Tr([g (g) g][(-Pg — .PJ )(§>!]) can be written in a similar fashion as an expectation value on a two-copy state. 

For the experimental measurement of these quantities, two beamsplitters are required (see Fig. 19). Con- 
sider first the left beamsplitter (BS) A. If two identical photons 1 and 2 arrive at the BS the photon statistics 
of the BS will result in a coalescence of the photons, i.e., they will either be both in the output port Di or 
-D2 [13]. More precisely, if the two photons are in a symmetric polarization state, they will coalesce, while 
if they are in an antisymmetric polarization state, they will anti-coalesce and there will be one photon in 
each of the output ports. Therefore, the probability of coalescence equals the mean value of PS, and the 
probability of anti-coalescence corresponds to P%. 



45 p or generalizations of such inequalities see Refs. 184 . l8a| . 

46 For a general scheme to measure nonlinear properties without state tomography see Ref. [4461 . |44' 
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Fig. 19. Schematic view of the setup of the experiment in Ref. [445l |. Si and S2 denote the sou rces of the entangled photons, 
A and B two beamsplitters and D\, ...,Dn the four output ports. The figure is taken from Ref. |445)l . 

In the experiment, the coincidence probabilities of coalescence and anti-coalescence for Alice and Bob, 
PcaPacPca and p aa have been measured. In this language, the inequality Tr(g 2 A ) > Tr(g 2 ) is equivalent to 
Pea > Paa- Experimentally, p ac = 0.0255 ± 0.008 and p aa — 0.2585 ± 0.008 have been found, clearly violating 
the inequalities for separable states. 

There are, however, some caveat s with the interpretation of such data. First, in the generation of entangled 
photon pairs in the setup of Ref. [445| it happens that the source Si emits four photons, while S2 emits 
none. This is not a rare event (it has the same probability as the case that each source emits two photons) 
and may alter the probabilities p C c,Pac,Pca and p aa . However, these events could be taken into account by 
varying the relative phase in the pump beams of the sources <Si and $2. 

A more fundamental problem lies in the usage of two copies of the state g [448j . From a theoretical 
point there is no problem to write the separability conditions as inequalities for two copies. Experimentally, 
however, it is not clear that the two instances of the state really are the same and uncorrelated. This can not 
be verified without state tomography [2lj |. Especially, classical correlations b etween th e two systems may 
lead to a systematic error in some entanglement detection schemes, see Refs. }448l l449j for a discussion. 

A sch eme simil ar to the one described can also be used for the investigation of entanglement in optical 
lattices [450l l45l| . Namely, any fully separable state fulfills inequalities similar to Eq. (176), and also a 
transformation acting as a beamsplitter can, in principle be implemented. It should be noted, however, that 
the resulting scheme is unable to detect genuine multipartite entanglement, only full separability can be 
excluded. 



7.2. Estimating the concurrence by using several copies 

Similar ideas as in the previous subsection can be used to investigate the entanglement in experiments 
quantitatively, by giving lower bounds on the concurrence. To start, recall that the squared concurrence of 
a pure state is given by 

C 2 (|</;» = 2[l-7>(^)]. (178) 
In a similar way as in the previous section one can write therefore the concurrence as 

C 2 m) = 4<V| ® (#P2<81]|V) ® IV»> = m ® (HPl^PlW ® \i>), (179) 

where the second equality comes from the fact that P^<E>Pg acts on the global antisymmetric s ubsp ace only, 
and its expectation value should therefore vanish, if the two copies of \ip) are indeed identical [452J. 

Therefore, by measuring the (local) projector onto the antisymmetric subspace on Alice's part only, one 
can get a quantitative statement about th e entanglement contained in the state. This has been implemented 
using hyper- entangled photons [452], |453| and a proposal for the implementation using ion traps or atoms 



in cavities exists 454J. It should be noted, however, that this approach requires the assumption that the 
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generate d state is pure, which is difficult to guarantee in any experiment without doing state tomography 



449] 
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Interestingly, however, expressions similar to Eq. (179) could be used to obtain lower bounds on the 
concurrence for mixed states. As shown in Ref. [4551 ] , a lower bound on the concurrence can be obtained by 

C{pf > Tr(p(g>pV) with V = <±{Pa®Pa ~ P%®P h S )- (180) 

The difference to the previous expression in Eq. (179) is the sign of the term P|c§)Pj, which therefore, 
could be interpreted as taking account of the mixedness of the states. For an experime ntal implementation, 
still assumptions about the measurements on two copies have to be discussed [4481 14491 ] . In theory, however, 
Eq. (180) already allows to estimate the convex roof of the concurrence, and the lo wer bounds are reaso nable 
for states that are close to pure states, for investigations of the tightness see Refs. 14191 14521 14551 1456| . Also 
upper bounds on the concurrence can be achieved from measurements on two copies |457| . 
Moreover, Eq. (180) can be generalized to 

C(p)C(a) > Tr(p $ aV), (181) 

where g and a are different states fl57t This inequality has been experimentally tested using photons [458| . 
In this experiment, two two-photon states g and a have been generated via down-conversion, then the 
nonlocal observable V has been measured using a controlled-phase gate between two of the four photons. 
Using C{g) < 1 one obtains bounds on C(cr) and vice versa, however, assumptions about the fidelity of the 
phase gate have to be discussed. Furthermore, Eq. (180) shows that 

W = - 75 ^ T Tr 2 (l ® &V) = -tJ-t(o-- 1®o-b) (182) 
C{a) C(a) 

is a quantitative entanglement witness that can be used for the estimation of the concurrence. This witness 
corresponds to the reduction criterion [see Eq. (13) and the discussion after Eq. (26)], hence a lower bound 
on the concurrence can be obtained, iff the state violates this criterion. 

Finally, note that one may define a multipartite generalization of the concurrence C m as the sum of the 
bipartite concurrence over all bipartitions, 

cim = i^iY, c U\^ ( l83 ) 

bp 

where the sum goes over all 2^ — 2 bipartitions [459j . This quantity does not measure genuine multipartite 
correlations, but again expre ssions sim ilar to Eqs. (179, 180) can be derived, since bounds for each C% 



*■ ■ ■ ■ ■ * \ ' / "'bp 

are available from Eq. (180) [460L l461 |. 



7.3. Variance based criteria and nonlinear witnesse 



Entanglement witnesses and Bell inequalities can be viewed as inequalities for mean values of obscrvables, 
where a violation indicates entanglement. Therefore, it is natural to ask whether it is possible to derive 
similar conditions using the variances of observables. These conditions would then also be nonlinear in the 
state g, hence one may consider them as nonlinear witnesses. One the one hand, it might be interesting to 
use nonlinear witnesses as an improvement of a given linear witness. On the other hand, nonlinear witnesses 
are also interesting from a geometrical point of view, as they approximate the convex set of separable states 
better. 

In the literature, entanglement c riter i a using variances were fir st known for the case of continuous variables 
and infinite dimensional systems ( |l!8 . 462l 463l I46l . l465l . l466] . for reviews see Refs. dlJil, [ll[l3 ) ■ For 
the case of discrete systems, the first variance based criteria have b een derived in Refs. 255j,|467J. A very 
interesting approach was then proposed by Hofmann and Takeuchi 1071 ]. These, so-called local uncertainty 
relations (LURs) turned later out to be a very interesting tool to investigate the separability problem, so 
we will explain them now. 



47 Moreover, since one explicitly assumes that the generated state is pure, one could also just determine qa by state tomography 
on Alice's side, and then compute the concurrence or the entanglement of formation. 
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7.3.1. Local uncertainty relations 

First, recall that the variance of an observable M in the state g is given by (AM) 2 :— ((M — (M) e ) 2 ) e = 
(M 2 ) g — (M) 2 and is concave in the state |l07| . i.e., if g = pgi + (1 — p)g2 is a convex combination of two 
states, then (AM) 2 g > p(AM) 2 gi + (l-p)(AM)^. If g = \ip)(tp\ describes a pure state, the variance of M is 
zero iff \ip) is an eigenstate of M. 

Let us assume that we have observables for Alice and Bi for Bob, so that the observables in each of 
these two sets do not share a common eigenstate. Then, there must be positive numbers Ua and Ub such 
that 

n n 

J2(&Ai) 2 >U A and ^(AB,) 2 >(7 B (184) 



holds, for all states on Alice's (resp. Bob's) system. Under this condition, it was shown in Ref. [I07j that for 
separable states and for the observables Mi := Ai <E) 1 + 1 <£) Bi the inequality 

n 

Y^i&Mi) 2 , >U A + U B (185) 

i=l 

holds and violation of it implies entanglement. Physically, this may be interpreted as if the separable states 
inherit the uncertainty relations from the reduced states. 

From an experimental point of view it is interesting, that the LURs can b e vie wed as nonlinear witnesses 
improving an optimal linear witness. To see this, let us consider an example [4681 ]. For a single-qubit system 
is is easy to see that for the Pauli matrices the uncertainty relation ^2 i=x y z {AGi) 2 > 2 holds. Defining 
Mi = <Tj <g> 1 + 1 <g> <Ti this yields the LUR ^2 i=x y z (AMi) 2 > 4. Writing down this explicitly, one obtains 
the condition 

(1 <g> 1 + a x <g) g x + G y <g) g v + a z <g> g z ) - i ^ (^i <S> 1 + 1 ® G t ) 2 > 0, (186) 

i— x,y ,2 

which has to hold for all separable states. This is a quite remarkable inequality for the following reason: 
The first part, which is linear in the expectation values is known to be an optimal entanglement witness [see 
Eq. (123)]. From this witness some quadratic terms are subtracted. Note that an experimental measurement 
of these terms would require just the same measurements as the measurements that are required anyway 
for the measurement of the witness in Eq. (123) . T hus, in this case, the LUR can be viewed as a nonlinear 
witness that improves an optimal linear witness F^l 

In general, the e ntanglement criteria resulting from the LURs are very strong: they can detec t bound 
entanglement 469], indeed, it can be shown that they are stronger than the CCNR criterion [llpf 



-1!) 



Finally, it turned out, that they are equivalent to the covariance matrix criterion described a bove [KM 1 114 ] 



(see Eq. 18). Also, statements similar to the LURs can be extended to nonl ocal obse rvables 108], or other 
formulations of the u ncert ainty principle, like entropic uncertainty relations [470ll47ll or the Landau-Pollak 



uncertainty relations [4721 ]. The local uncertainty relations and the varia nce based criteria c an also be used 
to give lower bounds on entanglement measures, similarly to Section 6.9 [116L 14731 . 1474L |475| . Conc ernin g an 
experimental implementation, some of the variance based criteria have been implemented in Ref. [476| . 



7.3.2. Nonlinear entanglement witnesses 

The LURs show that in some cases it is possible to construct nonlinear witnesses. However, it is easy to 
see that the recipe of the LURs can not be applied to arbitrary witnesses. Therefore, the question arises 
whether it is possible to find a general recipe to improve a witness by some nonlinear terms. 



Indeed, as shown in Ref. [4771 ] . for bipartite witnesses this is the case. The recipe is the following: Let us 



consider the standard witness from the PPT criterion 

W=|^)(^| Tb (187) 



48 It is further interesting that for any pure product state \ip) = \a)\b) equality holds in Eq. (186), which shows that the 
nonlinear witness is curved around the whole set of separable states. 



49 



This can be easily seen from the fact that LURs always allow to improve the witnesses for the CCNR criterion, see Eq. (27). 
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Let us define X = with an arbitrary \<p) and consider the functional 



He) = (W) - —r(x T -)((x T ^), (188) 

s((j>) 

where s(4>) denoted the maximal squared Schmidt coefficient of \<p), which is a bound on the maximal 
overlap between |</>) and product states. Then we have for a pure product state |?7) = \a)\b) that F(\ri){ri\) = 
|(?7|V>)| 2 [1 — | (^l^) \ 2 /s(4>)] > 0, where |^) = |a)|6*) is the vector corresponding to the projector \rj)(r]\ TB = 
— \ab*)(ab*\. By convexity, this positivity holds also for all separable states, proving that F{g) is 
indeed a nonlinear entanglement witness, which improves the witness W. 

There are also some other constructions of nonlinear witnesses bases on the same idea [473, Ezi and due 
to the Choi-Jamiolkowski isomorphism [see Section 2.5.3 and Eq. (35)] all these methods can be extended 
to arbitrary bipartite witnesses. This has also some co nsequen ces for the geometry of separable states, as 



it proves that the set of separable states has no facets [478jl 50 1 Also, one can define sequences of nonlinear 
witnesses by an iteration procedure, such that the nonlinear witnesses become stronger and stronger in each 
step and fi nally detect everything that is detected by the corresponding positive, but not completely positive 
map [122] . 

At the moment, there are not much results known concerning nonlinear witnesses for multipartite entan- 
glement. The method from above can not be extended to witnesses for genuine multipartite entangle ment . 



Some nonlinear witnesses for multipartite entanglement have already been proposed in 2002 by Uffink 255 1 
and recently some crite ria have been developed, which can also be viewed as some nonlinear conditions on 
the state (see Ref. |l90j and Section 3.5). 

7.4. Estimation of positive maps 

Given the importance of positive maps for the theoretical characterization of entanglement, one may 
wonder whether it is possible to use them directly for entanglement detection. Indeed, several proposals 
for this task exist. Although their implementation is beyond the current experimental capabilities, they 
offer interesting insights in measuring certain nonlinear properties of the density matrix without doing state 
tomography. 

As the transposition is not a com pletely positive map, it is not possible to perform it directly in a 
lab. However, as shown in Ref. [48C| one can circumvent this problem by a so-called structural physical 
approximation (SPA) of the transposition: if a sufficient amount of white noise is added, the total map 
becomes completely positive. Especially, it was shown in Ref. |480l ] that in a d x d system the map 

a(q) = ^TT lA ® 1b + d^Ti gTB (189) 



is completely positive and can be implemented. SPAs of other positive maps have been discussed in Ref. 481] . 
The noise term d/(d 3 + 1)-1a®1b only shifts the spectrum of the partial transposition q Tb , hence it remains 
to measure the eigenvalues of a = A(g) to evaluate the PPT criterion. 

To do so, it suffic es to mea s ure T r(a), Tr(a 2 ), ...,Tr(a d ), as from these d 2 quantities all d 2 eigenvalues 



can be determined [446l . I447L |482|. The quantity Trip ) can be determined similarly as in Section 7.1 
by measuring a multiple of the shift operator on a fc-fold copy. Alternatively, one may use the spectrum 
estimation technique from |483| . Note that knowing the spectrum of a, one can then determine the negativity 
as an entanglement mea sure. 

As shown in Ref. |484| one can determine the concurrence in a similar way: As explained in Section 4.2.3 
the concurrence for two qubits can be computed from the eigenvalues of gg with 

g = (a y <X> a y )g T {a v ® a y ). (190) 



50 A face of a convex set is a hyperplane at the boundary, and a facet is a face of maxi mal d imension. For example, a three- 
dimensional cube has its six surfaces as facets, while the edges are faces, but not facets 479] . Therefore, the set of separable 
states is not as simple as in Fig. 1(b). 
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To determine Tr[(gg) k ] requires 2k copies, where on half of them first the transposition is implemented via 
the SPA, and then the unitary transformation (a y ® a y ) is performed. Then, the eigenvalues of gg can be 
determined similar as before. 

The previous proposals try to ch aract erize entanglement by implementing an interesting positive map 
physically. However, as noted in Ref. [4851 ] this is not mandatory. The fc-th moments Tr\( p TB ) k ] are invariants 
under local unitary transformations, which can be measured using qua ntum net works [4861 ] . They can also 
be measured directly by making measurements on k copies of t he st ate [48 7l |488| . Similar ideas can be used 
to determine the concurrence or the 3-tangle of mixed states 14891 or to evaluate the separability criteria 
given by arbitrary positive, but not completely positive maps 490]. 

Finally, for the case of two qubits, similar met hods can be used to decide separability using a single, 



nonlinear entanglement witness. As shown in Ref. 49 1| . for two qubits one has that 



g is entangled det(g TB ) < 0. (191) 

This equivalence reflects the fact that for two-qubits the partially transposed state can have maximally one 
negative eigenvalue |47j |. The quantity det(g TB ) is a fourth order p olyno mial of the matrix elements of g and 
can hence be measured by a single witness on a fourfold copy of g 49 1| . Similarly, for rota tionally invariant 
states on a 2 x <i-system, the PPT criterion is necessary and s ufficient for separability [492{ , and if d is even, 
Eq. (191) can again be used for entanglement detection [4931 ]. 



8. Entanglement detection with collective measurements 



In the following Section we discuss how to detect entanglement in a many-qubit physical system if the 
qubits cannot be individually accessed. Fo r su ch systems entanglement criteria are needed that are based on 
the measurement of collective observables l 51 1 This includes spin squeezing inequalities as well as measure- 
ments of the magnetic susceptibility or the energy in spin models. We also discuss entanglement detection in 
optical lattices of cold atoms. In future, when larger and larger coherent quantum systems will be prepared, 
such criteria will surely gain large importance. It should be noted, however, that several criteria in this 
Section do not detect genuine multipartite entanglement, they only rule out full separability. 



8.1. Spin squeezing inequalities 



In this Section we describe entanglement detection with collective spin-i observables, which is relevant 
to multi-qubit systems in which the qubits cannot be accessed individually. Even for systems in which the 
qubits can be accessed individually, such criteria can be advantageous since they need few measurements. 

In a multi-qubit system, the quantities that can be measured collectively are the components of the 
collective angular momentum 

i N 

Jr-=lY,°i k) > ( 192 ) 

fe=i 

where I = x,y,z and Oj are the Pauli spin matrices, and their moments (<//") where m = 2,3,4,... 
This includes also the variances (A J*,) 2 = ( J|) — (Jk) 2 ■ Clearly, measuring only the expectation values of 
collective observables is not enough for entanglement detection: Such measurements do not give information 
on entanglement between the spins. The minimum requirement for detecting entanglement is measuring the 
first and second order moments of collective angular momenta [400] . 



8.1.1. Spin squeezing and entanglement. 

The first criterion based on such ideas was the spin squeezing criterion. On the one hand, Ref. [494| defined 
spin squeezing in analogy with squeezing in quantum optics. Spin squeezing, according to the definition of 

51 Note that in the literature, the notion of a "collective measurement" is sometimes differently used to describe measurements 
on several copies (as in Section 7.1). 



Rcf. 494], means the following: In general, the variances of the angular momentum components are bounded 



by the uncertainty relation 

{^Jyf{AJ Z f>\\{J X )\\ (193) 

If (AJ 2 ) 2 is smaller than the standard quantum limit \(J x )\/2 then the state is called spin squeezed ! 52 1 In 
practice this means that a spin squeezed state has a large mean sp in and a small variance in a direction 
orthogonal to the mean spin. On the other hand, in Refs. [495l |496] |. spin squeezing was defined from the 
point of view of spectroscopy. States fulfilling this definition are useful for reducing spectroscopic noise or 
increasing the accuracy of atomic clocks . 



53 



It has already been pointed out in Rcf. 494] that the collective spin can be squeezed due to quantum 
correlations between the particles. Indee d, sp in squeezing can be shown to be connected to entanglement: If 
an AT-qubit state violates the inequality [493] 

(AJz)2 >L (194) 



(J*) 2 + (Jy) ~ N 



then the state is entangled, i.e., not fully separable. The de gree of violation can also be used to characterize 
how strong the entanglement is [5001 ] . It was found in Ref. [50l| that if a symmetric state violates Eq. (194) 
then its two-qubit reduced density matrix is entangled, that is, the state is two-qubit entangled. However, in 
gene ral, the violation of Eq. (194) does not indicate two-qubit entanglement, as one can find counterexamples 

A typical spin squeezing experiment starts with a fully polarized sample |^') init := \^)^ N ■ Such a state 
saturates Eq. (194). During the entangling dynamics the larg e spin decreases only slightly, while the un- 
certainty in an orthogonal direction decreases considerably ! 54 1 Thus, the resulting state violates Eq. ( 194) . 
Entanglement can be generated by direct interaction between the particles, as in the proposal of Ref. 499] . 
In another scheme, the atomic ensemble interacts with light [503j . The light is then measured, projecting 
the atoms into an entangled state. For an ensemble strongly polarized into the x direction one has (J x ) ~ 
and (J y ) ss 0. Thus, the entanglement condition Eq. (194) can be reformulated as (AJ Z ) 2 > ±1. 

8.1.2. Singlet criterion 

After the discovery of the criterion Eq. (194), it was realized that other generalized spin squeezing criteria 
can also be constructed that detect various types of entangled stat es and still need only collective measure- 
ments. Such an entanglement condition was presented in Ref. [4001 ] . It contains the variances of all the three 
angular momentum components: For separable states 

(AJ x f + (AJ y ) 2 + (AJ Z ) 2 > ^ (195) 



holds. States violating this condition have small angular momentum: Eq. (195 ) detects entanglement in the 
vicinity of many-qubit singlet states, which were considered in Sec. 3.4.4 1 55 1 Such states give zero for the 
left hand side of Eq. (195). One of such states is the chain of two-qubit singlet states. Another example is 
the ground state of the isotropic anti-ferromagnetic Heisenberg chain. Interestingly, all pure product states 
saturate the inequality. Eq. (195) has a good noise tolerance, it detects a singlet state mixed with white noise 
as entangled if p no isc < f ■ Again, violation of Eq. (195) does not imply two-qubit entanglement. Moreover, 
it can de tect PPT entangled states (bound entangled states) in the thermal states of several spin models 



52 Note that for Eq. (193) we assume that the mean spin points to the x direction. Without this requirement, even a pure state 
in which all spins point into the same direction could violate Eq. (193) with an appropriate choice of the x, y, and z coordinate 



axes [494 

53 The definitions of sp in squeezing in Ref. I494H and in Refs. |495| . I496H are slightly different. For a discussion about that, see 
Sec. IV. A in Ref. 497]. For ananalysis concerning the two-qubit case, see Ref. [498l |. 

54 For an exhaustive survey of the topic and a list of references, see the excellent review in Ref. |497l| . 

55 Note that the c riteri on Eq. (195) can straightforwardly generalized for TV spin-j particles. In this case on the right hand side 
there must be Nj |400|| . 
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In addition, the violat ion of the criterion gives information about the number of spins that are unentangled 
with the rest 502L 15041 ] : Let us consider a pure state for which the first M qubits are not entangled with 
other qubits while the rest of the qubits are entangled with each other \9) = {®^ =1 \^k)) ® W )m i i n- F° r 
such a state, based on the theory of entanglement detection with uncertainties, one has [1081 ] 

(AJ,) 2 + (AJ y ) 2 + (AJ Z ) 2 > M (196) 

If a mixed state g := J2kP k \^ k )(^ k \ violates Eq. (196) then at least one of the components |^ r fc)(^fc| must 
have M or more spins that are entangled with other spins. If the left-hand side of Eq. (196) is smaller than 
i then the state cannot be created by mixing states that have one or more unentangled spins. 



8.1.3. Inequalities detecting entanglement close to Dicke states 

A third type of inequality for entanglement detection appeared in Ref . |213| : For separable states 



Uz) < 



N 



N 
~2 



(197) 



holds. This inequality detects entanglement close to the iV-qubit Dicke state with two excitations. It detects 
such a state as entangled if for the added white noise p n oise < jr- Thus, the robustness to noise is decreasing 
with N. Moreover, the violation of this inequality implies the presence of two-qubit entanglement. This can 
be seen as follows: The average two-qubit density matrix is 

1 x - 



£>av2 



N{N — 1) ^ 



Qkh 



(198) 



k^l 



Here Qki represents the two-qubit reduced density matrix obtained from the original one after tracing out 
all qubits but qubits k and I. Eq. (197) can be re-expressed with q w i as 



(4 X) 4 2) 



WV. < i. 



(199) 

Eq. (199) can be violated only if g av2 is entangled [see Eq. (151)]. 

Criteria, similar to Eq. (197) can also be used to detect genuine multi-qubit entanglement. Such criteria 
for N = 3 and 4 qubits, respectively, are 



(Ji 



K 2 ><2- 

7 



V$/2 « 3.12 for N = 3, 



<^ 2 ) + <J|)<x + V3« 5.23 for N 



(200) 
(201) 



Entanglement was detected base d on Eq. (197) around a four-qubit Dicke state with two excitations in a 
four-qubit photonic system [2151 ] . Although in such a system the qubits are individually accessible, using 
such a criterion was still advantageous since only two measurement settings (i.e., the {<J Xl o~ x ,a x , ...,cr x } and 
{<7y, o~ v , o~ v , o~ v } settings ) were needed for the detection. The experiment will be explained in Section 8.2. 



8.1.4. Detecting t wo- and three-q ubit entanglement 

Finally, in Refs. [iOTHHoBl |506| spin squeezing criteria were constructed detecting two-qubit entanglement 
of the reduced two-qubit state. It was proven that if 



> 



N 



{4) + { J l)-i 



+ (n - 1) 2 ( j x y 



is violated then the state is two-qubit entangled. This criterion can be simplified to 
1 4(J Z ) 2 ^ 4(AJ 2 ) 2 



> 



(202) 



(203) 



N 2 ~ N 

for symmetric states. A symmetric state is two-qubit entangled iff it violates Eq. (203) for some choice of 
the z-axisl 56 ! Moreover, criteria detecting genuine three-qubit entanglement of the W and the GHZ classes, 



Ref. [507] conjectures a formula with collective obscrvables for the two-qubit concurrence of symmetric states. 
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respective ly, w ere also presented that need the measurement of third order moments of the collective angular 
momenta [40 7j . T he p revious ideas have already been applied to investigate the data of a ion trap experiment 
(see Section 6.8) [4071 ]. Such criteria a re especia lly applicable to W states, as W states maximize two-qubit 
concurrence among symmetric states |l76l . Il84| . 



8.1.5. Optimal spin squeezing inequalities 

After presenting several generalized spin s quee zing inequalities the questions arises: can one find a unifying 
framework for all these inequalities? In Ref. [l38j this has been done. It turned out, that the first and second 
order moments of collective angular momenta for fully separable states fulfill simple constraints, namely 

(J 2 x ) + (J 2 y ) + (Jl)<^^, (204a) 
(AJ X ) 2 + (AJ y ) 2 + (AJ Z ) 2 > ^, (204b) 

(4) + (J?) y < (N- l)(AJ m ) 2 , (204c) 

(N 1) [(AJ,) 2 + (AJ,) 2 ] > (J 2 m ) + N{N - 2) . (204d) 

While Eq. (204a) is valid for all quantum states, Eqs. (204. b-d) can be violated. Any state violating them is 
entangled. It was also proven, that this set of generalized spin squeezing inequalities is complete: For certain 
cases (e.g., (J x ) — (J y ) — (J z ) — and N is even) or for large N, a state that does not violate Eqs. (204b-d) 
cannot be detected by any other entanglement condition based only on the expectation values of the first 
and second moments of Jfc. 

For any value of J := ((J x ), (Jy), (Jz)) these eight inequalities define a polytope in the three-dimensional 
((J 2 ), (Jy), (J 2 ))-space. Separable states lie inside this polytope. If one of the inequalities is violated, then 
the state is on the outside and is hence entangled. For the case J = and N = 6 the polytope is depicted 
in Fig. 20. Such a polytope is completely characterized by its extremal points. They are given by 



A x :— 
B~ := 



(205) 



where K := (N — 1)/N. The points A y / Z and B y / Z can be obtained from these by permuting the coordinates. 

Interestingly, violation of these spin squeezing inequalities does not necessarily imply entanglement of the 
reduced two-qubit states. As all the expectation values of first and second moments of Jk can be evaluated 
knowing the average two-qubit reduced density matrix g aV 2 defined in Eq. (198), one could naively expect 
that these criteria detect the entanglement of the reduced two-qubit density matrix. However, somehow 
surprisingly these criteria are also able to detect entangled states for which the reduced states are not 
entangled. 

This phenomenon is related to the representability problem |508l ]. i.e., to the problem of finding multipar- 
tite quantum states that have a given set of states as reduced states. For a given g aV 2, one has to answer 
the question whether there is a separable N— qubit state that has g aV 2 as the average two-qubit reduced 
state? If the answer is no then the system is entangled. Interestingly, it is also possible that g aV 2 is separable, 
however, there is not an ./V-qubit separable quantum state that has it as a reduced state. In this case we can 
conclude that the system is an entangled state even if p aV 2 is separable. If g a v2 is separable and of the form 

£av2 = ^PkQk <S> Qk (206) 

k 

then there always can be found an A^-qubit separable state that has this state as its reduced state qn = 
12kP k @f N ■ Note that separable states in the symmetric subspace can always be written in the form of 
Eq. (206). On the other hand, if p aV 2 is not symmetric then it is possible that it is separable, however, there 
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Fig. 20. The polytope of separable states corresponding 
to Eq. (204) for N = 6 and for J = 0. The point S 
corresponds to a man y body singlet state. The Figure is 
taken from Ref. KM- 




Fig. 21. Experimental setup for the analysis of the 
four-photon polarization-entangled state I.D24). The 
state is detected after the symmetric distribution of four 
photons onto the spatial modes a, b, c and d via non-po- 
larizing beam splitters (BS). The Figure is taken from 
Ref. [215II. 



is not an N— qubit symmetric separable state that has it as its reduced two-qubit state. One can write down 
the condition Eq. (204b) with the operator expectation values of g aV 2 as 



N 



(a x ® a x + <j y ® <j y + a z ® a z ) eav2 > — - ^ ^ (a k ®\)\ 



1 



k—x,y,z 



N - 1 



(207) 



If this condition is violated then Eq. (204b) detects the state as entangled. However, Eq. (20 7) is not a 
condition for detecting the entanglement of f? aV 2- That is, there are separable g aV 2 that violate it ! 57 1 



8.2. Experimental observation of a four-qubit symmetric Dicke state 

In this Section we explain an experim ent a iming to create a four-photon symmetric Dicke state with two 
excitations, I-D2.4) in a photonic system [2151 ]. On the one hand, the experiment is interesting since the state 
I-D2.4) is highly entangled while the experimental setup is quite simple. On the other hand, the generalized 
spin squeezing criteria can be used for the detection of entanglement with few measurement settings. 

The Dicke state |i?2,4) has the form 

\D 2 4) = -j=(\HHVV) + \HVHV) + \VHHV) + \HVVH) + \VHVH) + \VVHH)) 7 
V6 

where \H) and \ V) denote linear horizontal (H) and vertical (V) polarization of a photon in the four spatial 
modes, denoted by a, 6, c, and d. This is the symmetrization of a product state with two H and two V 
photons. 

The setup can be seen in Fig. 21. The four photons are created by the second order emission of collinear 
type-II parametric down-conversion. The BBO crystal is used for the down-conversion, while the UV-mirrors 
separate the down-converted photons from the pump beam. As a result, in the single-mode fiber four photons 
appear, two in state \V), the other two in state \H). Now, the four photons must be distributed into four 
spatial modes. This is done by beam splitters. Birefringent effects of these non-polarizing beam splitters 
is compensated by pairs of perpendicularly oriented YVO4 crystals in all the four modes. Note that the 



57 For example, for the T = ground state of the Hamiltonian H = + Jy + J% a simple calculation shows that 
(u x ® a x + cT y (g> (j y + a z ® cr z ) eav2 = — 4^ JV 3 +1 - ) and (a k ® l}g av2 = 0, but p av 2 can be shown to be unentangled for N > 4 for 
even N [Hoi . 
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Dicke state can also be p repared with a different setup, which flexibly allows to create an entire family of 



multi-photon states [217 1 



In the experiment, first the density matrix was determined by measuring all correlations of the type 
(<7fc <8> ctj ® Cm ® c„) where fc, l,ro,n6 {0, a;, y, z}. For that 3 4 = 81 measurement settings are needed. This 
makes it possible to compute the fidelity of the prepared state: F exp = 0.844 ±0.008. Note that the fidelity is 
high since the generation of the state was quite straightforward with parametric down-conversion and post- 
selection and no additional gates were needed. In comparison, for creating a cluster state, it is necessary to 
realize a phase gate . The larger complexity of the quantum circuit results in a decrease of the fidelity of the 
prepared state [232|. 

The fidelity was higher than the bound for biseparable states, 0.6 [see the witness Eq. (75)]. Thus, the 
fidelity computed from the density matrix already signals the presence of genuine multipartite entanglement. 
Instead of complete tomog raph y, 21 settings would also be enough to get the fidelity and detect genuine 
multipartite entanglement [2151 ]. However, in the experiment a more efficient criterion, namely Eq. (201) 
was also tested, which needed only two settings. For the left hand side of Eq. (201) a value of 5.58 ± 0.02 
was obtained that clearly demonstrated genuine multipartite entanglement. 

In the second part of the experiment, one of the photons was measured. This leads to various symmetric 
states, in particular states of the type 

\1>) = a\W 3 ) + 0\W 3 ), (208) 

where \W 3 ) is the three-qubit W state and \W 3 ) is a state obtained from \W 3 ) by inverting all the qubits. In 
the experiment, the criterion Eq. (200) was used to detect the genuine multipartite entanglement of these 
states. After the loss of one of the photons, from the \D 2 ^) state we get a mixture of \W 3 ) and \Wz) states. 
This state is also genuine multipartite entangled, that indicates the persistency of the entanglement of the 
|D2,4) state. Criterion Eq. (200) was then used to detect entanglement of this state. 

8.3. The Hamiltonian as a witness 



In this and the following Section we discuss how entanglement can be d e tect e d by measuring fundamental 
quantities of a physical system, e.g., the energy of a spin chain [Eol H3 

1 5 1 IB . This allows to connect 

thermodynamical quantities of condensed matter systems to their entanglement properties [2^, |512| . We will 
also discuss how these ideas can be extended to detect multipartite entanglement and how other fundamental 
observables beside energy can also be used. It should be noted, however, that in these approaches the 
conclusion that a given state was entangled very often depends on the assumption that the considered spin 
model indeed describes the physical system. 

There are two basic approaches in entanglement detection: (i) One would like to detect entanglement in the 
vicinity of a given quantum state. Then some entanglement condition, for instance an entanglement witness, 
is designed for this aim. (ii) It is also possible to proceed in the opposite way and look for a fundamental 
observable of a physical system that is easy to measure. Then one can examine what quantum states can 
easily be detected by such an observable. When following the second path, one constructs witness operators 
of the form 

W O -=0- inf [(* |0|tf }] , (209) 

where S is the set of separable states, "inf" denotes infimum, and O is a fundamental quantum operator 
of a spin system that is easy to measure. Note that in the general case inf^ e s(^E'|0|^') is difficult, if not 
impossible, to compute. 

A natural observable for the second approach is the Hamiltonian of the system. First of all, it is typically 
easy to measure. In some systems it can be measured directly. In other systems the expectation value of 
the Hamiltonian can be measured indirectly by two-body correlation measurements since the Hamiltonian is 
typically the sum of few two-body correlation terms. Second, it turns out that it is simple to find the minimum 
for separable states for many spin model Hamiltonians. This makes it possible to investigate entanglement 
in systems in thermal equilibrium, which appear many areas of physics. In thermal equilibrium the state of 
the system is given by 
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Fig. 22. (a) Heisenberg chain of 8 spins. Nearest-neighbor entanglement (quantified by the entanglement of formation) as a 
function of magnetic field B and temperature T. (b) The same for an Ising spin chain. Here kg is the Boltzmann constant, J 
and J x are coupling constants. Lig ht color indicates the region where entanglement is detected by Hamiltonian as a witness. 
The Figure is taken from Ref. [51011 . 



qt oc exp(— H/kgT), 



(210) 



where T is the temperature and fcs is the Boltzmann constant. For simplicity we will set fee = 1- Using 
Eq. (209) a temperature bound, T E , can be found such that when T < T E then the system is shown to be 
entangled. 

To give a simple example, consider the one-dimensional Heisenberg Hamiltonian on N spins with periodic 
boundary conditions, 



(211) 



With an argument as in Eq. (151) one can directly see that for a fully separable state (Hjj) > —N has to hold. 
As the ground state energy for large N for the Heisenberg model equals Eq = — [41n(2) — 1]N ss — 1.773A 
the thermal state at low temperatures (namely k E T < 3.18) has to be entangled [193J. 

Similar ideas have been applied for the detection of entanglement in various spin mode ls such as the XY- , 
Heisenberg- and Ising spin lattices with an external field in different various dimensions [HI EH, [Ho3, [509, 
Hill [Bill . l5ll l515l IBia l517j . An example is given in Fig. 22. The figures show the entanglement of formation 
for the reduced two-qubit state as a function of the external field and the temperature. Light color indicates 
that the energy of the system is smaller than the minimum for separable states thus the Hamiltonian as a 
witness detects it as entangled. Note that if the entanglement is not very small, the Hamiltonian as a witness 
detects the thermal state as entangled. Thus, in such systems measurement of a single observable, i.e., the 
energy, seems to be sufficient for entanglement detection. Moreover, for systems in thermal equilibrium it 
is po ssible to find witnesses that are different from the Hamiltonian and detect all entangled thermal states 
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Beside entanglement in general, one can also use similar ideas for detecting various forms of multipartite 
entanglement. As explained in Section 3.3 [E q. (53)] one can ask whether multipartite entanglement is 



necessary to form a given state [190l . Il92l Il93l |. Then, a state \ip) contains only /c-party entanglement (or is 



fc-producible) if one can write the state as a tensor product l^) = \(f>i} ® \4>2) ® ■■■ ® \<t>m)i where the 
states \<j>i) are states on maximally fc-qubits. In this definition, a two-producible state does not contain any 
multipartite entanglement, since it suffices to generate the two-qubit states \<fii) to arrive at the state l^). In 
addition, one says that a state contains genuine fc-party entanglement if it is not producible by (k — l)-party 
entanglement. This definition can be extended to mixed states via convex combinations. 

For spin chains of macroscopic size, it is in general very difficult to prove that the total state is genuine 
./V-partite entangled via energy measurements. This is due to the fact that the notion of genuine iV-partite 
entanglement is extremely sensitive to the properties of a single qubit. Indeed, in order to prove genuine 
multipartite entanglement, one has to exclude the possibility, that one single qubit can be separated from 
the remaining N — 1 qubits. However, multipartite entanglement in the reduced states of small numbers 
of qubits can easily be detected. Moreover, if the r educed state is mul tipartite entangled then the state is 
not two-producible. Based on these ideas, in Refs. [192 . EM 15181 l519j energy bounds for two-, three- and 
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four-producibility are presented for various spin models such as the Heisenberg- and the XY model, and also 
for two-dimensional models with frustration. 



8.4. Entanglement detection with susceptibility measurements 

Finally, let us consider the magnetic susceptibility as another thermodynamic quantity. In particula r, th e 
generalized spin squeezing inequality Eq. (195) can be reformulated using susceptibility measurements 520]. 
Let us consider a system with a Hamiltonian H and add a magnetic interaction Hj := ^2i =x y z BiJi to it, 
where B is the strength of the magnetic field. It has been shown in Ref. |520l ] that if we define the magnetic 
susceptibility along the directions I — x, y, z as 

( d{Ji) 
V dBt 

and if H commutes with J k for k = x,y, z then for separable states 

Xx + Xy + Xz > 2fcy. ( 213 ) 

Here it is assumed that the system is in thermal equilibrium. As one can directly check [using \i = 'Wri^Ji) 2 ] 
Eq. (213) is equivalent to Eq. (195), giving the spin squeezing inequality a new physical interpretation. 
There have be en ex periments for detecting entanglement with susceptibility measurements in solid state 
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(212) 

B=0 



systems. In Ref. 5211 ] . it is explained that the magnetic susceptibility of copper nitrate measured in 1963 
shows the presence of entanglement. Recently, experiments in low-dimensional spin systems have been carried 
out with pyroborate MgMnB205 and the warwickite MgTiOBOa, systems with spin 5/2 and 1/2, respectively 

E q. (21 3) has also been used for looking for entanglement in small molecular spin-clusters theoretically. 
Ref. 523] determined the bound temperature below which the system is en tang led. The four spins interact 
via nearest-neighbor, diagonal, and four-spin Heisenberg interactions. Ref. 524] considered thermal entan- 
glement in the nanotubular system Na2Vs07 and determined also the bound temperature for entanglement. 

8.5. Entanglement detection in optical lattices 



r,.s 



Optical lattices of cold atoms are important examples for large-scale quantum information processing 
In particular, in a three-dimensional lattice of 10 5 two-state bosonic atoms an array of one-dimensional 



cluster states has been generated [13[. The entangling-disentangling dynamics, which created the cluster 
state and then restored the initial product state, was observed through the appearance and disappearance 
of the interference fringes. In this experiment the lattice sites could not be addressed individually, thus 
entanglement detection in such systems is a relevant example of entanglement detection with collective 
measurements. 

The simplest model of such a one-dimensional array of trapped two-state bosonic atoms is the following. 
Let ak and b k be the destruction operators corresponding to the two species at site k. (Here, s pecie s "a" 



and "b" mean atoms in states \a) and respectively.) A typical Hamiltonian for this system is 527] 
^lattice := - J a ^2 i a k a l+i + a l a k+i) - Jb ^ (kfc & l+i + b l b k+i) ~ (a k bl + a[b k ) 

+ \Ua a a\a\a k a k + \U bb b\b\b k b k + U a ba k a\b k b\, (214) 

where J k are the tunnel couplings between the sites, £1 is the Rabi frequency for the transition \a) — > \b), 
and U k are the on-site interactions between atoms. The term a k a k+1 + a k a k+ i denotes the tunneling of an 

atom in state \a) between sites k and k + 1. Moreover, a k b k + a k b k denotes a process in which an atom at 
site k is turning from state \a) to state \b) and vice versa. This can be the result of manipulating atoms 



'For a review on quantum information processing in optical lattices of cold atoms see Refs. [525, 52 ! 
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with a laser pulse. The sign of the on-site interactions U a and Ub between atoms of the same specie control 
whether atoms favor to be at the same site or they tend to be distributed equally over the lattice. The sign 
of U a b determines whether atoms with different internal states favor to be at the same site. 

In the realization of Ref. [13j |. the atoms cannot be individually addressed. In the experiment, after 
producing the desired state, one lets the cloud expand to many times of its original size. Atoms with a 
large momentum can get further from their original position than atoms with a small momentum. Thus, 
from the distribution of the particles after the expansion, one can obtain information on the distribution in 
momentum space before the expansion. Knowing the distribution in the momentum space, one can obtain 
the distribution in real space by a Fourier transform. This way, in principle, one can obtain for the prepared 
state such quantities as 

{QD := 5>t,<W*> and (Q?) := Y,^l^+^ m+x b m ) ■ (215) 

m rn 

In one of the realizations of quantum information processing, most of the lattice sites have a single atom. 
The qubit is encoded in the internal states of the atoms. Schemes can aim to detect specific states in the 
experiment: 

(i) Cluster states. In a lattice of such two-state atoms, typically the atoms can interact with the near- 
est neighbors in a one-dimensional chain. While in some systems such an interaction is part of the 
Hamiltonian 528], in the experiment of Ref. fl3l it was efficiently engineered by applying two different 
optical lattices for trapping the two species [5291 ]. Thus, cluster states can be realized, since they could 
be created from product states with a nearest-neighbor Ising dynamics. However, first and second 
mome nts o f the collective observables Jk are not sufficient for detecting the entanglement of cluster 
state s [400j . One proposal in this case is to use dynamics before the detection of collective observables 
400] . In this way, a simple scheme can be constructed that can be used to study the decohering cluster 
state and find a good lower bound for its entanglement lifetime. However, since now the dynamics is 
part of the entanglement detection, its accuracy must be included (or assumed) in the calculations. 

(ii) Single t sta tes. Instead of cluster states, two-particle singlets can be generated in a lattice of double- 
wells [530j |. Such states can be readily detected by the criter ion Eq. (204b), based on second order 
moments, and its generalization for multi-qudit systems 400]. This is also true for various proposals 
for real izing multi-particle singlet states as ground states of Heisenberg spin lattices with spins \ or 
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Entanglement detections schemes not designed for particular states are also available. Such proposals are as 
follows: 



(i) In Ref. [5321 ]. an approach is described using two copies of a multi-qubit pure state to compute a 
multi-particle entanglement measure based on the single-qubit purities. The two copies can be in two 
one-dim ensio nal optical lattices that interact with each other during the entanglement detection. 

(ii) In Ref. [450l ] , a method is suggested that works if one possesses two copies of the mixed multi-qubit 
state of the lattice (see also Section 7.1). It is based on the fact that for separable states the purity of 
reduced states is larger than the purity of the state. If for a quantum state this is not the case, then 
it is ent angl ed. 

(iii) In Ref. [533j , it is shown how to obtain a bound for the fidelity of the two-qubit states with respect to 
some maximally entangled state, measuring correlations between the momentum distributions, that 
is, measuring (Q* ) for k,l = a,b. Knowing the fidelity with respect to a maximally entangled state 
allows to obtain a bound on the two-qubit entanglement of formation. 

In the other realization of quantum information processing one does not employ the internal stat es fo r 



storing information but the presence or absence of atoms at the lattice sites encodes the single qubits [534 1 . 
That is |0) corresponds to an empty site, while |1) describes a site with a single atom. The system can 
be described by the Bose-Hubbard model, that is, Eq. (214) with only specie "a". In the limit of very 
strong repulsive interaction between the particles, there is at most a single particle per site. Entanglement 
conditions for such systems are the following: 

(i) In Ref. [504| . a simple entanglement condition is described that needs the measurement of two-site 
correlations of the type {aka\ +1 + a^a^+i), where df. is the destruction operator corresponding to the 
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site k. This criterion detects entanglement close to the ground state of the system, 
(ii) In Ref . [533l | , an approach is presented that exploits particle number superselection rules. These prohibit 
to have single lattice site states that are the superposition of states with different particle numbers. 
Under such conditions, all states that can be prepared locally are diagonal in the basis of |0) and 
|1). Thus, if a state is nondiagonal, then it implies the presence of a resource, which is, however, 
less powerful than entanglement defined independently from superselection rules. Nevertheless, the 
state detected cannot be prepared locally with atoms, thus if the state is nondiagonal, the presence of 
quantum correlations are demonstrated. Any non-flat moment um distribution indicates entanglement 
in this sense, which has been observed in several experiments [533]. Moreover, knowing (Q x ) makes it 
possible to obtain a lower bound on the entanglement of formation of the two-qubit reduced state. 



9. Conclusion 

As we have seen, many methods for the detection of entanglement have been proposed, ranging from 
Bell inequalities to spin squeezing inequalities. Each of them has its advantages and disadvantages and the 
question which of them is preferable depends on the given experimental implementation. But one has always 
to be careful that the used entanglement detection scheme really detects the desired type of entanglement. 
Furthermore, possible hidden assumptions in an entanglement verification procedure have to be taken into 
account. 

For the future, it can be foreseen that due to the progress in the experimental techniques more and more 
qubits can be entangled. Especially in ion traps and in photonic experiments using hyper-entanglement 
significant progress can be expected. In these experiments the creation of genuine multipartite entanglement 
will probably remain the main goal, as this ensures that the experiment presents something qualitatively 
new compared to previous ones. However, it will become more and more important to design entanglement 
verification schemes that require less experimental effort, as measurements become more demanding. 

On the other hand, as such systems become larger, it can be expected that in many of these experiments 
the qubits will not be individually accessible anymore. Thus, entanglement detection schemes based only on 
collective measurements will be needed. In such experiments, the detection of genuine multipartite entan- 
glement of the total state is not in all cases realistic. However, as we have seen, one can aim to investigate 
how many qubits are entangled, or to quantify the entanglement in the state. 

In this way, quantum information science and in particular entanglement theory can play a crucial role in 
the technological development of quantum control and quantum engineering. Finally, quantum control over 
larger and larger systems might also help to answer fundamental questions concerning quantum theory, such 
as the appearance of a classical macroscopic world based on a quantum microworld. 
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